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COMPUTER PROGRAM FOR QUASI -ONE-DIMENSIONAL COMPRESSIBLE FLOW 
WITH AREA CHANGE AND FRICTION - APPLICATION TO GAS FILM SEALS 


by John Zuk and Patricia J. Smith 
Lewis Research Center 

SUMMARY 

The computer program, AREAX, presented in this report calculates the properties 
of compressible fluid flow with friction and area change. The program carries out a 
quasi -one -dimensional flow analysis which is valid for laminar and turbulent flows under 
both subsonic and choked flow conditions. The program was written to be applied to gas 
film seals. The area -change analysis should prove useful for choked flow conditions with 
a small mean film thickness, as well as for face seals where radial area change is sig- 
nificant. 

The program requires the seal geometry, the reservoir conditions, and the gas 
properties as input. It will then calculate seal mass leakage and volume leakage; force; 
center of pressure; and distributions of pressure, temperature, density, mean friction 
factor, Reynolds number, friction parameter, velocity, and Mach number across the 
seal for both laminar and turbulent flow regimes and for choked and subsonic flow. 

The program is written in FORTRAN IV. Input, internal calculations, and output 
can be in either the International System of Units (SI) or the U. S. customary system. 

This report includes a description of the program, a summary of the mathematical 
model on which the program is based, a complete description of the input variables, and 
sample input and output. 


INTRODUCTION 

Shaft seals in advanced rotating machinery, such as in advanced aircraft turbine en- 
gines, will operate at higher speeds, temperatures, and pressures than shaft seals in 
current use. Because of these severe operating conditions, a positive face separation 
(no rubbing contact) is required for long life and reliability (ref. 1). However, seal- 
face deformation is very likely to occur. 



These deformations may be caused by various distortions (thermal, centrifugal, 
pressure, etc.)- Seal-face distortions become more pronounced under severe operating 
conditions and are usually detrimental to seal performance. Hence, prediction of these 
face deformation effects on gas-film-seal performance is of paramount importance. 

This report presents an analysis and computer program to enable prediction of gas- 
film-face-seal performance when face deformations and/or radial area change is signif- 
icant. The analysis is especially useful for choked flow conditions. 

The computer program presented in this report, AREAX, calculates the properties 
of compressible fluid flow with friction and area change across shaft face seals. It ex- 
tends the capabilities of the program QUASC (ref. 2) to account for the change in flow 
properties resulting from radial area change and face deformations represented by 
linear distortions of the seal faces. The mathematical model for this analysis is given 
in the next section. The analysis includes viscous friction, fluid inertia, area change, 
and entrance losses. Rotational effects and body forces are neglected. The model is 
valid for subsonic flow and for choked flow. It is also valid for both laminar and turbu- 
lent flow regimes. 

Computer programs have proven useful in seal design, where much of the physical 
information of interest is difficult to determine experimentally. The program QUASC 
has proven to be a good model for most applications. It rapidly evaluates physical quan- 
tities of interest. However, QUASC is valid only for parallel sealing surfaces and 
constant -area flow. The program presented here should be used when the effects of 
seal -face distortions are desired and when the radial area change is significant. 

Some of the physical parameters of interest in designing a seal are the leakage, the 
pressure distribution across the seal, and the opening (separating) force. These and 
other parameters are determined by the program AREAX for specified seal geometries 
(flow length, gap, surface deformation), reservoir pressures and temperatures, and gas 
properties. The program also requires two additional parameters. One of these is the 
variation of a mean Fanning friction factor with Reynolds number. The other is the en- 
trance loss coefficient. The program input and output format are almost identical to 
QUASC (except for the accounting of the face deformation and radial area change) . 

This report is intended to serve three purposes: (1) to give a summary of the quasi - 
one -dimensional analysis of compressible flow with friction and area change, (2) to give 
a detailed description of the computer program AREAX which performs this analysis, 
and (3) to serve as a user’s guide for AREAX. 

QUASI -ONE -DIMENSIONAL FLOW ANALYSIS 

The program AREAX is based on a mathematical model for sealing surfaces sepa- 
rated by a narrow gap. One surface may be tilted with respect to the other. The model 
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is valid for both rectangular strip surfaces and coaxial ring surfaces. A pressure dif- 
ference exists across the sealing dam. The cavities on either side of the sealing dam 
are assumed to be constant -pressure reservoirs. 

The flow is assumed to be quasi -one -dimensional. It has been shown that the ef- 
fects of rotation for a circular ring seal geometry are negligible for most gas film seal- 
ing applications (ref. 3). It is assumed that the rectangular surfaces do not slide with 
respect to one another. 

The flow analysis can be separated into two parts. The first is the flow in the pas- 
sage itself, where the flow is assumed to behave as adiabatic flow with friction and area 
change. The second part is flow in the entrance region. These two parts are discussed 
separately . 


Seal Leakage Passage Flow 

It is assumed that the flow in the seal leakage flow region behaves as a variable- 
area adiabatic flow with friction. A qua si -one -dimensional approximation is made 
wherein it is assumed that the flow properties can be described in terms of their cross- 
j sectional averages. 

The following assumptions have been made in the analysis: (1) the effects of rotation 
| are negligible; (2) the flow is adiabatic; (3) no shaft work is done on or by the system; 

(4) no potential energy gradient is present, such as caused by evaluation difference, etc; 

(5) the fluid behaves as a perfect gas. In addition to these assumptions, a simplification 
( to the area change is made in the momentum equation where P dA and A dP » dA dP, 

! or A » dA. This assumption should be satisfactory for most sealing applications. 

The control volume is shown in figure 1. The governing equations with area changes 
reduce to the following differential forms (all symbols are defined in appendix A) : 

Conservation of mass: 

d£ + I^ + dA = 0 (1) 

p 2 u 2 A 

Conservation of energy: 


dT 

T 




( 2 ) 
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Equation of state: 


™ = + (3) 

P p T 

Conservation of momentum (for a small area change): 

-A dP - t w dA w = M du (4) 

With the introduction of the mean Fanning friction factor f and the hydraulic 
diameter D 

dP 

f = ^ = f(x) for radial flow 

2pu 2 
D 

4A 

D = — = 2h for radial flow between parallel disks and plates 
w 

equations (1) to (4) can be combined into a single equation (5) 
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Equation (5) will be referred to as the Mach number equation. This equation is identi- 
cal to the equation obtained from the Table of Influence Coefficients for generalized one- 
dimensional flow in references 4 and 5. 

The other dependent variables can be found in a similar way and are 

du _ -1 dA + yM 2 

u 1 - M 2 A 2(1 - M 2 ) 


4f dx 


( 6 ) 
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(7) 


dT _ (y - 1)M 2 dA y(y - 1)M 4 4fdx 
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Solving equation (5) for 4f dx/D and substituting that into equations (6) to (9) give 


du 
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2M 


dM 

1 +7L_^1m 2 


dT = -fy 
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P A 
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(10) 


( 11 ) 


( 12 ) 


(13) 


Equations (10) to (13) can be integrated directly from the point of interest (M ,A ) 

x X A 

to the point of choking (M = 1,A ) since the variables are separable . The integration 
gives the following equations: 
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(14) 
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( 15 ) 
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(16) 


(17) 


Hence, once the Mach number distribution is known, the physical quantities of interest 
can be found . 

With the introduction of the definition of the flow area A, the radius r, and the film 
thickness h 


A = 2nrh for the coaxial ring geometry 


A = Wh for the rectangular plate geometry 


(18) 


r = Rj ± x 


h = hj + x sin a 


(19) 

(20) 


equation (5) can be written as 


dM* 


2M 


* F(x)sln d! 


to h(l-M 2 ) I F(x) 


(21) 


where 
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x for flow with no radial expansion 


F(x) = 1 + x for radial outward flow 




Rg - x for radial inward flow 


and x = 0 is defined as the flow entrance whether the flow is radially inward or outward. 
The mean friction factor varies with Reynolds number according to the relation (ref. 6) 

f = — (22) 

Re 11 

Consequently, f is an implicit function of x. Equation (21) can be integrated numeri- 
cally from the point of interest (x,M) to the point of choking (x* , M* = 1). 

Other equations that are needed are those for Reynolds number and Mach number : 

Re = 2^- (23) 

V 

M = (24) 

c ^/yStl 


Entrance Flow 


The entrance flow region can be assumed to be an isentropic and adiabatic expan- 
sion (ref. 6). Since real entrance flows are not isentropic because of viscous friction, 
turning losses, and so forth, the entrance loss is accounted for by introducing an em- 
pirically determined entrance velocity loss coefficient C L into the isentropic equations. 
The resulting entrance equations for temperature and pressure become 


and 



(25) 
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When = 1, the flow is isentropic. When is large (e.g. , 30), the entrance 
losses become negligible. 


SOLUTION OF EQUATIONS 

The equations developed in the preceding sections, plus the equation of state written 
as 


— = 5PT (27) 

P 

must now be solved. The independent variable is chosen to be x, the distance across 
the seal face (fig. 2), with the x origin always located at the flow entrance. The known 
parameters are 

Pq sealed-gas pressure (upstream reservoir pressure) 

Pg ambient pressure (downstream reservoir pressure) 

Tq sealed-gas temperature (upstream reservoir temperature, also the stagnation 
temperature) 

entrance velocity loss coefficient 
Ar distance across face of sealing dam 

h(x) film thickness distribution 

k,n constants in friction factor - Reynolds number relation 

These two constants (k,n) apply to laminar flow (Re < Re^) and turbulent flow (Re ^ Re t ), 
respectively. For flow in the transition region (Re^ < Re < Rep, variation of the fric- 
tion factor with Reynolds number is determined by an interpolation method which is de- 
scribed in appendix B of reference 2 . 

There are two types of flow that must be considered: (1) choked flow, when the exit 
Mach number is 1 and the exit pressure is greater than the ambient pressure; and 
(2) subsonic flow, when the exit Mach number is less than 1 and the exit pressure equals 
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the ambient pressure. The method of solution is the same for both types of flow. How- 
ever, the boundary conditions are different. (The lengths and stations used in the analy- 
sis are illustrated for the subsonic -flow case in fig. 3.) 

In both cases, equation (21), the Mach number equation which relates Mach number, 
flow distance, and friction is a differential equation which must be solved numerically. 

The fourth-order Runge-Kutta method is used. Since equation (21) is a first-order dif- 

2 2 

ferential equation in x and M , one boundary condition on M is needed. The only 

. o 2 

known condition is at x = x* , M =1. Equation (21) is solved by guessing a value Mf' 

2 2 jjc l 

for M at x = 0 and integrating until M =1. The final value of x is x . Since the 
Mach number is known to vary rapidly as x approaches x* , the interval of integration 
is divided into small subintervals in that region. 

For choked flow, the solution of the equations for this case is as follows: a value of 

2 * * 

M. is guessed, the Mach number equation is solved, and x is noted. If x = AR, 

1 * / 2 
the solution is complete. If x ? AR, new values of M£ are guessed until the condition 

x* = AR is satisfied. 

For subsonic flow, the solution of the equations for this case is as follows: a value 

2 * * 
of M. is guessed, the Mach number equation is solved, and x is noted. If x is 

1 2 

less than AR, new values of M^ are guessed until the solution of the Mach number 
gives a value of x* greater than AR, When a satisfactory x* is found, Pg is calcu- 
lated. If Pg and Pg are equal, the solution is complete. If Pg and Pg are not 
equal, new values of M^ are guessed until the condition P 2 = Pg is satisfied. 


COMPUTER PROGRAM 

The program AREAX which performs the analysis of quasi -one -dimensional flow 
with friction and area change across shaft face seals is written in FORTRAN IV for the 
IBM 7094H/7044 direct couple computer at the Lewis Research Center. 

The program must be supplied with the geometry of the seal, the gas properties, 
the reservoir conditions, the constants for determining the variation of mean friction 
factor with Reynolds number, and certain logical variables which control output. Input 
and output can be in either the International System (SI) or the U.S. customary system 
of units. 

In general, AREAX performs the following operations in analyzing the flow across 
a seal: It reads the input data and checks that these data are consistent. For instance, 
there are three input variables which can determine the flow length. Since only two are 
necessary, the third must be made consistent with the other two. When the input data 
have been read, AREAX analyzes the flow for each combination of film thickness and 
tilt angle . 


9 



The program first solves the Mach number equation and determines the Mach num- 
ber distribution across the seal face. AREAX then determines the distributions across 
the seal face of pressure; temperature; density; velocity; mean friction factor; 

Reynolds number; mass and volume flow rates; Knudsen number; seal opening force; 
center of pressure; and where appropriate, rotational Reynolds number, variables asso- 
ciated with power dissipation, and axial film stiffness. 

When these data have been calculated, the Mach number distribution across the seal 
face is punched on cards. Since the running time of the program can be fairly long, it is 
convenient to restart it and not have to rerun the cases that are complete. When all the 
data for all the film thicknesses have been calculated, AREAX writes the output data with 
appropriate labels and headings. 

To increase program efficiency and to facilitate program writing, AREAX includes 
a number of subprograms. Figure 4 shows the heirarchy of the subprogram calls. 
Variables are transmitted between the main program and the subprograms through 
labeled COMMON storage. A more detailed description of AREAX and descriptions of 
the subprograms are given in appendix B. Since the programs AREAX and QUASC are 
so similar, the numerical methods described in reference 2 apply to both programs. 

The formulas for the parameters calculated by AREAX are listed in table I. When 
the flow is in either the transitional or turbulent flow regime, the power loss due to ro- 
tation is not calculated. Also, there is no Reynolds number criterion for determining 
turbulence due to rotation. A complete list of the variable names used by the program 
is given in the program listing in appendix C. Flow charts of AREAX and the subpro- 
grams are also in appendix C. 


Input Data 

Input to AREAX is by punched cards. The NAMELIST feature of FORTRAN IV is 
used to read the data. This feature allows the individual variables to be named and 
eliminates complicated card formats. 

The first card required by AREAX is a title card. The title identifies the data and 
uses columns 1 to 72 of one card. It is read by format (12A6). 

The next cards required by AREAX contain the parameters in NAME LIST /SEAL/. 
These parameters define the seal geometry, the gas properties, and the logical varia- 
bles. The variables in /SEAL/ are listed in table II. The next cards required contain 
the parameters in NA ME LIST /HDA T A / which define the gap. These parameters are 
listed in table III. The last cards required contain the parameters in NAMELIST/ 
RESDAT/ which define the reservoir properties. These variables are listed in table IV. 
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Output 


Computer output consists of the input data and calculated parameters. If input is in 
SI units, output is also in SI units. If input is in U. S. units, output is also in U.S. units. 
The printed output parameters are identified, and the units of each are also printed. A 
sample of the output data appears in appendix D. 

The first page of the output contains the title which identifies the data; the input data 
as they are read from cards; the checked input data; the calculated parameters - mean 
flow width and gas constant; and a list of what optional parameters will be calculated, 
namely power, normalized center of pressure, pressure profile factor, and distributions 
across the seal face. The key at the bottom of the page identifies the flow regime asso- 
ciated with a particular film thickness. The key is as follows: 

/ choked flow 

+ flow in laminar -turbulent transitional regime 

T flow in turbulent regime 

The second page contains lists of parameters that vary only with mean film thick- 

ness. These parameters are mass flow rate, standard volume flow rate, mean Knudsen 
number, mean free path of tne gas molecules, axial film stiffness, sealing dam force, 
center of pressure, normalized center of pressure, pressure profile factor, rotational 
flow Reynolds number, choking pressure, choking distance, power, heat generation due 
to viscous shearing, apparent temperature rise due to power dissipation, and torque. 
Power and parameters based on power dissipation are not printed for the transitional or 
turbulent flow regimes since they are not calculated. 

The third and following pages contain lists of parameters that vary with radial dis- 
tance as well as with mean film thickness. These parameters are film thickness, Mach 
number, velocity, density, pressure, temperature, Reynolds number, and mean friction 
factor. The maximum and minimum film thicknesses, the mean film thickness, the rela- 
tive area change, the area change with respect to radial distance, and the tilt angle for 
each film thickness are printed above each set of lists. 


SAMPLE PROBLEM 

An example of the use of the computer program is given here with the following 

2 

conditions: Air at 45 N/cm abs (65.0 psia) is to be sealed from ambient air at 
2 

10.3 N/cm abs (15.0 psia) by a coaxial ring seal operating in the externally pressur- 
ized mode. The mean temperature of the pressurized air is 311 K (100° F). The seal- 
ing dam outside diameter is 16.84 centimeters (6.630 in.), and the inside diameter is 
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16. 58 centimeters (6. 530 in.). The seal faces are separated by a negative linear tilt 
of 1 milliradian. (The smallest gap is located at the outer radius.) The design surface 
speed is 61 meters per second (200 ft/sec). 

It is desired to find a design film thickness which is large enough so that power dis- 
sipation and viscous heating temperature rise are sufficiently low, yet small enough so 
that the mass leakage is tolerable. From our experience, the best method is to try 
mean film thickness inputs of 7.62 to 25. 4 micrometers (0. 3 to 1.0 mil) in increments 
of 2. 54 micrometers (0. 1 mil) . However, to give a sample output for transitional flow 
and turbulent flow, the range of film thicknesses has been increased to 40.64 microm- 


eters (1.6 mils). For this study, isentropic entrance conditions are assumed. Thus, 
the program input will include 

Mean rotational velocity , V, m/sec (ft /sec) 61 (200) 

Molecular weight of gas, m 28.9660 

Reservoir temperature, T n , K (°F) 311 (100) 

Reservoir pressure (highest pressure), Pq, N/cm^ abs (psia) 45 (65.0) 

Ambient pressure (lowest pressure), Pg, N/cm^ abs (psia) 10.3(15.0) 

Specific heat at constant pressure, c , J/kg-K (Btu/lbm-°R) 10^ (0.24) 

r 

Film thickness (increase in increments of 2. 54 pm (0. 1 mil)), 

h, pm (mil) 7.62 to 40.64 (0.3 to 1.6) 

Tilt angle, a, rad -0.001 

Inner radius, R^, cm (in.) 8.300 (3.265) 

Flow length, AR, cm (in.) 0.127 (0.050) 

Specific -heat ratio, y 1.4 

Numerical constant in laminar friction factor - Reynolds number relation, k^ 24 

Exponent in laminar friction factor - Reynolds number relation, n^ 1.00 

Numerical constant in turbulent friction factor - Reynolds number 

relation, k^. 0.079 

Exponent in turbulent friction factor - Reynolds number relation, n^. 0.25 

Loss coefficient, 1-00 


The input data sheet for this sample problem is given in table V, in both SI and U.S. 
units. Plots of profiles of pressure, temperature, density, Mach number, Reynolds 
number, and mean friction factor for a mean film thickness of 12.7 micrometers 
(0. 5 mil) are shown in figure 5. 

REGION OF VALIDITY AND USAGE AND PROGRAM EFFICIENCY 
Results obtained using the AREAX computer program are compared with known 
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solutions to check their validity. The first case considered is a parallel-film seal with 
a relatively large radial area change . 

Figure 6 shows both the pressure distribution for pure radial viscous flow found 
from this variable-area analysis (AREAX) and the analytical solution obtained from the 
classical viscous flow model (ref. 3), which is 



o 

The conditions are representative of aircraft idle operation: P Q = 45 N/cm abs 
(65 psia), Pg = 10.3 N/cm'* abs (15 psia), T Q = 311 K (100° F), = 5.880 centimeters 

(2.315 in.), and R 2 = 8.410 centimeters (3.315 in.) (AA/A = 0.43). The parallel- 
surface case of 12. 7 -micrometer (0.5-mil) film thickness was solved. The variable- 
area analysis shows excellent agreement with the analytical solution. Also shown in fig- 
ure 6 is the pressure profile obtained when a constant friction factor calculated at the 
mean radius is used. This constant-friction-factor approach slightly underestimates 
the pressure along the seal passage length. 

Figure 7 compares the pressure profiles obtained by using the area expansion 
analysis with those from the viscous flow solution for the case of a positive linear tilt of 
1 milliradian. The pressure profile predicted by the viscous compressible flow solution 
is (ref. 7) 


1/2 
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~^\ -i 
ip 
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xh 2 (2h 1 + ax) 


< E 2 ' R l )2h m( h l + ax )‘ 


(29) 


The conditions are the same as for the case shown in figure 6, except that R^ equals 
8.300 centimeters (3.265 in.). Hence, the radial area change is negligible. 

At a mean film thickness of 5.1 micrometers (0.2 mil) and a positive linear tilt of 
1 milliradian, there is excellent agreement between the present analysis and the small- 
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tilt analysis of reference 7. Also shown in figure 7 is the corresponding case of a nega- 
tive linear tilt of 1 milliradian but a mean film thickness of 4. 4 micrometers 
(0.175 mil). Again, excellent agreement is found. 

Figure 8 shows pressure distribution results obtained from the present analysis for 
a divergent tilt of 2 milliradians. Distributions for mean film thicknesses of 2. 5, 5.1, 
7.6, and 12.7 micrometers (0.1, 0.2, 0.3, and 0.5 mil) are presented. The other con- 
ditions were P Q = 148 N/cm 2 abs (215 psia), Pg = 10.3 N/cm 2 abs (15 psia), T Q = 700 K 
(800° F), R 1 = 8. 300 centimeters (3. 265 in. ), and Rg = 8. 410 centimeters (3. 315 in. ). 
These conditions are representative of advanced aircraft cruise conditions (ref. 1). 
These conditions represent subcritical (subsonic), critical (Pg = Pg and Mg = 1), and 
supercritical (choked) flow conditions. Also shown is the parallel -film pressure pro- 
file for 2. 5-micrometer (0. 1-mil) film thickness. This is the classical parabolic pro- 
file for viscous compressible flow. In addition, figure 8 shows a supercritical flow 
pressure profile for parallel sealing surfaces and for a film thickness of 12. 7 microm- 
eters (0. 5 mil) which was obtained using the constant -area analysis (QUASC) of refer- 
ence 2. The present analysis shows excellent agreement with this parallel-film profile 
with a 12. 7-micrometer (0. 5-mil) film thickness. Similar results were found for the 
case of 2 -milliradian convergent tilt. 

Experience has shown that QUASC (ref. 2) performs each film -thickness case six to 
ten times faster than AREAX. Hence, for parallel surfaces, QUASC is more efficient 
and should be used. The small-seal -face-deformation cases using AREAX on the IBM 
7094/7044 direct couple computer required about 0. 16 minute for each film thickness 
where the flow was choked and about 0. 5 minute for each film thickness where the flow 
was subsonic. 

For small face deformations and subsonic flow conditions, the analytical solution 
(ref. 6) should be used. For deformed surfaces with a relatively large mean film thick- 
ness and choked flow, QUASC (ref. 2) may give satisfactory results (faster computer 
solution time) . For a more complete solution, AREAX should prove useful for choked 
flow conditions with a relatively small mean film thickness, as well as for face seals 
where radial area change is significant. 


CONCLUDING REMARKS 

A summary of the quasi -one -dimensional analysis of compressible flow with friction 
and area change has been presented. Also, a detailed description of the computer pro- 
gram AREAX, which performs this analysis, is given. This program has proven useful 
in extending the capabilities of the computer program QUASC by including the area 
change due to (1) change in radius and (2) deformed seal faces represented by a small 
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linear tilt. Results obtained using AREAX showed excellent agreement with analytical 
solutions for pure radial viscous flow and for viscous flow with a small tilt of the sealing 
surfaces. Favorable agreement was also found between AREAX and QUASC (ref. 2) for 
a parallel film under choked flow conditions. An example mainshaft seal problem is also 
given. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 7, 1973, 

501-24. 
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APPENDIX A 


SYMBOLS 

A area 

Cj^ entrance velocity loss coefficient 

c speed of sound 

c specific heat at constant pressure 

P 

D hydraulic diameter (D = 2h for coaxial rings and narrow slots) 

F seal opening force 

F dimensionless seal -opening -force pressure profile factor 

f mean Fanning friction factor 

H shear heat 

h film thickness (gap) 

Kn mean Knudsen number, X/h m 

k numerical constant in friction factor - Reynolds number relation 


M 

M 


Mach number , u 

i 

mass flow 



m molecular weight of gas 

n exponent in friction factor - Reynolds number relation 

P pressure 

wetted perimeter 
Q volume leakage flow rate 

R^ inner radius 

Rg outer radius 

AR sealing dam radial width (physical flow length), Rg - R^ 
£ gas constant (universal gas constant/molecular weight) 

£ u universal gas constant 

Re Reynolds number 

Re(p) Reynolds number of leakage flow 
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Re(r) Reynolds number of rotational flow 
r radius 

S axial film stiffness 

T temperature 

u leakage flow mean velocity (x -direction) 

V mean rotational velocity 

W flow width 

x flow direction coordinate 

x c center of pressure 

x dimensionless center of pressure 

a relative inclination angle of surface (positive a designates hg > h^) 

y specific -heat ratio 

X mean free path of gas molecules 

Id absolute (dynamic) viscosity 

p density 

t shear stress 

Subscripts : 

l laminar flow 

m mean 

max maximum 

t turbulent flow 

w wetted surface area or wall 

x location along flow leakage length from entrance 

0 sealed (reservoir) conditions 

1 flow entrance conditions 

2 flow exit conditions 

3 ambient sump conditions 

Superscript: 

* 


critical flow condition 


appendix b 


DETAILED DESCRIPTION OF PROGRAM 

This appendix contains details regarding checking the consistency of the input data, 
traps for invalid data, and the subprograms. The variables mentioned in this appendix 
are defined and listed either in the main-text section Input Data or in the program list- 
ings (appendix C) . Further details of numerical methods and subroutines that are used 
by both AREAX and QUASC can be found in reference 2. 


Main Program AREAX 

The main program AREAX controls the complete quasi -one -dimensional analysis. 

It defines the labeled COMMON storage for transmission of data among the subprograms. 
It defines some constants and labels for the output. It reads the input data, checks that 
these data are consistent, and rejects cases for which the input data are invalid. It 
calls subprograms to solve the Mach number equation and to determine the distributions 
across the seal face of pressure, temperature, density, velocity, mean friction factor, 
and Reynolds number. It calculates the other parameters associated with the flow, such 
as mass flow rate, and writes the output. AREAX then transfers to read new input data. 

The labeled COMMON storage defined by AREAX contains constants needed by the 
subprograms and the output data. Most of the variables in the COMMON blocks are 
listed in tables VI to IX. The COMMON blocks not described contain variables which 
are useful in the operation of the program but are not necessary for the analysis. The 
variable names and the array names are not listed because the names are not the same 
in all the subprograms. 

There are two kinds of parameters in the COMMON block /ARRAYS/ (table VI): one 
is the distributions across the face of the seal, the other is the integrated flow param- 
eters that vary only with mean film thickness. The arrays are dimensioned for 20 film 
thicknesses and 11 points across the seal face. 

The parameters in COMMON block /CONSTS/ (table VII) are constants needed by the 
program for internal calculations. The array CNVT has dimension (11,2). The first 
column of the array contains constants for calculations in SI units. The second column 
of the array contains constants for calculations in U.S. customary units. The variable 
NSI is 1 for SI units and 2 for U.S. units. Table VIII lists the parameters for which each 
element of CNVT is used. 

The arrays in COMMON block /TRAYS/ (table IX) contain the parameters used in the 
solution of the Mach number equation for the film thickness under consideration. 
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COMMON blocks /CTITLE/ and /PRNT/ contain the title, which identifies the data, and 
the labels for identifying which parameters are calculated. COMMON block /LOGICL/ 
contains logical variables. 

The program AREAX first reads the input data and checks that these data are con- 
sistent. One parameter it checks for consistency is the flow length, RDIF. There are 
three input variables that can be used to determine RDIF. These variables are RINNER, 
ROUTER, and RDIFIN. Since only two are necessary, the third must be made consis- 
tent with the other two. The check is made as follows: If RINNER 4 0 and ROUTER 4 0, 
AREAX sets R1 = RINNER and R2 = ROUTER and then calculates RDIF = R2 - Rl. If 
RINNER 4 o, ROUTER = 0, and RDIFIN 4 0, AREAX sets Rl = RINNER and 
RDIF = RDIFIN and then calculates R2 = Rl + RDIF. If RINNER = 0, ROUTER 4 0 , and 
RDIFIN 4 0, AREAX sets R2 = ROUTER and RDIF = RDIFIN and then calculates 
Rl = R2 - RDIF. These three conditions are for the coaxial ring geometry. The varia- 
ble WIDTH is calculated as 27 tR. If RINNER = 0, ROUTER = 0, and RDIFIN 4 0, the 
geometry is for the rectangular plate. Then AREAX sets RDIF = RDIFIN, and WIDTH 
must be nonzero. Any other combination of RINNER, ROUTER, RDIFIN, and WIDTH is 
considered an error since there is not enough information available to determine a non- 
zero flow length and flow width. If this error condition exists, AREAX writes a mes- 
sage and transfers to read new data from NAMELIST/SEAL/. 

The second parameter that is checked for consistency is the pressure ratio. Three 
input variables are available for determining the pressure ratio, but only two are neces- 
sary. The check is made the same way as the check for flow length. If P1IN 4 0 and 
P2IN 4 0, AREAX sets PI = P1IN and P2 = P2IN and then calculates PRAT = P1/P2. If 
P1IN 4 0, P2IN = 0, and PRIN 4 0, AREAX sets PI = P1IN and PRAT = PRIN and then 
calculates P2 = Pl/PRAT. If PUN = 0, P2IN 4 0, and PRIN 4 0, AREAX sets P2 = P2IN 
and PRAT = PRIN and then calculates PI = PRAT x P2. Any other combination of P1IN, 
P2IN, and PRIN is considered an error. In that case, AREAX writes an error message 
and transfers to read new data according to the input code. Since the flow may be 
radially inward or outward for the coaxial ring geometry. AREAX chooses the larger 
value of PI and P2 to be the sealed-gas pressure P0 and the smaller value to be the am- 
bient pressure P3. 

The third parameter that must be checked for consistency is the rotational velocity . 
If SPEED 4 0, CAPV is calculated from SPEED. If CAPV 4 0 and SPEED = 0, SPEED 
is calculated from CAPV. If both are 0, the system is considered static, and the logical 
variable PWRSKP is set to . TRUE . to omit calculations involving power . 

For each film thickness and tilt angle combination, AREAX calls subroutine 
NCFLOW to solve the Mach number equation for the given film thickness. Subroutine 
NCFLOW also punches the solution of the Mach number equation on cards for restarting 
the program if it rims out of time. AREAX then calls subroutine DISTS to calculate the 
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distributions across the seal face. Function subprogram SIMPS1 is used for the numer- 
ical integrations in the calculation of force and center of pressure. When all the data 
for all the film thicknesses have been calculated, AREAX call subroutine STFNSS to 
determine the axial film stiffness. 

When all the calculations are complete, AREAX writes the input data, the 
"checked" input data, and the output data with appropriate headings and labels. The 
final command in the program is a transfer to read new input data. 


Subprogram NCFLOW 

Subprogram NCFLOW controls the solution of the Mach number equation, for any 
given film thickness and tilt angle combination. The subprogram calls subprogram RK1 
to solve the Mach number equation. NCFLOW iterates first on the boundary condition 
x* > AR. If x* > AR, NCFLOW then iterates on the boundary condition P 2 = Pg. 
When the solution for each film thickness is complete, NCFLOW punches the x distri- 
tion and the Mach number distribution on cards for restarting the program. 


Subprogram DISTS 

Subprogram DISTS determines the distributions across the seal face of pressure, 
temperature, Mach number, density, velocity, Reynolds number, and mean friction 
factor. It uses a three -point Lagrange interpolation on the data from the solution of the 
Mach number equation to get the Mach number distribution at the given x grid points . 

It then calculates the other parameters as follows: temperature from equation (15), 
velocity from the definition of Mach number, density from the continuity equation, pres- 
sure from the equation of state, Reynolds number from the definition of Reynolds num- 
ber, and mean friction factor from the friction factor - Reynolds number relation. 


Subprogram PRESS 

Function subprogram PRESS uses a three -point Lagrange interpolation to determine 
the Mach number and film thickness at distance x across the seal face. Then pressure 
is calculated from equation (17). 
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Subprogram RK1 


Subprogram RK1 solves the Mach number equation by the fourth order Runge-Kutta 
solution. Since the flow must remain subsonic, there are traps built into the subpro- 
gram to ensure that the Mach number remains less than 1. The initial guess for M^ 
is supplied by the calling program. The equation is considered solved when the final 
M is 1. 0. All variables and operations in this subprogram are in double precision. 


Subprogram SIMPS1 

Function subprogram SIMPS1 uses Simpson's rule to evaluate the integrals in the 
force and center -of -pressure equations. Reference 2 provides a detailed description of 
this subprogram. 


Subprogram STFNSS 

Subprogram STFNSS performs a numerical differentiation to determine the axial 
film stiffness. Reference 2 provides a detailed description of the numerical technique 
used. 


Subprogram FRFUNC 

Function subprogram FRFUNC determines the mean friction factor for a given 
Reynolds number. FRFUNC examines the Reynolds number to determine whether the 
flow is laminar, turbulent, or transitional. It then uses the input constants to deter- 
mine f . If the flow is transitional, FRFUNC calculates f by an interpolation formula 
developed in reference 2. 


Subprograms FFUNC and XCFUNC 

Function subprograms FFUNC and XCFUNC evaluate the integrands in the force and 
center -of -pressure integrals. Both FFUNC and XCFUNC call subroutine PRESS, which 
evaluates the pressure at any given x. 
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Subprogram SORTXY 


Subprogram SORTXY rearranges the ordered pairs of numbers (X(l), Y(l)), 

(X(2), Y(2)), . . . , (X(N), Y(N)), which are stored in arrays X and Y, in order of ascend- 
ing X. The (X,Y) pairing is preserved. 


Subprogram GRAFIC 

Subprogram GRAFIC makes the following plots: 

(1) Power against mean film thickness 

(2) Center of pressure against mean film thickness 

(3) Force against mean film thickness 

(4) Pressure against x 

(5) Temperature against x 

(6) Density against x 

(7) Mach number against x 

(8) Reynolds number against x 

(9) Mean friction factor against x 

The subprogram presented here is a dummy routine to satisfy the call to GRAFIC 
from the main program. Each user must write his own plotting subroutine, using the 
dummy subprogram and the flow chart as a guide, that is appropriate to the plotting 
devices available. 


Subprogram BLOCK 

Subprogram BLOCK is a block data subprogram to load constants into labeled 
COMMON. 
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APPENDIX C 


FLOW CHARTS AND PROGRAM LISTING 


AREAX 
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[ Write data identification , 


(Write flow length and reservoir data 



Write consistent input data and array of labels! 


< JDONE 

J > ~ ( 201> 


'Read solution of Mach number equation^ 
.from previous running of program J 




Skip calculations based^-v. No 
^n power dissipation? s' 


Calculate power, shear heat, 
temperature, and torque 


Calculate sealing dam force and center of pressure^ 


24 











































Comment: Runge-Kutta solution of Mach number equation 
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OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


1 

QUASI-CNE DIMENSIONAL COMPRESSIBLE FLOW WITH FRICTICN AND AREA 2 

CHANGE. AREA CHANGE MAY BE DUE TO RADIAL EXPANSION AND/OR 3 

RELATIVE TILTING OF THE SEALING SURFACES. A 

5 

INPUT VARIABLES 6 

*++**♦+*$****** 7 

8 

TITLE - ALPHANUMERIC IDENTIFICATION OF THE DATA 9 

IC 

NS I - NUMERICAL SIGNAL WHICH INDICATES WHICH SYSTEM CF UMTS IS 11 

USED FOR INPUT, OUTPUT, AND INTERNAL CALCULATIONS 12 

(NSI = 1 MEANS SI UNITS - NS I = 2 MEANS US UNITS I 13 

RINNER - INNER RADIUS OF CIRCULAR SEAL 14 

ROUTER - CUTER RADIUS OF CIRCULAR SEAL 15 

RDIFIN - DISTANCE ACROSS SEAL 16 

WIDTH - FLOW WIDTH FOR SEALS WITH NO RADIAL EXPANSION 17 

MOLWT - MOLECULAR WEIGHT 18 

CP - SPECIFIC HEAT OF GAS IS 

MUIN - RESERVOIR VISCOSITY FOR GAS OTHER THAN AIR 20 

GAMMA - RATIO OF SPECIFIC HEATS 21 

SPEED - ROTATIONAL VELOCITY 22 

CAP V - SURFACE SPEED 23 

XL AM - EXPONENT IN FORMULA FOR VARIATION OF FRICTION FACTOR 24 

WITH REYNOLDS NUMBER FOR LAMINAR FLCW 25 

XTURB - EXPONENT IN FORMULA FOR VARIATION OF FRICTION FACT CR 26 

WITH REYNOLDS NUMBER FOR TURBULENT FLOW 27 

CONLAM - CONSTANT IN FORMULA FOR VARIATION OF FRICTION FACTOR 28 

WITH REYNOLDS NUMBER FOR LAMINAR FLOW 29 

CONTRB - CONSTANT IN FORMULA FOR VARIATION OF FRICTION FACTOR 30 

WITH REYNOLDS NUMBER FOR TURBULENT FLOW 31 

RELAM - UPPER LIMIT OF REYNOLDS NUMBER FOR LAMINAR FLCW 32 

RETURB - LOWER LIMIT OF REYNOLDS NUMBER FOR TURBULENT FLOW 33 

PWRSKP - LOGICAL VARIABLE - IF TRUE, SKIP CALCULATIONS 34 

INVOLVING POWER 35 

PRSSKP - LOGICAL VARIABLE - IF TRUE, SKIP PRINTOUT CF 36 

DISTRIBUTIONS OF PRESSURE, TEMPERATURE, DENSITY, 37 

MACH NUMBER, VELOCITY, AND FRICTION PARAMETER 3e 

ACROSS THE FACE OF THE SEAL 39 

NRMSKP - LOGICAL VARIABLE - IF TRUE, SKIP NORMALIZATION OF 4C 

FORCE AND CENTER OF PRESSURE 41 

PLTSKP - LOGICAL ARRAY - IF ANY ELEMENT IS TRUE, SKIP THE 42 

PLOT OF THE CORRESPONDING VARIABLE 43 

PLTSKP(l) - POWER VS CHARACTERISTIC FILM THICKNESS 44 

P L T SK P ( 2 ) - CENTER OF PRESSURE VS CHARACTERISTIC 45 

FILM THICKNESS 46 

PLTSKPI3) - FORCE VS CHARACTERISTIC FILM THICKNESS 47 

PLTSKPI 4) - PRESSURE VS X 48 

P LT SKP ( 5) - TEMPERATURE VS X 49 

PLTSKPI 6) - DENSITY VS X 50 

P L T SKP ( 7) - MACH NUMBER VS X 5 1 
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c 

c 

r 

C 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PL7 SKP( 8) - REYNOLDS NUMBER VS X 
PLTSKP< 9) - MEAN FRICTION FACTOR VS X 


HMtAN - MEAN FILM THICKNESS 
ALPHA - TILT ANGLE 

NJ - NUMBER OF FILM THICKNESS 
JDONE - NUMBER OF CASES COMPLETED 
OF THE PROGRAM 


IN PREVIOUS RUNNING 


C EXIT OF SEAL WITH NO RADIAL EXPANSION 

C PRIN - PO/P3 

C TOIN - SEALED GAS TEMPERATURE 

C LOSS - ENTRANCE VELOCITY LOSS COEFFICIENT 

C INCODE - INPUT CODE 


H - FILM THICKNESS THAT CORRESPONDS TO X 
PARAMETERS THAT VARY ONLY WITH FILM THICKNESS 

44444444444444444444444444444444444444444444 * 

MOOT - MASS FLOW RATE 
Q - VOLUME FLOW RATE 
KN - KNUDSEN NUMBER 

LAMCA - MEAN FREE PATH OF GAS MOLECULES 
RER - ROTATIONAL FLOW REYNOLDS NUM8ER 
POWER - POWER DISSIPATED DUE TO VISCOUS SHEARING 
DELT - APPARENT TEMPERATURE RISE DUE TO VISCOUS SHEARING 
H SHEAR - SHEAR HEAT 
TORQUE - TORQUE 
FORCE - SEAL ING DAM FORCE 
F PAR - NORMALIZED SEALING DAM FORCE 
XC - CENTER OF PRESSURE 
XCBAR - NORMALIZED CENTER OF PRESSURE 
STIFF - AXIAL FILM STIFFNESS 

HI - FILM THICKNESS AT INNER RADIUS CF SEAL 
H2 - FILM THICKNESS AT OUTER RADIUS OF SEAL 

PARAMETERS THAT VARY WITH FLOW LENGTH AND FILM THICKNESS 

444444 4444444444 44444444444444444 444444444 44 44 444 4 44 4 4 44 

P - PRESSURE 
MACH - MACH NUMBER 

U - RADIAL VELOCITY 


52 

53 

54 

55 

56 

57 
56 
59 
6C 


C 

P1IN- PRESSURE 

AT 

INNER RADI 

US 

OF 

CIRCULAR 

SEAL 

CR 

PRESSURE 

AT 

61 

c 

ENTRANCE 

OF 

SEAL WITH 

NO 

RADIAL EXPANSION 




62 

c 

P2IN - PRESSURE 

AT 

OUTER RADI 

US 

OF 

CIRCULAR 

SEAL 

CR 

PRESSURE 

AT 

63 


64 

65 

66 

67 

68 


C 

I NCODE = 1» 

PROGRAM 

TRANSFERS 

TC 

READ 

NEW 

TITLE CARD 

69 

C 

INCODE = 2» 

PROGRAM 

TRANSFERS 

TO 

RE AO 

NEW 

SEAL DATA 

7 C 

c 

I NCODE=3* 

PROGRAM 

TRANSFERS 

TC 

READ 

NEW 

FILM THICKNESS 

71 

c 


DATA 






72 

c 

INC0DE = 4, 

PROGRAM 

TRANSFERS 

TO 

READ 

NEW 

KESERVCIR 

73 

c 


DATA 






74 

c 








75 

c 

OUTPUT VAR IABLES 







76 

c 

444444 4 444444 444 







77 

c 








78 

c 

X - DISTANCE ACROSS SEAL 

FACE 





79 


80 

81 

82 

83 

84 

85 

86 

87 

88 
89 
9C 

91 

92 

93 

94 

95 

96 

97 
9e 
99 

IOC 
10 1 
102 

103 

104 
1C 5 
IC6 
1C7 
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ooo 000000000000000000000000000000 


T - TEMPERATURE ICE 

RHO - DENSITY ICS 

REP - PRESSURE FLOW REYNOLDS NUMBER 1IC 

FRICT - FRICTION FACTOR 111 

112 

PROGRAM VARIABLES - SCME OF THESE ARE OUTPUT, ALSO 113 

************************************************** 1 1 A 

115 

RC - RADIUS OF SEAL AT INLET (RO = R1 FOR INTERNALLY lit 

PRESSURIZED CIRCULAR SEALS, RO = R2 FOR EXTERNALLY 117 

PRESSURIZED CIRCULAR SEALS, AND RO = 0 FOR SEALS WITH NO 118 

RADIAL EXPANSION) US 

R GAS - GAS CONSTANT 12C 

FAREA - SURFACE AREA 121 

MU - RESERVOIR VISCOSITY 122 

TO - SEALED GAS TEMPERATURE 123 

PTOL - TOLERANCE FOR CUT-CFF ON PRESSURE ITERATION 124 

HMIN - MINIMUM FILM THICKNESS 125 

CNVT - ARRAY OF CONSTANTS NEEDED IN THE CALCULATIONS 126 

( CM VT ( I ,11 ARE CONSTANTS IN SI UNITS) 127 

(CNVT (1,2) ARE CONSTANTS IN US UNITS) 128 

129 

PARAMETEPS CALCULATED AT POINT OF CHOKING 13C 

***************************************** 131 

132 

XSTAR - DISTANCE FROM ENTRANCE TO POINT OF CHOKING 133 

PSTAR - PRESSURE 134 

ASTAR - CROSS-SECTIONAL AREA 135 

TSTAR - TEMPERATURE 136 

137 

REAL MACH, MOLWT, LOSS, KN , L AMOA , MUI N , MU, MOOT 138 

DOUBLE PRECISION XX,FM,HH 139 

LOGICAL PWRSKP,PRSSKP,NRMSKP,PLTSKP 140 

EXTERNAL FFUNC.XCFUNC 141 

DIMENS ION CALC ( 8, 4),M00T( 20) ,Q(20) ,KN(20) ,LAMDA (20 I ,RER(20 ), 142 

1 XC ( 20 ) , DELT( 20) ,HSHEAR(20) , TORQUE (20) ,HMIN(20 ) 143 

COMMON /ARR A YS/X(11),P( 1 1 , 20 * , MACH ( L l , 20 > ,U( 11 ,20) ,T (LI ,20) , 144 

1 PFO( 11, 20), REP( 11,20) ,FRICT( 11,20) ,H(11 ,20) ,XCBAR(20 ) , 145 

2 POWER ( 20 ) , FORCE (20) ,STIFF(20) ,FBAR(20),H1(20) ,H2(20) , 146 

1 ALP HA ( 20 ) ,HME AN (20) ,XSTAR<20) , PSTAR (20) ,ASTAR(20) 147 

COMMON /CONSTS /GAM MA,RD IF ,R0, SIGN, XL AM, XTURB ,C CNLAM ,CONT RB, RELAM, 148 

1 RETURB, T0,P0,P3,PT0L,RGAS,L0SS,MU,PI ,RUNI V(2) , CNVT (11, 2), NS I 149 

COMMON/TRAYS/N,XX( 201 ),FM( 201) ,HH( 201) ,J, TSTAR 15 C 

C0MM0N/CTITLE/TITLE( 12) 151 

COMMON /PRNT/C 1( 4) ,C2(4),C3(4),C4(4) ,C5(4I,C6(4) ,C7(4),C8(4), BLANK 15 2 

COMMON /LOGICL/PWRSKP, NRMSKP ,PR S SKP , PLTSKP (9 ) 153 

NAM EL IST/SEAL/RINNER,ROUTER,RDIFIN, WI DTH, MOLWT , CP, MUIN, GAMMA, CAPV, 15 4 

1 SPEED, XLAM, XTURB, CONLAM,CONTRB, RELAM, RETURB, PWRSKP, PRSSKP, 155 

2 NRMSKP, PLTSKP ,NSI 156 

NAMEL IST/HDATA/ HMEAN , ALPHA ,NJ, JDONE 157 

NAM EL IST/RESDAT/P 1IN,P2IN,PRIN, TCI N, LOS S , I NCO DE 15 e 

159 

READ TITLE CARD AND SEAL DATA 16C 

161 

100 R EAD( 5, L> TITLE 162 

110 R EAC( 5, SEAL ) 163 
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164 

CHECK CCNSISTENCY OF FLOW LENGTH PARAMETERS 165 

166 

IF (R INNER .EQ .0.0) GO TO 111 167 

IF (ROUTER. EQ. 0.0) GO TO 112 168 

R 1 = R INNER 165 

R 2 = RCLTER 17C 

ROIF = R 2-R 1 171 

GO TO 116 172 

17 3 

111 IF (ROUTER. EQ. 0.0) GO TO 113 174 

IF (RDIFIN.EQ.0.0) GO TO 114 175 

RDIF = RO IF IN 176 

R2 = ROUTER 177 

R 1 = R2-RDIF 176 

GO TO 116 17<3 

18C 

112 IF (RDIFIN.EQ.0.0) GO TO 114 181 

RDIF = RDIFIN 182 

R 1 = RINNER 183 

R 2 = Rl+RDIF 184 

GO TO 116 185 

186 

113 IF (RDIFIN.EQ.0.0) GO TO 114 187 

RDIF = RDIFIN 188 

R1 = 0.0 185 

R 2 = O.C 1<3C 

IF (W ICTH.NE.O.C) GO TO 116 191 

192 

ERROR IN FLOW LENGTH DATA 192 

194 

114 WRITE (6,10) TITLE, RINNER, ROUTER, RDIFIN, WIDTH 195 

115 R EAD ( 5.RESDAT) 196 

GO TO ( 100, 110, 115, 115) , INCODE 197 

19 8 

FLOW LENGTH DATA CONSISTANT 195 

CHECK CONSISTENCY OF PLOT CONTROL PARAMETERS 20C 

20 1 

116 IF (NRMSKP) PLT SKP( 3 » = . TR IE . 202 

IF ( .NCT.PRSSKP) GO TO 118 203 

DO 117 1=4,10 204 

117 PLTSKP(I) = .TRUE. 205 

118 DO 119 1=1,11 206 

119 X(I) = FLOAT( I- l)*RDI F/10.0 207 

2 C 8 

READ FILM THICKNESS DATA 205 

2 1 C 

120 R EAC( 5 , HCATA ) 211 

212 

READ RESERVOIR DATA ANO CHECK CONSISTENCY OF PRESSURE DATA 213 

214 

130 REAC( 5 ,RESDAT ) 215 

IF (P 1IN.EQ.0.0) GO TO 134 216 

IF (P2IN.E0.0.0) GO TO 132 217 

131 PI = P1IN 218 

P 2 = P 2 IN 219 
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PR = P 1/P2 220 

GO TO 140 221 

C 222 

132 IF (PR IN.EQ.O .0 > GO TO 133 223 

PR = PRIN 224 

PI = PI IN 22* 

P 2 = P 1/PR 226 

GO TO 140 227 

C 228 

133 PI = P1IN 229 

P 2 = P2 IN 230 

PR * C. 231 

GO TO 140 232 

C 233 

134 IF (P21N.NE.0.0) GO TO 131 234 

235 

ERROR IN PRESSURE OATA 236 

237 

135 WRITEU, 12) TITLE 238 

GOTO ( 100, 110, 12C, 1301, INCOOE 239 

24C 

PRESSURE CATA CONSISTANT 241 

242 

140 PO * AMAX 1 ( P l , P 2 ) 243 

P 3 = AMIN1(P1,P2) 244 

245 

WRITE INPUT OATA 246 

247 

WRITE (6,20) TITLE 248 

IF (NSI.EQ.1I WRITE (6,22) RI NNER , P1I N , MOLWT , RCUT ER , P2I N ,CP, 249 

1 CONLAM, RDIFIN, PRIN, GAMMA, XLAM, WIDTH, TOI N,MUlN,RtLAM,LCSS, 250 

2 SPEEO,CONTRB,CAPV,XTURB,RETURB 251 

IF (NSI.EQ.2) WRITE (6,241 RI NNER , P l I N , MOLWT , ROUT ER , P2IN , CP, 252 

1 CONLAM.RDIF IN, PRIN, GAMMA, XLAM, WIDTH, TOI N.MUIN , RELAM, LCSS, 253 

2 SPEED, CONTRB, CAPV, XTURB ,RETURB 254 

255 

CALCULATE PROGRAM CONSTANTS 256 

257 

RGAS = RUNIV< NSII/MOLWT 258 

SIGN = 1.0 259 

IF (WICTH.NE.O.O) GO TO 141 260 

W IDTF = P I*(R 1*R2) 261 

FAREA = PI*(R2**2-R1**2) 262 

IF (P2.GT.P1) SIGN = -1.0 263 

RO = R 1 264 

IF (P2.GT.P1) RO=R2 265 

GO TO 142 266 

14 1 RO x C .0 267 

FAREA * W I0TH4RDI F 268 

CAPV *0.0 269 

SPEED * 0.0 27C 

142 TO = T C IN 271 

IF (NSI.EQ.2) TO=TO I N+460 • 272 

MU * ML IN 273 

IF (GAMMA. EQ. 1.4) MU*CNVT ( 1 , N SI )*TO**1.5/ (TO+CNVT (2 , NS I ) ) 274 

N = AL OG 10( P 3 ) 275 
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PTOL = 5.*10.**(N-4> 276 

IF (SPEED. EQ.O.C) GO TO 143 277 

CAP V = PI*SPEED*(R 1+R2)/CNVT(3,NSI ) 278 

GO TO 145 279 

143 IF (CAPV.EQ.O .0 > GO TO 144 28C 

SPEED = CNVT( 3, NSI )*CAPV/PI/(R1«-R2) 281 

GO TO 145 282 

144 PWRSKP = .TRUE. 283 

CP = 0.0 284 

145 IF (PWRSKP I PLTSKPI 1 )=.TRLE. 285 

286 

CALCULATE ENTRANCE, EXIT, AND MINIMUM FILM THICKNESSES - ALSO 287 

CALCULATE RADIAL DISTRIBUTION OF FILM THICKNESS 288 

289 

DO 147 J=1,NJ 290 

Hl(J) * HMEAN ( J 1- . 50*RDIF*S I N(ALPHA(J!) 291 

H2( J ) = 2 .0* HMEAN < J l-H 1( J 1 292 

HM INI J 1 = AMINKHK J) ,H2( J) > 293 

00 146 1*1,11 294 

146 H( I ,J ) = H1(J )+( R0-R1 +X ( I )* SIGN) *SI M AL PHA ( J) I 295 

147 CONTINUE 296 

297 

SET UP ARRAY OF LABELS 298 

299 

DO 15 C 1 = 1,8 300 

00 150 J = 1, 4 301 

150 CAL C( I, J I = BLANK 302 

DO 160 1 = 1,8 303 

GOTO (151,152,153,154,155,156,157,1581,1 304 

151 IF (PWRSKP1 GO TO 16C 305 

DO 1151 J= 1, 4 306 

1151 CALC( I,J 1 = C1IJI 307 

GO TO 160 308 

152 IF (NRMSKP1 GO TO 160 309 

DO 1152 J= 1, 4 310 

1152 CALC( I , J ) = C2< J) 311 

GO TO 160 312 

153 IF (PRSSKP1 GO TO 160 313 

DO 1153 J = 1, 4 314 

1153 CALCI I, J ) = C 3( J ) 315 

GO TO 160 316 

154 IF ( PLT SKP (111 GO TO 1154 317 

CALC( I, 1 1 = C4( 11 3ie 

CAL C< 1,2) = C4( 2) 319 

1154 IF (PLTSKPI 21 1 GO TO 160 32C 

CAL C( 1,3 ) = C4( 3) 321 

CALC( 1,4) = C 4< 4) 322 

GO TO 16 C 323 

155 IF (PLT SKP (3)1 GO TO 1155 324 

CAL C( I, 1 ) = C 5( 1) 325 

CALC( 1,2) = C5( 2) 326 

1155 IF (PLTSKPI 4) I GO TO 160 327 

CALCI 1,3 ) = C 5( 3) 328 

CALC( 1,4 > = C 5( 4) 329 

GO TO 160 330 

156 IF (PLT SKP ( 5 ) ) GO TO 1156 331 
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CAL C( It 1) = C6< 1) 332 

CAL C ( 1,2) * C6< 2) 333 

1156 IF ( P L T SKP ( 6 ) > GO TO 160 334 

CALC( 1,3 ) = C6( 2) 335 

CAL C( 1,4) = C 6(4) 336 

GO TO 160 337 

157 IF ( P L T SKP ( 7 ) ) GO TO 11 57 338 

CAL C( 1,1) = C 7( l) 335 

CALC< 1,2 ) = C7( 2) 34C 

1157 IF (PL7 SKP( 8) ) GO TO 16C 341 

CALC( 1,3 ) = C7( 3) 342 

CALC( 1,4 ) = C 7 ( 4) 343 

GO TO 160 344 

158 IF ( P L7 SKP( 5 ) ) GO TO 16C 345 

C AL C ( I, 1 ) = C8( 1) 346 

CAL C( 1,2) = C8( 2) 347 

160 CONTINUE 348 

345 

WRITE CHECKED INPUT DATA 350 

351 

IF (NSI.EQ.ll WRITE (6,21) Rl»PO»MOLWT»(CALC(l»l ) ,1=1,4) ,R2,P3, 352 

1 CP, (CALC (2, I), 1 = 1,4 ),RDIF, PR, GAMMA, (CALC (3, I) ,1=1,4) .WIDTH, 353 

2 TO, MU, (CALC (4, I), 1=1, 4) ,RGAS, SPEED, (C AL C ( 5 , I ) , 1= 1 ,4) , LOSS , 354 

3 CAP V, ( C ALC ( 6, I ) , I = 1 , 4 ) , ( C ALC ( 7 , I ) , I = 1 ,4 ) , (CALC ( 8 , I ) , I * l, 4 ) 35 5 

IF (N S I .EO. 2) WRITE (6,23) R1 , PO , MOL WT , (C AL C ( 1 , 1 ) , I* 1 ,4) , P2, P 3, 356 

1 CP, ( CALC ( 2, I), 1 = 1,4) ,RDIF, PR, GAMMA, (CALC (3,1 > ,1 = 1,4) .WIDTH, 357 

2 TO, MU, (CALC (4, I >,1=1,4) ,RG AS, SPEED, (CALC < 5, I ) ,1 = 1 ,4) , LOSS , 358 

3 CAP V, < CALC (6,1) , 1=1,4) , (CALC ( 7, I ) , I = 1 ,4 ) , (CALC ( 8, 1) , 1 = 1, 4 ) 355 

WRITE (6,19) 360 

36 1 

BEGIN MAIN CALCULATION 362 

363 

200 DO 27 C J*1,NJ 364 

365 

CHECK FOR RESTARTING WHEN SOME CASES ARE FINISHED BY PREVIOUS 366 

RUNNING OF THE PROGRAM. IF JDONE IS GREATER THAN ZERO, READ X AND 367 

MACH NUMBER DISTRIBUTIONS CALCULATED BY PREVIOUS RUNNING OF THE 368 

PROGRAM. CALCULATE H DISTRIBUTION AND A*. 0 1 MENS I ON AL IZE X, X*, 365 

P*, ANC T*. 37C 

371 

IF (J.C-T.JOONE) GO TO 21C 372 

READ (5,2) JJ,N,PSTAR( JJ) ,XSTAR( JJI ,TSTAR 373 

IF (J.EQ.JJ ) GO TO 201 374 

WRITE (6,15) 375 

STOP 376 

C 377 

201 READ (5,3) ( X X ( I ) , I = 1 , N) 378 

READ ( 5,3) ( F M ( I) ,I = 1,N) 375 

DO 202 1=1, N 38C 

X X ( I ) = XX ( I )*RDIF 38 1 

HH( I ) = H1(J )+(RO-Rl«-XX(I )*SIGN)*SI N(ALPHA( J) ) 382 

202 CONTINUE 383 

PSTAR(J) = PSTAR(J)*PC 384 

XSTAR(J) = XSTAR( J >*RDIF 385 

TSTAR = TSTAR *T0 386 

ASTAR ( J ) = HH (N ) 387 
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IF (RO.Nfc.O.O) ASTARI J ) =2.0*P I *( R0+ XX ( N ) *SI GN ) *HH (N ) 388 

GO TO 2 AO 389 

390 

CALL SUBROUTINES TO DETERMINE SOLUTION OF MACH NUMBER EQUATION AND 391 

DISTRIBUTIONS OF PARAMETERS THAT VARY WITH X 392 

393 

2 10 CALL NCFLOW 394 

240 CALL 0 ISTS 395 

396 

CALCULATE PARAMETERS THAT VARY ONLY WITH FILM THICKNESS 397 

398 

245 M DOT! J I = WIDTH*H ( 6, J » *RHO (6, Jl *U(6 , Jl *CN VT (4 ,NSI 1 399 

0 ( J ) = CNVT( 5 ,NSI )*MDOT( J ) 400 

KN(J) = 2 .96* MACH (6»J ) /REP ( 6, J) 401 

LAMCAIJ) = KN( JI*HMEAN(J> 402 

FMU = ML 403 

IF (GAMMA. EQ. 1.41 F MU=CNVT ( 1 , NSI >*T( 6, J I **1 .5 / ( T (6 , J> +CNVT ( 2, NS 1 11 404 

RER(J) = RHO(6,J)*CAPV*HMEAN(J)/FMU/CNVT(6,NSI) 405 

IF(PWRSKP) GO TO 260 406 

407 

CALCULATE POWER > SHEAR HEAT, TORQUE, AND TEMPERATURE RISE DUE 406 

TO POWER DISSIPATION 409 

41C 

250 POW ER ( J ) = 0. 411 

IF (REP(6,J>.LE.RELAMI POWER! J » =F MiJ*FAREA*C APV**2 /H ME AN( J 1/ 412 

1 CN VT ( 7 , NSI 1 413 

HSHEAR ( J I * CNVT( 8, NSI )*POWER( Jl 414 

DELT(J) = HSHEAR ( J)/ABS(MDOT(J) 1/CP 415 

TORQU E ( J I * POWER ( J )*CNVT ( 9 ,NSI ) /SPEED 416 

417 

CALCULATE FORCE AND CENTER OF PRESSURE 416 

419 

26C K=0 42C 

KK=0 421 

FORCE ( J I = SIMPS1(0.,RDIF ,FFUNC , K) * WIDTH 422 

XC( J ) = SIMPS1(0.,RDIF, XCFUNC , KK)* WIDTH /FORCE ( J I 423 

265 IF IK.NE.O) WRITE (6,13) 424 

IF(KK.NE.O) WRITE(6»14) 425 

IF (NRMSKP) GO TO 270 426 

FBAR(J) = FORCE ( J ) /RDI F/ ( P0-P3) /WIDTH 427 

XCBAR(J) = XC(J)/RDIF 428 

270 CONTINUE 425 

275 CALL STFNSS 430 

431 

WRITE OUTPUT DATA AND PLOT THEM (SUBROUTINE GRAFIC) 432 

433 

IF (NSI.EQ.l) WRITE ( 6,25) 434 

IF (NSI.EQ.2) WRITE (6,26) 435 

DO 710 J=1,NJ 436 

WRITE (6,43) HMEAN( J ),HMIN( J) ,ALPHA(J) ,MDCT(J ) ,C( J) ,RER(J» , 437 

1 STIFF(J ), FORCE ( J ) , XC ( J) 43e 

IF (ABS(MACH( 11, J »- 1. CI.LT. l.E-51 WRITE (6,45) 435 

IF <REP(6,JJ.GT.RELAM.AN0.REP<6, J) .LT.RETURB) WRITE (6,46) 440 

IF (R EP( 6,J ).GE .RETURB ) WRITE (6,47) 441 

710 CONTINUE 442 

C 443 
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WRITE 16,271 444 

IF C .NCT.NRMSKP) WRITE (6,281 445 

IF (NSI.EO.l) WRITE (6,291 446 

IF (NSI.E0.2I WRITE (6,30 447 

IF ( .NOT .NRMSKP ) WRITE (6,311 448 

00 720 J = 1 , N J 445 

WRITE (6,481 HMEAN( J),HMIN( J) .ALPHA ( J) , KN ( J >, LAMDA ( Ji , 450 

1 XSTAR( J >,PSTAR( J) 451 

IF (.NCT.NRMSKP) WRITE (6,44) FBAR< J) , XCBAR I J ) 452 

IF <ABS(MACH< 11, J )- I. 0) .LT. l.E-5) WRITE (6,45) 453 

IF (REP( fc,J ).GT .RELAM. AND. REP(6, J) . L T. RETURB) WRITE (6,46) 454 

IF (REP(6,J).GE.RETURB) WRITE (6,47) 455 

720 CONTINUE 456 

C 457 

IF (PWRSKP) GO TO 740 458 

WRITE (6,32) 459 

IF (NSI.EO.l) WRITE (6,33) 460 

IF (NSI.E0.2) WRITE (6,34) 461 

00 730 J=1,NJ 462 

WRITE (6,43) HMEAN( J) ,HMIN( J) .ALPHA (J) ,PCWER( J) ,HSHEAR( J) ,DELT< J ), 463 

l TORQUE! J ) 464 

IF (ABSIMACHl 11, J )-1.0).LT. 1.6-5) WRITE (6,45) 465 

IF (P EP(6,J ).GT. RELAM. AND. REP(6,J) .LT. RETURB) WRITE (6,46) 466 

I F (R EP(6,J ) .GE .RETURB ) WRITE (6,47) 467 

730 CONTINUE 46e 

C 469 

740 IF (PRSSKP ) GO TO 800 470 

DO 750 J= 1, N J 471 

00 745 1 = 1, 11 472 

IF (NSI.E0.2) T ( I , J ) = T( I , J) —46 0. 473 

745 CONTINLE 474 

IF (MOC( J,3).E0.1) WRITE (6,35) 475 

Al = WICTH*HI( J » 476 

A 2 = W I0TH*H2< J ) 477 

IF (RO.NE.O.O) A 1=2 . 0*P I *R 1 *H 1 ( J ) 478 

IF (RO.NE.O.O) A 2=2 . 0*PI*R2*H2(J) 475 

0 ADR = ABS(A1-A2)/RDIF 480 

DELA = A8S(A1-A2)/AMIN1(A1,A2) 481 

IF (NSI.EO.l) WRITE (6,36) HME AN ( J) ,H 1 ( J) ,D ADR , AL PHA ( J ) , 482 

1 F2( J ) ,DELA 483 

IF (NSI.E0.2I WRITE (6,37) HMEAN( J ) , H 1 ( J ) ,D ADR , AL FHA ( J ) , 484 

1 F2( J I, DELA 485 

IF (ABS(MACH( 11, J )-l.O).GT. l.E-5) WRITE (6, 38 ) 486 

IF (ABS(P ( 11, J )-P 3) .GT.PTOL) WRITE (6,39) 487 

WRITE (6,401 488 

IF (NSI.EO.l) WRITE (6,41) 489 

IF (NSI.E0.2) WRITE (6,42) 490 

WRITE (6,43) (X(I),H(I,JI, MACH (I , J) , U ( I , J ) ,RH 0 ( I , J ) ,P(I,J), 49 1 

1 T< I, J),REP( I,J) ,FRICT(I,J),I=1,11) 492 

750 CONTINUE 493 

C 494 

800 CALL CRAFIC 495 

GO TO (100,110,120,130), INCODE 496 

49 7 

END OF PROGRAM 498 

495 
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1 FORMAT ( 12A6) 500 

2 FORMAT <2X, 13, 3X, I 3.5X.EI 4. 7, 5X , E14.7 ,5 X,fc 14. 7 ) 501 

3 FORMAT (5E14.7) 502 

10 FORMAT ( 1H1, 22HAREA EXPANSION PROGRAM, 5X, 12A6 ,/ ,LHO , 502 

1 35HFLOW LENGTH PARAMETERS INCONSI STANT, 5X ,4HR 1 =,G10.3,5X, 504 

2 4HR 2 =» G10.3, 5X, 7HR2— R1 = ,G1 0.3 ,5X ,7HWI DTH =,G10.3) 505 

12 FORMAT ( 1H1, 22HAREA EXPANSION PROGRAM ,5X , 12 A6 1H0 , 506 

1 26H INSUFF IC IENT PRESSURE DATA,5X,11HPG = P3 = 0 ) 507 

13 FORMAT < 1H0,53HINACCURATE NUMERICAL INTEGRATION IN FORCE CALCULATI 508 

ION I 509 

14 FORMAT ( 1H0.66HINACCURATE NUMERICAL INTEGRATION IN CENTER CF PRESS 510 

1URE CALCULATION! 511 

15 FORMAT ( 1F0,26HDISTRIBUTICNS OUT OF ORDER) 512 

19 FORMAT ( IH0,44X, 30H************************** **** , / ,1H ,44X,1H*, 513 

1 26X, IH* , /, IH ,44X,30H* / - CHOKED FLOW *, /,1H , 514 

2 44X.30H* ♦ - TRANSITION REGION *,/,lH ,44X, 515 

3 3 OH* T - TURBULENT FLOW *,/»lH ,44X , IH* , 28X , 1H*, / , 516 

4 IH , 44X , 3 OH **************************** ** ) 517 

20 FORMAT ( 1H1, 7 1HPR0GRAM FOR QUASI-ONE DIMENSIONAL FLOW WITH FRICTIO 518 

IN AND AREA EXPANSION, /, 1H0 ,2CX , 12 A6 , /> 519 

22 FORMAT ( IHO, 1 2H INPU T DATA - , / , 1HC, 1 OX ,9HR1 , ME TE RS ,22X , 7HP1 , NSMA , 520 

1 18X, 16HM0LECULAR WEIGHT, 1 7X , 9HF=K/RE**N , / ,IH ,E18.3, E30.3, 521 

2 F29. 3, 22X.9H*********,/, IHO, 10X.9HR2, METERS, 22X.7HP2, NSMA, 522 

3 17X,18HCP , JOULE S /KG— DEG K , 16X , 10HK ( LAMI NA P) , / , l H ,E18.3, 522 

4 E30.3,E 28.3,F30.3»/»1H0,6X,1 8HFL0W LENGTH .METERS , 18X , 524 

5 5HP 1/P2 » 24X » 5HGAMMA ,23X,10HN(LAMINAR),/ ,IH , E 18 .3 ,E29 .3, 525 

6 F29.3»F30.2»/»1H0»7X»1 7HF LOW WI DTH , METE RS , 16X ,8HT C, DEG K, 526 

7 17X,18HVISC0SITY,N-SEC/M2»9X» 24 HUP PER LIMIT RE (LAMINAR),/, 527 

8 IH ,EL8.3,F29.1,G32.4,F28.1,/,1HO,39X,10H LOSS C0EF..21X, 528 

9 SH SPEED ,RPS»19X,12HK( TURBULENT) ,/, IH ,F47 .2, E32 .4 ,G30.4, / , 529 

X 1H0.70X ,7HV, M/SEC, 21X,12HN(TURBULENT) ,/ ,1H ,F77. 2.G32.4, / , 530 

1 1H0.91X , 26HL0WE R LIMIT RE ( TURBULENT ),/, IH ,F107.1) 531 

21 FORMAT ( 1H0, 1 3H0UTPUT DATA -,/, IHO, 1 OX, 9HR1 .METERS , 22 X , 7HP C/NSMA, 532 

1 18X, 16HM0LECULAR WE IGHT ,17X,9HCALCULATE ,/ ,1H ,£18.3, E30. 3, 533 

2 F29 .3, 2 2X ,9H* ******** ,/,lH0»10X,9HR2,METERS»2 2X, 7HP3 /NSMA, 534 

3 17X, 18HCP .JOULE S/KG-DEG K , 9X, 4A6, / , IH , E 18.3 , E30 . 3, E28.3, /, 535 

4 IH ,92X ,4A6,/,1H ,6X,18HFL0W LENGTH , MET ERS , l 8 X, 5HP0/ P3, 24X, 536 

5 5 FGAMMA , / , IH ,E 1 8.3 , 2F29. 3 , 16 X,4A6 , / , IHO ,4X , 537 

6 22HMEAN FLOW WIDTH, METERS, 14X ,8HT0, DEG K,17X, 538 

7 18HVISC0SITY,N— SEC/ M2 , 19X.4HPL0T,/, IH , E 18.3 , F29 . 3.G32.4, 539 

8 14X, 22H******** ************** ,/,lHO ,30X , 540 

9 28HGAS CONSTANT, JOULE S/KG-DEG K , 12 X ,9HS PEED, RPS , 14X , 4A6, / , 541 

X IH , E48 .3 ,G32 .4 , / , IH ,93X ,4A6 , / ,IH ,39X,10HL0SS CGEF., 542 

1 21X.7HV, M/SEC, / , IH ,F 47. 2 , E31 . 3 ,15 X ,4 A5 , / , IHO ,93X ,4A6 ,/ , 543 

2 1H0.93X , 4A6, / ) 544 

23 FORMAT ( IHO, 1 3H0UTPUT DATA - , /, IHO, 10X ,9HR1 ,1 NCHES , 22X , 7HP C, PS I A, 545 

1 18X, 16H MOLECULAR WEIGHT, 17X.9HCALCULATE ,/ ,1H ,F 18 .4, F30. 3, 546 

2 F29.3.22X , 9 H* ******** ,/,lH0,10X,9HR2,INCHES,22X,7HP3,PSIA, 547 

3 18X, 16HCP,BTU/LBM-DEG R , 10X.4A6 , / , IH ,F 18 .4, F 30 . 3 ,F2 8 .3, / , 548 

4 IF ,92X ,4A6, /, IH ,6X,18HFLCW LENGTH, I NCHES.18X, 5HP0/ P3, 24X, 549 

5 5HGAMMA , / , IH ,F 18. 4 ,2F 29. 3 , 16 X ,4A6 , / , IHO ,4X, 550 

6 22HMEAN FLOW WI DTH, INCHES, 14X.8HT0, DEG R,16X, 551 

7 20HVISCOSITY, LB— SEC /IN2.18X, 4HPL0T , / , IH , F18 . 4, F29. 3 , G32.3, 552 

8 14X, 22H**********************, / , IHO ,32X , 553 

9 24HGAS CONSTANT, LB-FT/L8M— R,14X, 9HSPEE0 ,RPM, 14X , 4A6,/,1H , 554 

X F47.5,G32.4,/,1H ,93X ,4A6 , / , IH ,39X,10HLOSS CCEF., 555 
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1 2 IX » 8HV »FT/SEC»/»IH ,F A7. 2 ,F3 0. 2 , 16 X ,AA6 , / , 1H0 ,93X,A A6, / , 556 

2 1F0, 93X , AA6 , / ) 557 

2 A FORMAT < IHO, 12HINPUT DATA - , / , IHO, 10X ,9HR 1 , I NCHES ,22X , 7HP1 , PS I A, 556 

1 18X, 16HM0LECULAR WE IGHT,17X,9HF=K/RE**N ,/ ,IH ,F 18 .A, F30 . 3, 559 

2 F29.3,22X,9H*********,/,1HG, 10X ,9HR2 , INCHES ,2 2X , 7HP2 , PSI A, 560 

3 18X, 16HCP,BTU/LBM-0EG R , 1 7X, l OHM LAMI NA P) ,/ ,1H ,F18. A, F30.3, 561 

A F28.3,F3C.3»/»1H0»6X»18HF LOW LENGTH , I NC HES ,18 X, 5HPI/ P2, 2AX, 562 

5 5FGAMMA , 23X , 10HM LA MI NAR ) ,/,lH ,F 18 .A ,F 29 .3, F29 . 3 ,F30 .2/ , 563 

6 IHO, 7X , 17HF LOW WI DTH , INCHES , 16 X,8HTC ,DE G F,17X, 56A 

7 20FVISCOSITY, LB-SEC / 1 N2 , 8X »2AHUPPER LIMIT RE (LAMINAR),/, 565 

8 1H ,F18.A,F29.1,G32.A,F28.1,/,1H0,39X,10HL0SS C0EF.,21X, 566 

9 SHSPEED.RPM , 19X , 12H K( TURBULENT) ,/,lH ,FA7 .2 , E32 . A ,G3 0 . A, / , 567 

X 1F0,70X,8HV, FT/SEC, 20X , 12HN< TURBULENT ),/, IH , F77. 2, G32. A/ , 568 

l 1H0,91X , 26HLOWE R LIMIT Rfc ( TURBULENT) ,/ ,1H .F107.1) 569 

25 FORMAT ( 1H1, 5X , 7HH(CHAR ) , 7X ,6HH ( MI N) , 7X , 5HALPHA ,7X ,6HM ( DOT ) ,9X , 570 

1 IKO, 11X , 5HRE(R) , 6X, 9HSTIFFNESS,6X,5HF0RCE ,10X ,2HXC,/ , 1H ,5X, 571 

2 6HMETER S, 8X , 6HME TERS ,6X , 7HRADI ANS, 6X,6HKG/SEC,8X, AHSCMS , 2 3X , 57 2 

3 3HN/M , 1 IX, 1HN,10X,6H METERS) 573 

26 FORMAT ( 1H1, 5X , 7HH( ME AN) , 7X , 6HH (MIN) ,7X,5HALPHA,7X,6HM( DOT ) , 9X , 57A 

1 1HQ, 11X ,5HRE(R) , 6X, 9H STIFFNESS, 6X.5HF0RCE ,10X ,2HXC,/ , 1H , 6X , 575 

2 6HNCHE S, 7X,6HINCHES,7X,7HRADIANS,5X,6HLB/MIN,8X,AHSCFM, 2 2X, 576 

3 5FLB/IN, 10X , 2HLB , 9X , 6HI NCHE S) 577 

27 FORMAT ( IHO, 5X , 7HH( ME AN) , 7X ,6HH ( MI N) , 7X ,5HALPH A ,7X ,7HKNUDS EN, 6X , 578 

1 6HLAMBD A , 9X » AHX ( * ) » 9X , AHP ( *) ) 575 

28 FORMAT ( 1H*,97X , 5HF0RCE , ICX,2HXC) 58C 

29 FORMAT ( IH , 5 X, 6HMETE R S, 8X , 6HMETERS ,6 X , 7HRADI ANS ,6X ,6HNUMBER, 7X , 58 1 

l 6FMETER S, 8X,6HMETERS, 8X,AHN/M2) 582 

30 FORMAT ( IH , 5X , 6H INCHE S, 8 X , 6HINCHES ,6X , 7HRADI ANS, 6X ,6HNUMBER, 7X , 583 

l 6HNCHES,8X, 6HI NCHE S,8X,AHPSIA) 58A 

31 FORMAT ( 1H>,97X,5H(BAR),9X,5H(BAR) ) 585 

32 FORMAT ( 1H0,5X,7HH(MEAN), 7X,6HH(MIN) ,7X , 5HALPHA ,8X , 5HP0WER , 6X, 586 

l 8FH( SHEAR ),5X,8HDELTA(T),7X,6H TORQUE) 587 

33 FORMAT ( IH , 5 X, 6HME TE RS, 8X , 6HME TERS ,6X, 7HR ADI ANS, 7 X ,5HW ATTS , 7X , 588 

1 5 HWATTS , 8X, 5HDE G K,10X,3HN— M) 585 

3A FORMAT ( IH ,5X,6H INCHES, 7X, 6HINCHES,7X,7HRADI ANS, 9X ,2HHP,8X, 590 

1 7HBTU/M I N , 7X, 5HDEG F , 9X , 5HFT-LB » 591 

35 FORMAT ( INI) 592 

36 FORMAT < ihO,9HH(MEAN) =,G11.3,7H ME TERS, 5X , AH HI =,G11.3,7H METERS, 593 

1 5X , 1 1HDA /DR =,G11.3 ,5X,7H ALPHA =,Gll.3,8H RADIANS,/, IH , 59A 

2 32X.AHH2 =,Gll.3,7H METERS ,5X , 1 1HDA/A (M I N ) = , G1 1 . 3) 595 

37 FORMAT ( IH0»9HH(MEAN ) =,Gll.3,7H INCHES, 5X,AHH1 =,Gll.3,7H INCHES, 596 

1 5X, 1 1 HD A /OR =,G11.3,5X,7HALPHA =,G11.3,8H RADIANS,/, IH , 597 

2 32X , AHH2 =,Gll.3,7H I NCHE S ,5X ,1 1 HDA/A (M I N ) =,G11.3> 598 

38 FORMAT ( IE + , 9 2X , 1 5H( SUBSONI C FLOW)) 595 

39 FORMAT ( 1H*,9 3X , 13H(CH0KED FLOW)) 60C 

AO FORMAT ( IHO, 8 X , 1HX , 11 X, 1HH , 1 IX , AHMACH ,5 X , 8HVE LOCIT Y ,8 X , 7HDENS ITY, 601 

l 3X, 8HPR E SSURE , AX , 1 IHTEMPt RATURE »7X,5HRE ( P ) ,6 X ,8HFRICTI0N) 602 

A 1 FORMAT ( IH , 6X , 6HMETERS, 6X , 6HME TERS , 7X ,6HNUMB ER ,6X , 5HM/SEC , 10X , 603 

1 5FKG/M3, 6X, AHN/M2 , 9X , 5HDEG K,22 X,6HF ACT OR » 60A 

A2 FORMAT ( IH , 5X , 6HINCHE S, 7X , 6HI NCHE S , 7 X , 6HNUMB ER ,5X ,6HFT/SE C, 7X , 605 

l 11HLB-SEC 2/FT A, 3X,AHPSIA,9X,5HDEG F ,22X ,6HFACT0R ) 606 

A3 FORMAT ( IH ,9G13.3) 607 

AA FORMAT ( IH*, 9 IX , 2G 13. 3 ) 608 

A 5 FORMAT ( IH*, IH/ ) 605 

A6 FORMAT ( lh«, 2H ♦> 61C 

A 7 FORMAT ( 1H + .2H T) 611 

A8 FORMAT ( IH ,7G13.3) 612 

C 613 

END 6 l A 
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DETERMINE THE MACH NUMBER DISTRIBUTION SUCH THAT M=1.0 AT X=X*. 
FOR CHOKED FLOW, X*=RDIF - FOR SUBSONIC FLOW, P2=P3. 

PROGRAM VARIABLES 
* * * * * * * * * * * * * * * * * 


ENTRANCE CONDITIONS 
***** * ****** * ***** * 

HI - FILM THICKNESS 
A 1 - CROSS-SECTIONAL AREA 
FM ( 1 1 - INITIAL GUESS OF MACH NUMBER 
TT1 - TEMPERATURE 
PI - PRESSURE 

EXIT CONDITIONS 
*************** 


H2 - FILM THICKNESS 
A 2 - CROSS-SECTIONAL AREA 
P 2 — PRESSURE 
FM 2 — MACH NUMBER 

CONDITIONS AT CHOKING 

**** ** ****** ********* 

XSTAR - FLOW LENGTH 
HSTER - FILM THICKNESS 
ASTER - CROSS-SECTIONAL AREA 
PSTER - PRESSURE 
TSTAR - TEMPERATURE 

PROGRAM VARIABLES 
***************** 

DEL X - STEP SIZE 

N - NUMBER OF POINTS IN SOLUTION OF MACH NUMBER EQUATION 
ARITHMETIC FUNCTION FLGRNG DOES A 3 POINT LAGRANGE INTERPOLATION 
SUBROUTINE NCFLOW 

DOUBLE PRECISION X ,FM ,H , XO , XI , X2 , YO , Yl , Y2 , XX( 2 > ,PP ( 2) , F M2, P2, 

1 ASTER, PSTER , PI, A1,A2 

REAL LOSS, MU 

COMMON/ARRAYS/E 1( 231) ,FMM (1 1 ,20) ,E2( IA20) ,H U 20 ) ,H2 (20 ) , ALPHA! 20), 
1 FMEAN( 20), XSTAR (20) ,P STAR (2 0) , A STAR (20) 

COM MON /TR AYS/ N, X( 201) ,FM( 201) ,H ( 201 ) , J , TSTAR 

COM MON /CONSTS /GAMMA ,RQ IF, RO ,SIGN,REF(6),T0, PO , P3, PT CL , RGAS ,LOSS, 

1 MU,PI,RUNIV(2),CNVT(I1,2) ,NSl 

DIMENSION SX( 201), SY( 201) 

FLGRNG (X,X0,X1,X2,Y0,YI,Y2) = YC* ( X-Xl) * I X-X2 ) / (XO-XL )/ (X0-X2 ) ♦ 

1 Yl*( X-XO l*( X-X2) /( Xl-XO ) /( X1-X2)* Y2*( X-XO) *(X-X 1 ) /( X2-X0) / 

2 (X2-XL) 

c 

C SET BOUNDARY CONDITIONS - DEFINE CONSTANTS NEEDED IN SOLUTION OF 

C MACH NUMBER EQUATION 

C 


1 

2 

2 

A 

c 

6 

7 

8 
9 

IC 

11 

12 
13 
1 A 

15 

16 

17 

18 
19 
2C 
21 
22 
23 
2A 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3 A 

35 

36 

37 

38 

39 
AC 
A 1 
A2 
A3 
AA 
A5 
A6 
A7 
Ae 
A9 
5C 

51 

52 

53 
5A 

55 

56 

57 

se 


45 


o o o ooo oooo 


POLO = C.O 59 

EX = .50 6C 

NN = ALOGIO(RDIF) 61 

X TOL = POIF* 10 .0**<NN-5) 62 

FM( 1 I = .2500 63 

IF (J.NE.l) FM< 1 >=FMM(1,J-1)*H1< J)/H1( J — 1 } 64 

CON = S IN ( ALP HA ( J ) I 65 

X ( 1 ) = 0. 66 

Hill * HU) 67 

IF (SIGN.LT.O.C) H ( l ) =H 2 ( J » 68 

IP (RO.EO.O.O) GO TO 50 65 

A 1 = 2.C*P I*R0*H( 1) 7C 

IF (S ICN .GT.0.0 ) A2=2.0*PI*(RQ*RDIF )*H2 < J) 71 

IF (SIGN.LT.O.C) A2=2.0*PI*( RO-R DI F ) *H1 ( J ) 72 

GO TO 55 73 

50 Al * Ml) 74 

A2 = H ( 1 ) +RDI F*CON 75 

95 XX( 1) = C.O 76 

X X < 2 I = 0.0 77 

PP( 1) =0.0 78 

PP< 2) = C.O 79 

8C 

GUESS A STARTING /ALUE FOR THE ENTRANCE MACH NUMBER - CALL 81 

SUBROUTINE RK 1 TO FIND THE SOLUTION OF THE MACH NUMBER EQUATION 82 

83 

100 CALL RKl(CON) 84 

FMOLD = FM( 1) 85 

X STAR ( J ) = X ( N ) 86 

IF (XSIAR(J).LT.ROIF) GO TO 145 87 

8 8 

CALCULATE P2 89 

90 

120 TT1 = TO/( 1 . 0+ . 50* (GAMMA-l.O)M(FMll ) /LOSS)* *2 ) 91 

TSTAR = TT1*< l.*.5*<GAMMA-l.)*FM{L)**2)/.5/(GAMMA*l.) 92 

HSTER * MINI 93 

ASTER = FSTER 94 

IF (RC.NE.O.O) ASTER=2. C*PI*( RO+- XSTARIJ)*SIGN)*HSTER 95 

PI = PC/( 1. + .5*<GAMMA-1.)*(FM( l )/LOSS)**2) **(GAMMA/ (GAMMA-1. I ) 96 

PSTER = P 1*A1*FM< 1 ) /ASTER*SQRT< < 1. 0+. 50* ( GAMMA- 1 .0 ) *FM( 1)**2>/ 97 

l ,5/( GAMMA* 1. ) ) 98 

95 

13C FM 2 = 1 .DO IOC 

11 = 1 10 1 
IF (ABS(RCIF-XSTAR(J) I.LE.XTOL) GO TO 143 102 

135 DO 140 1 = 1, N 103 

II = I 104 

IF (ROIF-X(I)) 141,142,140 105 

140 CONTINUE 106 

141 IF (II.Lfc.2) 11=2 107 

IF ( 1 1 .GE.N ) II =N-i 108 

FM2 = FLGRNG( RDIF,X(I 1-1) ,X(II I.XlIIfl) ,FM( II -1 ) , FM (I I ),FM( II + l)) 109 

GO TO 143 11C 

142 FM2 = FM< II ) 111 

143 P 2 = ASTER *PSTER/A2/FM2*SQRT( . 5 0* ( G AMMA+1 .0 ) / ( 1 .0+ . 50 *< GAMMA-1 . 0) 112 

1 *FM2**2)I 113 

114 

FOR CFCKED FLCW, CMIT ITERATION FOR P2=P3 115 

116 

IF ( ABS ( XSTAR ( J l-RO IF ) .L6 -XTQL. AND. P2 .Gt «P3 ) GO T C 160 117 

IF (ABS(P2-P3 I .LT.PTGL) GO TO 160 118 
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119 

FOR SUESONIC FLOW, COMPARE P2 AND P3. IF THE Y ARE NOT EQUAL, SAVE 120 

THE VALUES OF P2 AND Ml. WHEN 2 SUCH POINTS HAVE BEEN SAVED, DC A 121 

LINEAR INTERPOLATION TO DETERMINE THE Ml THAT CCRRESPCNCS TO P2=P3 122 

123 

IF (PSTER.LE.P3I GO TO 146 124 

IF (P STER .LT . . 5DO*P2l EX=EX/2.0 125 

IF (ABSIPOLD-P2 l.LT.10.0*PT0L) EX=EX/2.0 126 

IF (EX.LT..1) E X=EX* 2. 0 127 

144 POLC = P 2 128 

145 FM ( 1 1 * FM(l)*(XSTAR< J)/ROIF)**EX 129 

GO TO 153 130 

146 IF (P2.GT.P3I GO TO 150 131 

XX( 1) = FM( 1) 132 

PP( 11 = P2 133 

IF ( XX ( 2 J.NE.O.d GO TO 152 134 

IF (SIGN.GT.O.O.OR.ABS(XSTAR( J l-RDI F > . LT. XTOL ) GOTO 151 135 

FM( II x FM( ll*(RDIF/X STAR (Jl)**.25 136 

GO TO 153 137 

150 IF ( P 2 .CT.PP (2)1 EX=EX/2. C 138 

XX( 21 = FM( 1 1 139 

PP( 2) = P2 HC 

IF (XX(l).NE.O.Ol GO TO 152 141 

151 FM( 1 ) = FM(1)*( P2/P3)**EX 142 

GO TO 153 143 

152 FM( 1) x (P3-PPI 21 1*(XX( 11-XXI2) 1 /< PP( l)-PP( 2) )+XX(2 ) 144 

153 IF <FM( ll.LT. 1.D01 GO TO 100 145 

FM(1) = ( 2.D04FMOLD+I.D01 /3.D0 146 

GO TO 100 14 7 

148 

SOLUTICN COMPLETE (SATISFACTORY X* FOUND AND P2=P3 FOR SUBSONIC 149 

FLOW 1 - PUNCH X AND M ON CARDS FOR RESTARTING 150 

151 

160 ASTAR( J 1 = ASTER 152 

PSTAR(J) = PSTER 153 

IND = C 154 

DO 170 Ixl, N 155 

IF (I.EQ.ll GO TO 165 156 

IF (I .EQ.N1 GO TO 166 157 

IF (X( I )-SX(IN01.LT.l.E-8> GO TO 170 158 

IF (FM( ll-SYt IND1.LT.1.E-8) GO TO 170 159 

165 IND * lND+1 160 

166 SX< IND l x x(I I 161 

9 Y ( IND J = FM( I » 162 

170 CONTINUE 163 

DO 175 1=1, IND 164 

SX< I) = SXm/RDIF 165 

175 CONTINUE 166 

PPUNCH = P STAR ( J I /PC 167 

XPUNCH * XSTAR( JI/ROIF 168 

TPUNCF = TSTAR/TO 169 

WRITE (6,101 J, IND, PPUNCH, XPUNCH, TPUNCH 170 

K = 0 171 

NNN = ( IND/5)* 5 1 ? 2 

NN * NNN+1 I 73 

00 180 1=1, NNN, 5 174 

II = 175 

WRITE (6, 111 (SX( JI,J=I,II» ,K 176 

K = K ♦ 1 177 
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18G CON T IN It 170 

II = IND-NNN 179 

IF (II.EQ.l) WRITE (6,12) ( SX ( J) , J=NN , I ND) , K 180 

IF (II.EQ.2) WRITE (6,13) ( SX( J) , J=NN ,1 ND) , K 181 

IF (II.EQ.3) WRITE (6,19) < SX ( J) , J = NN , I ND ) , K 182 

IF (II. ED. 9) WRITE (6,15) ( SX ( J) , J=NN , I ND ) , K 183 

K = K ♦ 1 189 

DO 185 1=1, NNN, 5 185 

II = 1+9 186 

WRITE (6,11) ( S V( J ) ,J=I ,1 I ) ,K 187 

K = K ♦ 1 188 

185 CONTINUE 189 

II = INC-NNN 190 

IF (II.EO.l) WRITE (6,12) ( SY( J) , J=NN , I ND) , K 191 

IF (II.EQ.2) WRITE (6,13) ( SY ( J) , J=NN , I ND ) , K 192 

IF (II.EQ.3) WRITE (6,19) ( SY ( J) , J=NN , I ND) , K 193 

IF (II.EQ.9) WRITE (6,15) ( SY ( J) , J=NN , I ND ) , K 199 

RETURN 195 

C 196 

10 FORMAT (2H* ,2HJ=,I3,3H N=,I3,5H P* =,E19.7,5H X* =,E19.7,5H T* =, 197 

1 E19.7) 198 

11 FORMAT (2H* ,5E19.7,I10) 199 

12 FORMAT (2H* , £19.7,166) 200 

13 FORMAT (2M- , 2E 19 . 7, I 52) 201 

19 FORMAT (2H* ,3E19.7,I 38) 202 

15 FORMAT (2E* ,9E19.7,I29> 203 

C 209 

END 205 


1 

CIS Tk I E L T IONS OF PRESSURE, IE MPEPA TUR E , DENSITY, MACH NUMBER, 2 

VELOCITY, REYNOLDS NUMBER, AND MEAN FRICTICN FACTOR 2 

9 

PROGRAM VARIABLES 5 

***************** ^ 

7 

N - NLMBER OF POINTS IN SCLLTICN CF MACH NUMBER ECUATICN t 

XX - ARRAY OF INDEPENDENT VARIABLE (DISTANCE) IN SOLUTION CF 5 

MACH NUMBER EQUATION 1C 

F v - maCH NUMBER ARRAY FROM SCLLTICN CF MACH NUMeER EQUATION 11 

12 

X - ARRAV CF EQUALLY SPACED DISTANCES 12 

MAO - MACH NUMBER AT X 19 

T - TEMPERATURE 15 

U - RADIAL VELOC I TV 16 

RFC - CENSITY IT 

F - PRESSURE 16 

REP - REYNOLDS NUMBER 15 

FRICT - MEAN FRICTION FACTOR 2C 

2 1 

ARITHMETIC FUNCTION FLGRNG DQcS A 3 FCINT LAGRANGE INT E P PC L AT 1 C N 22 

C 22 

S L P F n L TINE CISTS 29 

R EAL MACH, LOSS, ML, MLX 25 


48 


DOUBLE PRECISION FR FUNC , X X , F M , HH , XC , XI , X2 , VC , Y 1 , Y2 
COMMON /ARRAYS / X ( II ) ,P I 1 1 , 2C 1 , MACH U 1 , 20 ) , It 11 ,20) ,T (1 l ,20 J , 

1 R 1-0(11, 20, RE Pill, 2 0 ,FRICT (11 ,20) ,H(11 ,20) ,XC8AR <20 ), 

2 PC W£R ( 2 C ) , FORCE ( 2C >, S 1 1 FF < 20) ,F BAR ( 20 ), E X TR A ( 80 ), XST A P ( 20 ) , 

3 PSTAR(2C),ASTAR(2Ct 

COMMON /TRAYS/N,XX(2C1),FM(2C1),HH(2C1),J,TSTA» 

COMMON / CONSTS/CAMMA,ROIF , RC»SIGN»EX(8) , P3 , FTO L , EGAS ,L CSS, MU, PI, 

1 RLNIV( 2),CNVT (11,2) ,NSI 

FLGRNC-(X,X0,X1,X2,Y0,Y1,Y2) = Y C * ( X- X 1 ) * ( X- X2 ) / t X C- X 1 ) / (X0-X2>+ 

1 V 1* ( X— XO ) * ( X-X2)/{X1-XC)/(*1-X2) + Y2*(X-XC)*<X-X1)/(X2-XC)/ 

2 ( X 2— X 1 ) 

AM = l.C 
KK = 1 

ICC CO 16 C 1 = 1,11 
DO 12 C K - KK , N 
KK = K 

IF (X ( ))-XX(K ) ) 121 »122,12C 

1 2 C C0N7IMLE 
12 1 IF (KK .LI .2) KK = 2 
IF (KK.CF.N) KK =N— 1 

MACH 1,2, = ,'LGRNG(X( I ) ,XX(KK-1) ,XX(KK) ,XX(KK + 1),FM (KK-1), 

1 FM(KK),FM(KKM) ) 

GO TO lie 

122 MACH I , J ) = FM(KK) 

13C T ( I , J ) = .5C*(GAMMA+l.C)*TSTAR/(l.C+.50*(GAMMA— 1.0)* M ACH(I,J)**2) 
U(I,J) = MACH ( I , J )*SORT(CNVT(ll,NSl )*GAMMA*PGAS«T( I , J) ) 

IF (RO.NE.O.O) AN = 2.C*PI* (RC+X(I )*SIGN) 

P(I,J) = PSTARI J)*4$TAR(J)/AN/H (I , J ) * SORT < .50 * < C-AMM A+ 1 .0 1/ 

1 < 1 .C*. 50* (GAMMA- 1.0 )*MACH (I , J ) + *2 ) ) / W ACH ( 1 , J ) 

IF (KK.CT.2 50) GO TO 135 

IF ( I . L T . 1 1 ) GO TC 135 

IF ( A E S ( P ( 1,0 )-P 3 ) .LT.PTOU GC TC 135 

KK = 2 c C 

MACH I , J ) = l.C 

GO Tr 1 3 C 

135 R HO ( I , J ) = P( I,J)*CNV7( 1C,NSI ) /RGAS/T (I ,J> 

MUX = ML 

IF (GAMMA. EQ. 1.4) MLX =CNV T ( 1 , N S I ) * T (I , J ) + *1 .5 / ( T ( I , 2 ) +C NVT ( 2 , NS I ) ) 
R EP ( I , J ) = RHO( I, J )*U( I , J >*2. C*H (I , J) /MLX/CNVTI 6,N$I) 

FR I CT ( ] , J ) = FRFLNC <REP( I , J) ) 

16C CON TIME 

RETLPN 

END 


2 A 
27 
2 E 

2 S 

3 C 
3 1 
32 
3 3 

2 4 
; c 

36 

3 7 
3 E 

3 5 

4 C 

41 

42 

43 

44 

45 

46 
4 7 
46 

4 S 

5 C 
5 1 


54 

c c 

J „ 

56 
5 7 
56 

5 S 

6 C 
6 1 
62 
6 3 

64 
6 5 
66 

6 7 
6 6 

65 

7 C 
71 
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C 1 

PRESSURE CALCULATION FOR NL V ERI CAL INTEGRATIONS 2 

PROGRAM VARIABLES A 

****** i ** ****** ** * 

t 

X - GISTANCE ACROSS FACE CF SEAL 7 

XX - STOREC X OISTRIBLTICN FROM 0 TC X* 6 

F M - MACH NUMBER OISTRIBLTICN WITH XX S 

FF - FILM THICK NcSS 0 I S TP I B L T I CN WITH XX 1C 

1 1 

FMM - MACH NLMBER AT X 12 

F - F 1LM TFICKNE SS AT X 17 

AN - AREA AT X 14 

PRESS - PRESSLRE AT X 15 

It 

ARITFMET1C FUNCTION FLGRNG GOES A 3 PCINT LAGRANGE INT E RPC L AT I C N 17 

IE 

FUNCTION PRES S ( X ) IS 

DOIIEL E PRECISION XX ,F M ,HH , XC , XI , X2 , VC , Y1 , Y2 2C 

R EAL LC SS 2 1 

COMMON /ARRAYS/EXTRA (1971) ,PSTAR (20) ,ASTAR(20) 22 

CCMM0N/TRAYS/N,XX(2C1),FM(2CI),HH(2C1),J,TSTAR 23 

COMMON /C0NSTS/CAMMA,R0IF,R0,SIGN,DI6),TC,PC,P3,FT,RCAS,lCSS,MU,Fl, 2 < 

1 PLN I V ( 2 ),CNVT(1 1,2) ,NSI 25 

FLGPNG(X,XO,Xl f X2,YC,Vl f Y2) = V C* ( X- X 1 ) *1 X- X2 ) f ( XC- X 1 ) / (XC- X 2 1 ♦ 2t 

1 Yl*( x-xn )*( X-X2 ) /(Xl-XO) / (X1-X2) FY2* ( X-XC )* (X-Xl )/<X2-XC)/ 27 

2 (X2-X1) 2 E 

C 2 c . 

AN a 1 .C ‘ C 

DO 12 C 1 = 1 « N 3 1 

KK = I 32 

IF (X-XX( I ) ) 122. 121.12C 23 

12C CONTINUE 34 

GO TO 122 35 

12 1 FMM = FM (KK ) 3f 

H = FMKK ) 37 

GO TO 135 36 

122 IF (KK.LE.2) KK = 2 3S 

IF (KK .CE .N ) KK =N— 1 4C 

C 41 

1 3 C FMM = FL GPNG( X ,XX (KK- 1 ) ,XX (KK > , XX( KK+1 > ,F R( KK-1 ) ,FM (KK ) ,FM ( KK* 1 ) ) 42 

H = FL CPNC-( X, XX (KK- 1) ,XX( KK ) ,XX (KKF l) ,HH( KK-1 ) ,FH( KK) ,HF(KK + 1 ) ) 43 

C 44 

135 IF (RC.NE.C.O) A N =2 . C*P I * ( R C+ X * S I GN ) 45 

PRESS - PSTAP ( J )*ASTAR ( J) /AN/H* SCRT (. 50* (GAMMA* 1.0 ) / 4t 

1 ( 1 .C* .5C* (GAMMA- i.o *FMM**2) ) /FMM 47 

C 4 E 

RETURN 4 S 

END 5 C 
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OOURLE PRECISION RUNGE— KU1TA SCLLTICN CF THE MACH NL M B E P ECUATICN 
FOR FlCt* KITH AREA CHANGE AND FRICTION 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FOR FLCK KITH AREA CHANGE AND FRICTION 

BOUNCAPY CONDITION IS - Y=i.C AT X = X4. DECRcASf STEP SIZE UNTIL A 
SATISFACTORY X IS FOUND. 


WORKING VARIABLES ARE IN OCLRLt PRECISION 


PRO GR AN VARIABLES 

4 4 4 4444*444*4 


X - C I STANCE FROM INLET 
E - F ILM THICKNE SS 

Y - SOLUTION OF MACH NLMBEP 6QLATI GN 
R - RADIOS THAT CORRESPONDS TO X 
T - LOCAL TEMPERATURE 
T 1 - ENTRANCE TEMPERATURE 
Ml - VISCOSITY 
RE - REYNOLDS UNMBER 
F - FRICTION FACTOR 

CON - S INI TILT ANGLE) 

G - RATIO CF SPECIFIC HEATS 

GAMMA - RATIO CF SPECIFIC HEATS 

R 1 - RADIUS AT ENTRANCE 


YO 

R ECNST 
FACTOR 

D EL X 
HI 
DEL Y 
HR 


SQUARE OF ENTRANCE MACH NUMBER - INITIAL ESTIMATE 

SUPPLIED BY CALLING PROGRAM 

RHO*U*C 

MULTIPLICATIVE FACTOR FCR CHANGING YC TC ACCELERATE THE 
SOLUTION 

LENGTH OF X INTERVAL (STEP SIZE) 

FILM THICKNESS AT MIDPOINT CF X INTERVAL 
LENGTH OF Y INTERVAL 

RATIO OF FILM THICKNESS TO RADIUS FOR SEALS WITH RADIAL 
EXPANSION 


444444 

FK I * 
FK2 * 
FK2 * 
FKA * 
44444* 


INTERMEDIATE VALUES OF THE DERIVATIVE CY/CX IN THE 
RUNGE-KUTTA FORMULA 


ARITHMETIC FUNCTION DYDX DEFINES THt DERIVATIVE OF THE SCUAFE CF 
THE MACF NLMBER KITH RESPECT TC X 

ARITHMETIC FUNCTION HFUNC DEFINES THE FILM THICKNESS AS A FUNCTION 
OF X 


SUBROL T IN E RK 1 ( SCON ) 

COUPL E PR EC IS ION DYDX ,CON,F ,G , R , YG ,DELX , R 1 , FK I , FK2 , f K3 , FKA , H I , 

1 C ELY, FRF INC, HFUNC , S IGN ,HR , X TGL , X , Y , H 

REAL LCSS, MU, ML IN 

COMMON /TR A YS/N,X( 2C1),Y(2C1),H(2C1),J,TSTAR 

COMMON /CON STS /GAMMA ,RDIF,SR1,SSIGN,EX(6),TC,PC,F3,PT,RGAS,LCSS, 

1 ml IN, PI ,RLNI V ( 2 ) ,CMVT(II ,2 ) , NSI 

DYCX( Y,F,CON, F,G,H» ) = 2.0C4Y*(1 .DO.5D0* (G-l • CC ) * Y )4 (F*G4Y-CCN- 
1 )P )/H/( l.DO-Y) 

HFUNC ( X , H , CON ) = H + X4CCN 


1 

2 

•3 

A 

c 

t 

i 

E 

c 

1C 
1 1 
12 
13 

1 A 
15 
It 
17 
IE 
IS 

2 C 
2 1 
2 2 
2 2 

2 A 
25 
2t 
21 
28 
2S 

3 C 

31 

32 
2 2 
3 A 
35 
3 < 
37 
3 E 
3 S 
AC 
A 1 
A 2 
A 3 
AA 
A 5 
At 
A 7 
AE 
AS 
5 C 

51 

52 

C 1 

5 A 

C C 

5 1 
5 7 
5 E 
5S 

ec 
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61 

INITIAL STEPS - TRANSFER °ERMANENTLY STCREC SINGLE FRECISICN 62 

CONSTANTS TO DOLBLE PRECISICN VARIABLES 62 

64 

G = gamma £5 

R 1= SRI C 6 

CON = SCCN tl 

SIGN = SSIC-N fce 

FACTOR = IC.O f <; 

X( 1 ) * C .CC 7 c 

YO = YU 1**2 7 i 

HR = C.CC 72 

IF (R 1 .EQ.C.OC I SR 1=1. C 72 

74 

IOC Y< 1) = YO 7 * 

N ST = 11 76 

N 1 = 2 77 

N 2 = 11 76 

CEL X = R C IF /FACTOR 7 $ 

NN = ALCGIC(RCIF) eC 

X TOL = RC IF*10.C0**(NN-4> 61 

NOTE = C 62 

62 

RECNST = 2.0*CNVT( 1C, NSI l*PC*SRl*H (1»*SCRT< Y( II «GAMNA*CNVT ( 11, NS I I 64 

1 /RG AS/TO I /CNVTI 6 ,NS I ) /( 1. C+. 5C* (GAMMA-1 .0 )*Y ( 1) /l CSS **2 I ** 86 

2 ((GAMMA+l.C)/2.C/(GAMMA-i.CII 86 

T 1 = TC/( l.C*.EC*(GAMMA-l.CI*Y(l) /LCSS**2) 67 

86 

RUNC6-KIT7A SOLLTICN 6 S 

SC 

20C DO 20C I=N 1, N 2 c;i 

II = I 9 2 

XIII = X ( I- 1 l+DELX S3 

R = R HX( I- 1 1* S IGN s 4 

IF (Rl.NE.O.OC) HR=H(I-1J/R S 5 

T = T 1*( 1 .C+. 5C*( GAMMA-1. C)*Y < 1 } ) / ( 1. 0+.50* (G AMMA-1 .0 ) *Y( 1-1 1 I 96 

MIJ = MUN S7 

IF (GAMMA. EO. 1.41 M L=CN V T (1 , N S I )*T**1.5/ ( T* C N VT (2 ,NSI )) 96 

RE = PECNST/ml s 9 

IF (Pl.Nt.C.DC) RE=RE/R ICC 

F = FR FLNC(RE I 101 

FK1= EELX*DYCX(Y( 1-1 J,H( I — 1 > f CCN,F,G,HR) 102 

R = (X( 1-1I + .50C*0£LXI*SIGN+P1 102 

HI = F F INC ( X ( I— ll*.EDC*D 6 LX,H(l| ,C CM 1C4 

IF (Rl.NE.C.DC) HR =H I /R ICE 

T = T 1* ( 1 .CM ECMGAMMA-l. C) *Y ( l) ) /( 1.0U50* (G AMNA-i .0 I* 1C6 

l ( Y ( I— 1 ) ♦ • f 0 C*F« 1 ) ) 107 

MU = ml in XC 6 

IF (GAMMA .EQ . 1 .4) M U=C N V T ( 1 ,N 5 1) * T* * l . 5 / ( T* C N VT ( 2 ,NS I M 109 

RE = RECNST /ML IK 

IF (PI .NE.C.PC ) RE = RE/R HI 

F = FR F LNC ( P E ) U2 

FK 2 = CELX*PYnx(Y( I-i | + .5DC*FK1 ,HI ,CCN,F,G,HR ) 112 

T = T 1* ( 1 .C + . EC*< GAMMA-1. Cl *Y ( 1 }> /( 1. 0+.50* (GAMMA-1 .0 I ♦ 114 

l (Y( I- 1 )+ .EO0*FK2) I 11 * 

MU = Ml IN ik 

IF ( GAMMA .EO . 1 . 41 M L =CN VT ( 1 , N S I I * T* * 1 . 5 / { T+ C N VT ( 2 ,NS I 1 1 111 

RE = RKNST/Ml IK 

IF (Rl.NE.C.DC) RE = RE/R US 


52 


r> o r> r> no - ) r> o o o r> 


F = FR f INC (RE ) 12C 

FK3 = CELX*DYDX(Y< 1-1 )4 . 50 C*F K2 ,H I ,CCN,F,G,HR ) 121 

P = X( 11* SIGN 4R 1 12 2 

HI I ) = EFLNC( >( 1),H(1),CCN) 122 

IF (Rl.NE.C.DC) HR=H(I)/R 124 

T = Tl*{ l.C4.5C*(GAMMA-l.CHY(l))/(l.o .5CMGAPEA-1.0M 125 

l (Y( I-D + FK2)) 126 

MU = PUN 12? 

IF (GAMMA.EQ. 1.4) ML=CNVT (1,NSI >*T**1.5/ (T«-CNVT(2 ,NSI )) 126 

RE = RECNST/ML 12S 

IF (Rl.NE.C.DC) RE=RE/R 13C 

F = FRFLNCIRE ) 121 

FK4 = CELX*DYDX(Y(I-1)4FK3,H(I) ,CCN,F,G,HR) 132 

DEL Y= ( FK 1-*2.DC*( FK2+FK3) 4FK4) /6.D0 133 

Y( I ) = Y ( I-D+DELY 134 

2BC IF ( CA6S( l.DO-Y( I ) ) .LE. l.D-E) GC 10 600 135 

136 

FLOW ML SI REMAIN SUBSONIC LNTIL EXIT - IF Y(I) EXCEEDS 1.0 AT ANY 13? 

POINT, LOWER DE LX 136 

135 

IF < ( Y< I- 1 )*.5CC*FK1.GT. 1.0 C.CR. Y ( I - 1 H-.5C0+F K2 .GT . 1. CO ) .OP . 14C 

1 (Y( I— 1 )-*FK3.Gl.l.DC.0R.Y(I l.GT.l.DO) 1 GC TC 320 141 

IF (DELY.LE.O.DC) GO TO 41C 142 

IF (0ELY.GT..1DC) GO TO 32C 142 

IF (Y( I J.LT.O.COC) GO TO 32C 144 

IF (CAES( X( D-ROIF J.GT.XTOL) GC TC 290 14* 

II = I 1- 1 146 

IF (NOTE.EO.O) GO TO 315 14? 

29C IF ( I . ( E . 20 1 ) GO TO 4C0 146 

20C CONTINLE 145 

15C 

2 1C N 1 * N 2 ■* 1 151 

GO TO 2 2 C 152 

215 NOTE = 1 152 

22C IF ( II.CE.ll) N 1 = 11 154 

DEL X = CELX/iC.CO 155 

33 C N 2 = N 149 „ 156 

GO TO 2CC 15? 

156 

IF INITIAL GUESS CF YO WAS TCC LCW, RAISE IT 155 

16C 

4CC IF ( l.CC-Y( II ) .LE .. 5DC) GC TC 42C 16 1 

4 1C YO = YC* 1 .05 162 

GO TO ICC 162 

164 

IF NUMBER OF STEPS EXCEEDS THE SIZE CF THE ARRAY, ELIPIMTE POINTS 165 

WHERE TEE VALLE OF THE FUNCTION IS NOT CHANGING RAPIDLY 166 

16 7 

42C CCELX = RCIF/50.C 166 

INC = N ST 16 S 

N = 2 C 1 17C 

IF (N2.LT.2C1) N=N2 171 

DO 440 1=NST, N 172 

IF (CAES(X( II — RDIFJ.LE.XTCL) GC TC 430 172 

IF ( I.EG.l.OP .1 .EQ.N) GC TO 43C 174 

IF ( X ( 1 ) . C-T .X < IND- 1 l+DDE L X ) GC TO 430 175 

IF (Y( 1 ) .LE .Y( IND-1 )4. 20-1) GC TO 440 176 

4 3 C X( IND ) = X ( I ) 17 7 

H( IND 1 = E ( I ) 17 6 
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Y( INC ) = Y< I ) 17? 

inc = in c •* i iec 

??C CONI I N L E 161 

N ST = INC- 1 18 5 

N 1 = INC 182 

N2 = N H? 18? 

IF (NST.LE.90 ) GO TO 2 C C 18? 

FACTOR = FACTOR / 2 . 0 188 

GO TO ICC iei 

C 188 

C STAISFACTOPY X* FOUND - SCLLTICN COMPLETE 18? 

C 19C 

6CC N = I I 1? 1 

CO 8 5 C 1 = 1, N 192 

8 5 C Y( I) = CSQRTI Y( II ) 192 

IF (R 1 .EC .C.DC ) SR 1 = C.G 19? 

C 19? 

RETURN 198 

END 197 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NUMERICAL INTEGRATION BY SIMPSONS RULE 

INTEGRAL CF Y DX FROM XMI N TO XKAX = (H/31MY0+4Y1+Y2) 

WFERE YC = Y EVALUATED AT XMI N 

Y 1 = Y EVALUATED ATA < X MI N+ X * AX ) / 2 
Y2 = Y EYALLATEO AT XMAX 
H = (XMAX— XMIM/2 (STEP SI2E) 

CALL Y t C TOR VARIABLES 
♦ ♦jm***************** 

J - INDEX OF FILM THICKNESS 
XMA X - LOWiER LIMIT CF INTEGRATION 
XMAX - UPPER L I M I T OF INTEGRATION 
FUNC2 - FUNCTION SUBPROGRAM TC EVALUATE Y 

KEP - NUMERICAL CONSTANT TC INOICATE IF INTEGRATION IS ACCURATE 

PRCCRAV VARIABLES 
**>>** + < *****♦**♦<= 

V - INDEPENDENT VARIABLE IN INTEGRATION 
E - STEP SIZE 
A - YC 
8 - Y 1 
C - Y2 

P - 2* ( VA L LL CF INTEGRAND) 

** * 

E* - DIFFERENCE BETWEEN ANSWERS USING DIFFERENT STEF SIZES 
NE* 

*** 

/NS - ACCUMULATED ANSWER FOR MANY SLBINTERVALS 

N - SUBINTET VAL COUNTER (N Cc 20G M^ANS INTEGRAL IS 
INACCURATE ) 

T - ERROR TOLERANCE 


1 

2 

■a 

? 

c 

8 

7 

8 
9 

1C 
1 1 
12 
12 
1? 
1? 
18 

17 

18 
1? 
2C 
2 1 
2 2 
22 
2 ? 
2 ? 
28 
27 
26 
2 ? 
3C 
2 1 

22 
a a 

2 ? 
a c 

36 
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FPAC 


****** 
TEST* 
N T E ST* 
****** 
0 


FRACTION OF ERROR TOLERANCE APPLICABLE TC MTF) 
SU8CIVI SION 

TEST VALLE C CP ERRCR IN INTEGRAL 


TEST VALUE OF FINAL AN S PE R 


FUNCT IT N SIMPS 1 I XMI N, XNAX ,F LNC 1 » KER ) 

DIMENS ION V( 2CC J,H( 2CO) , A (2 CO) ,8 (2 CCI ,C(200) , F I 200 I ,E (200) , NEI 200) 
DIMENS ION V( 20C ),H< 2C0) ,A (2 CO) ,B (2 CC) ,C(200> , P { 200 ) ,E (2001 ,NE( ZOO I 
EQUIVALENCE ( E ,NE ) , (TE ST , NTEST) 


DEFINE STARTING VALUES 


T=3.0E-5 
V( 1 J=Xf> IN 

H( 1)=C ,5*(XMAX-XMIN1 
A ( 1 >=FLNC 1IXMIN > 

B( 1 )= FLNC II XM IN*H( l )) 

C( 1 1= F LNC 1( XMAX I 

p( i i=m 1 1*( Ai i ha.c*b i n+cm i 
E( 1 I=P(1 I 
ANS=P ( 1 1 
N = 1 

FRAC=2 . C* T 
1 FRAC=C .5*FRAC 


BEGIN INTEGRATION USING MORE SOeiNTERVALS PHERE VALLE CF INTEGRANC 
IS CHANCING RAPIDLY 

2 TEST=AES( FRAC*AN$ ) 

K =N 

DO 7 I = 1 , K 

IF (NTEST. GT.IAESINEI I)I» GC TC 7 
5 N = N •* I 

V ( N ) = V I I UFIII 
H(N )=C .5*M I I 
A(N) = B( 1 J 

B(N) = fLNCl(V(N )*H(N )) 

C(N )=C ( 1 1 

PIN )= M N )*(A(N 1 -» A . C*B (NI+CIN) ) 

Q=P I I ) 

HI I)=MN ) 

B I I )= FLNC 1( VI I HHI III 
C( I I — AIN I 

P( I )*M 1 J*( AI I )4A.C*B( I H-C I III 
0=P( I HPINI-Q 
ANS=AN S *Q 
E I I )=Q 
EIN )=Q 

IF (N.CE.2C0) GO TO 13 
7 CONTINUE 

IF (N.CT .K I CC TO 2 

0 = o.c 

DO 11 1=1, N 

11 Q=Q *E ( 1 ) 

12 IF (ABSIQ l-T*AeS( ANSI I IA,1A,1 

13 KER=K ER ♦ 1 


3 7 
36 
2 c 
AC 
A 1 
A 2 
A3 
AA 
A 5 
A 6 
A 7 
AE 
AS 
5 C 

51 

52 
5 3 
5 A 
5 5 

56 

57 

5 E 
5S 

6 C 
61 
6 2 
63 

6 A 

65 

66 
67 
6 6 
6S 

7 C 
71 
7 2 
73 
7 A 

75 

76 

77 

7 £ 
7S 

ec 

ei 

£2 

63 

8 A 

85 

86 
67 
EE 
8S 
SC 
SI 
S 2 

53 
SA 
95 
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sf 

ACCU M LL ATE FINAL ANSWER 97 

9 £ 

14 ANS=G.C 9 2 

CO 16 1=1 ,N ICC 

If ANS = ANHP( I) 1C1 

S IMPS l=(ANS+Q/2C.C)/3. C 1C 2 

1C 2 

17 RETURN 104 

ENC ICE 


PROGRAM VARIABLES 1 

#**<+**** fc****** 1* 2 


XX - CHARACTERISTIC FILM THICKNESS 
YY - SEALING DAW FCRCE 
EDY - AXIAL FILM STIFFNESS 
VAX - NUMBER OF FILM THICKNESSES 
MACE - MACH NUMBER 
PEP - PRESSURE FLOW RFYNCLDS NUMBER 

* 

X * 

Y * 

CY * 

A *** WORKING VARIABLES IN NUMERICAL D I FFE RE NT I AT I CN 
SI * 

SI * 

P 2 * 

KV * 

*** ** * 

LACRANCE NUMERICAL Cl FFEREN TI AT ICN CVER MAXIMUM CF * FCINTS 
TO CETERMINE AXIAL FILM STIFFNESS 

SUBROUTINE STFNSS 
REAL MACE 

DIMENSION X ( 2 C ) » Y<20), 0 Y ( 2C ) , A ( 5 , 5 ) 

COMMPN/CONSTS/C II 8) , RELAM , RETURB ,D2 110) ,CNVT<231 

COMMON /ARP AYS AC 2 ( 231) , M AC H ( 1 1 , 2 C ) , U A ( 66 C) , R 5P ( 1 1 ,20 ) , C5 ( 48 0 ) , 

1 YY(2C),DDY(2C),Dft8C>,XX(2C)»D7<6C) 

COMMON /TRAYS/C8 ( 1207) , MAX , T STAR 

ELIMINATE INVALID POINTS AND ARRANGE VALID POINTS IN ASCENC1NG 
ORDER. IF THERE ARE LESS THAN 2 VALID POINTS, NC C I F F E R cNT I AT I CN 
IS POSS 1BLE. 


DO SC M = 1 , max 

9 C DCYIM) = c. 

DO ACC 1 N C = 1 , 4 

MM i C 

ICC DO lie M = 1 , ma X 

GO TO ( 1 , 2, 2, 4 ) , I NO 
1 IF (MACH l 1,M ) .CE . l.C .ANO.PEP <6 , v) 


L=. RtLAM) 


4 

£ 

f 

7 

E 

c 

1C 
1 1 
12 
12 

14 
12 
If 
17 
16 

15 
2 C 
2 1 
2 2 
22 
24 
22 
2 f 
27 
2 E 
2 S 

2 C 

3 1 
32 
■> 

34 
2 2 

2 f 
27 

3 E 
2 S 

4 C 
4 I 
4 7 
42 

4 4 


GC TC 102 


o r> o o ooo ooo ooo 



GG Ttl nc 




A 5 

2 

IF (MACKll,M).GE.l.C.AND.REP(fc,M).GE.RETLRe) 

GC 

TC 

1C5 

At 


GO 70 11 0 




A 1 


IF (MACK ll,M).LT.1.0.AND.REP(fc,M».LE.RtLAM) 

CC 

TC 

105 

At 


GO TO 11C 




A c . 

A 

IF (MACK 11»M ).LT.1.C.AND.REP(6,M) .Gc.RETUPB) 

GC 

TC 

1C5 

5 C 


GO 70 11C 




5 1 

1C 5 

MM = MM 4 1 
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X ( M M ) = X X ( M ) 




5 2 


Y (MM ) = YY (M 1 




5 * 

11C 

CONTINUE 




c c 


IF (MM.L7.21 GO TC ACC 




5 6 

13C 

CALL SCR T X Y ( X » Y »MM ) 




5 7 
e F 


SET UP MATRIX OF X DIFFERENCES F CR EACH PCINT 

X (K) 



5 c 






t C 

20C 

N = M 1NC(MM,5) 




6 1 


DO 25 C K = 1 » MM 




6 2 


1ST = M A X 0 ( K— 2 » 1) 




t 2 


1ST = M 1 N C ( MM— N + 1 f I ST 1 




t A 


IN = ISl-tN-1 




£5 


DO 211 1 I = I ST , IN 




££ 


I = I I- 1ST+1 




£ 7 


CO 2 1 C J J= 1ST, IN 




£ £ 


J = J J- I ST * 1 




£5 


A( I,J 1 = X( I I )-X( JJ) 




7 C 

2 1C 

CONTINUE 




7 1 

21 1 

CONTINUE 




72 

7 - 


FORM SIMS AND PRODUCTS FOR DERIVATIVE FORMULA 




74 
1 c 

22C 

SI = c. 




It 


S2 = C . 




7 7 


P 2 = 1. 




7 E 


DO 23 1 11= 1ST, IN 




7 S 


IF (II.EO.K) GO TC 221 




ec 


I = I I- 1ST *1 




e i 


PI = X( 11 )— X( K ) 




82 


S 2 = 5 2-1. /PL 




62 


P 2 = P2XP l 




64 


DO 23 C JJ=1,N 




65 


IF ( I.NE.JJ1 P1=P1*A(I»JJ) 




66 

23C 

CONTINLE 




8 7 


SI = S]4Y( III/P1 




66 

22 1 

CONTINLE 




65 


IF (( N/2 »*2.NE.N) S2=-S2 




SC 






S 1 


CER IV A T 1 V E 




S 2 
S 2 


KY = S 2=*Y(K )+P 2*S 1 




S 4 


0Y( K 1 = KY 




<; c 

25C 

CONTINLE 




S 6 






S 7 


PUT CALCULATED DERIVATIVES IN CPCER TC CCRRESFCMC 

TC 

INFUT XX 

SE 


ARRAY 




SS 






ICC 

20C 

DO 22 C M = 1 , MA X 




1C 1 


DO 2 1 C 1 I = 1 , M M 




1 C 2 


IF (XX (M ) .NE. X( III ) GO TO 3 1C 




1C 2 
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IF <( N/2 )*2.NE.N) GO TO 211 1CA 

CCY(M ) = — 0 Y( II) 1C5 

GO TO ;2C ICC 

3 1 1 C CY ( M ) = CY ( I 1 ) 1 C 7 

GO TO 220 lCF 

3 I C CON T INI E ICS 

32C CONTINLE 1 1 C 

40C CQNTINIE 111 

C 112 

RETURN 112 

ENC 11A 


C 1 

C FUNCTIfN FOR FINDING LOCAL FRICTION FACTOR AS A FUNCTION CF 2 

C REYNOLCS NUMBER 2 

C A 

COUEt E PRECISION FUNCTION FRFUNC I RE ) - 

COMMON/CONSTS/EXK 4), XL AM , XT CRB ,C CNLA M , CCNT RB »RELAR ,RETUR8»cX2( 33) t 

DOUBLE PRECISION X1,X2,X3,Y2 ^ 

c e 

C LAMINAR FLOW S 

C IC 

IF (RE «GT .RELAR I GO TC ICC 11 

FRFUNC = CONLAM /RE**XLAM 12 

RETURN 13 

C 1A 

C TRANSITION FLOW 15 

C U 

IOC IF (RE .CE .RETURB) GO TO 11C IT 

XI = A LC G( REL AM ) 16 

X 2 = ALCG( RETUR B) IS 

X 3 = A LC G ( RE ) 2C 

Y 2 = ALCG(CONTRB)-XTURe*X2-2.CDC*(ALCG(CCNTRB/CCNLAN)«-XLAM*Xl- 21 

1 XT UR B*X2 )*< ( X3*( X3*( X3-1. 5 CO* ( X 1+ X2 ) 1 +3. ODD* X 1 *X2 ) 50 CO*X 2**2 22 

2 *< 3 .CCC*X l->2 ) ) 1/ ( X2- Xl)**3 22 

FRFUNC = CEXP ( Y2 ) 2A 

RETURN 25 

C 2 C 

C TURBULENT FLOW 2T 

C 2 E 

11C FRFUNC = C0NTR8/RE**XTURB 2S 

RETURN 3C 

C 21 

END 32 
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C 1 

C INTEGRAND OF FORCE INTtGRAL 2 

C - 

FUNCTION FFLNCm 4 

COMMON /CON STS /HOI 12),P3,ODD(3C» 5 

FFUNC = PPESSt > )-P2 6 

RETURN 1 

en n t 


c 1 

C INTEORANC of CENTER OF PRESSURE INTEGRAL 2 

C - 

FLNCTICN XCFUNC(X) 4 

COMNCN/CONSTS/CCt 12 ) , P 3 ,DOO t 3C ) - 

XCFLNC = (PRFSSt X J-P3 )*X * 

RETURN 1 

ENO £ 


c ] 

c ARRANGE POINTS (X,Y) IN CROER CF ASCENDING X 2 

C 2 

SUBROUTINE SO R T X Y ( X , Y , N ) 4 

DIMENS ION X( ICC », Y(10C) £ 

NN= N- 1 £ 

C 1 

DO IIC 1=1, NN £ 

I I = I £ 

CO IOC J = 2 , N 1C 

I Ft X ( J 1 .CE. X( I) ) GO TC IOC 1 1 

T= X( J ) 12 

X ( J )= X ( I ) 12 

X< 11= T 14 

T = Y( J ) 1£ 

Y ( j )= XII 1£ 

y t i >= t n 

IOC CONTINLE 1£ 

lie CONTINLE 1£ 

C 2C 

PETLRN 21 

ENC 22 
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i 

DUMMY PLOTTING ROLTINE - EACH USER VLST V.RITF HIS CWN SLHRCUT IKE 2 

TO FIT Tt-E PLOTTING EQUIPMENT AVAILABLE TC HI K 3 

4 

SU8ROL T IN E CP AF IC 5 

LOGICAL PWRSKP,PR$SKP,NRM$KP,PUSKP 6 

REAL MACE*LOSS,MU 7 

COMMON / ARRAYS /X( 11 ) ,P ( 11,20 ,MACH(l 1 ,20) ,L(Il ,20 ) ,T <1 1 , 20) , £ 

1 PEO( II, 20 >,RE°( 11,20 ,FRI CT<11 ,20 ,H(U ,20) ,XCHAP (20 1, S 

2 PCWER(2C), FORCE (20 .STIFF (20) »FBAR(20),E1 (607 ,HR E AM 20), 1C 

2 E 2 ( £ C ) 11 

CGMMPN /CCN STS /GAMMA ,RD I F , R 1 f P£CC M 6 1 , TC , F C , P3 ,PT , RG AS , L CSS , HJ , F I, 12 

1 PLNIV(2),CNVT(ll,2) ,NSf 12 

COMMON /CT ITLE/T ITLE ( 12 ) 14 

COMMON /LCGICL /PVR SKP, NPMSKP ,PRSSKP, PL 1SKP (0) IE 

16 

WRITE (6,10 17 

1C FORMAT < 1H1, 47EFL0T S CANNCT RE MADE BY THE SUPPCUTINt SUPPLIED) 1£ 

PETLPN IS 

END 2 C 


1 

BLOCK [ATA SUBRCL TINE FCR CONSTANTS 2 

2 

ELCICK [ATA 4 

COMM IN /CON STS /C( 17) ,P I ,RUNI V(2) ,CNVT(ll ,2 ) ,NS I 5 

COM MON /PPN T /C 1( 4) ,C2(4 ),C 2(4) ,C4(4) ,C5(4) »C6(4) ,C7 (4) ,C8(4), BLANK 6 

DATA P J , ( RUN IV ( 11,1=1,2) / 3.L415S27, 8.31436E3, 1545.4 / 7 

OATA (CNVTU, 11,1 = 1,111/ .145S1E-5, 110.3, 2*1., .5051554E2, £ 

16*1./ S 

DATA ( CN VT ( 1, 2 1,1=1,1 1 )/ 1.5763SE-1C, 1R8.6, 720.0 , 13.405823, 1C 

1 13.C83, 1728., 45.833333, 42.42, 33COO.O, 4.4756636, 32.174/ 11 

DATA ( C )( I 1, 1= 1,4 1/6H , 6 H PO V , 6 HE R , 6 E / 12 

DATA (C2( 1 1,1=1,41/ 6H0 I *F N S , 6 H I C N LE S , 6 H S CU AN ,6 FT IT 1 ES / 12 

DATA ( C 2 ( 1 1,1= 1,4 1 / 6 H PRESS, 6 HLRE I'll , 6 HSTPI BU ,61-TI OS / 14 

DATA (C4( 11,1=1, 41/6HPCVER , 6 H , 6 H XC(, 6 EBAR) / 1? 

OATA (C5( I 1,1=1, 41/ 6H FORCE ,£H , 6 H PRES, 6 *-SURE / 16 

DATA ( C6( l 1,1 = 1, 41/ 6H TEMPER 16 .HATLRE , 6 H DEN, 6 )-SITY / 11 

DATA (C7< I 1, 1 = 1, 41/ 6H MACH N. 6 HLMBER , 6 H RE NU, 6 H W BEP / 1£ 

OATA (C£( 11,1=1, 41/ 6HFRICT , 6 HFACTCP, 6 H , 6 R / 15 

DATA PLANK/ 6 H / 2C 

C 2 I 

END 22 
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APPENDIX D 


SAMPLE OUTPUT 


FHCuhAM hn" i»t, j S i -i ,.f- ■) I H l« .•> lu'.- 1 I U,n a 1 T n h -< 1 1_ T I L i'+ ANL Js:- '.ArTo ILi, 


I NPli I l.o T t 


KkuriLt H - A kc<* txrAhSIu'* Pi\ . UNITS 


.MFTf KS 

2*f-ci 


■ ) •‘itOt * bo 


■'III lCul A 1 - n: ImiiT 
2 o • ^c. u 


F = */•< -. **N 


K/'.frT Fl*S F2. I 

t,*iU/p-bl JdoJC^bU 


LPi JiJULtS/kw.to h 
U« 1 U J L » U ' t 


K ( LA* INAr. I 

^.000 


1 tM.Th.MF 1 f M 5 

0. 


U. 


PI /P2 


••'id. 4.-; t,\,4 5 » 
1. JO 


Ml Ik IfkS 

0. 


LOSS ULI- 
1- Ju 


VlM.CS I TV ( N-o'.C/ v -< 
U 


SP^ Eu, RF» 

J . 


V, K/SEC 
Oj.Uu 


LIMIT :ir ( LA > 

Zi J J. J 


K ( TjrfeJLL.'STl 
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HlHEAM 

Hi HI N) 

ALPHA 

PCdER 

HI $HtAR> 

DELTAIT i 

TORQUE 


INCHES 

INCHES 

RAC I AN S 

HP 

BTu/hlN 

DEG F 

FT-LB 


C • 300F-O 3 

C-27dc-Cj 

-0. 1CCE-G2 

0.825E-C2 

0.350 

13.96 

0.391E-01 


C.4QCF-03 

0. 375E-G3 

— C • iGOt- J2 

J.olbE-02 

3.261 

4.851 

0. 292E-01 

/ 

C. *>006-03 

Q • 47 56 — 03 

-C. 1006-02 

O.409E-C2 

0.2ud 

2.284 

0.232E-01 

/ 

C.&0CE-03 

C. 5/5E-C3 

-C • 1006-02 

J. 4056-02 

0.172 

1.293 

D.192E-01 

/ 

C.7CCE-t)3 

0*6? 5E — G3 

-C. 100E-02 

0. 3446-02 

0.146 

0.622 

0. 163E-01 

/ 

C •tt0CF-03 

0. 7756-03 

-C. A00t-02 

3. 299 E- 02 

0.12/ 

0.367 

0. 142E-01 

/ 

C. 9006-03 

C. 8756-03 

-C. 1GCE-02 

0-264E-02 

0.112 

0.413 

0.125E-01 

/ 

C. 1006-02 

0.97 56—03 

-c.iooe-02 

J.237E-02 

0.100 

0.315 

0. 112E-01 

/ + 

C.l ICE-02 

C. 1C7E-G2 

-C. 1CGE-U2 

0 

0 

0 

0 

/♦ 

c.iioe-o? 

0. 1 1 76-02 

— C • 1GGE-U2 

0 

0 

0 

0 

/♦ 

C.13QE-J2 

0- 128E-02 

-C. 1C0E-02 

0 

0 

0 

0 

/T 

C .1406-02 

C. 137E-G2 

-C. 1CCE-02 

o 

0 

0 

0 

/T 

C.150E-02 

C. 1476-02 

-C.1GGE-02 

0 

0 

0 

0 

n 

C.16CE-Q2 

C. 158E-C2 

-C. lCGE-02 

o 

0 

0 

0 


H 1 HE AM * i). 

3006-03 INCHES 

HI * 

0.3256-33 INCHES 

DA/UR 


* 0. 18dE-Q 1 

ALPHA 

= -0. 100E-02 

RADIANS 



H2 = 

O.275E-03 INCHES 

LA/ A t M IN 1 

* 0.16 4 


(SUBSONIC FLOW) 


X 

H 

► ACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE ( P ) 


FRICTION 

INCHES 

INCHES 

N0H8EH 

FI /SEC 

L8-SEC2/FT4 


PS I A 

DEG F 



FACTOR 

C 

0. 325E-03 

C.l 04 

120.7 

C. 9bb£— 02 


64.51 

98.79 

159.9 


0. 150 

C.500F-02 

C.3206-G3 

0.110 

12 7.1 

0.933L-J2 


62. U 

98.66 

159.7 


0.150 

C .100E-01 

0.31SE-C3 

C. 116 

134.6 

0.8936-02 


59.47 

98.49 

159.5 


0. 150 

C.L50E-01 

0.3106-03 

0.124 

143.3 

C.83oE-02 


56. 5tt 

98.29 

159.3 


0. 151 

C.200F-01 

0. 305E-03 

0*133 

154.3 

0.802E-02 


53.38 

98.02 

159. L 


0.151 

C.25QF-01 

C. 300F-U3 

0.145 

167.8 

0.749E-02 


49.79 

97.66 

159.0 


0.151 

C.300F-01 

0.2956-03 

0.160 

185.4 

0.68dE-02 


45.70 

97.14 

158. 8 


0. 151 

C.35OE-01 

C.290E-03 

C.182 

210.0 

O.ol 76-02 


40.93 

96.33 

158.0 


0.151 

C.400E-J1 

0.265E-03 

C.214 

247.0 

0.532E-U2 


35.17 

94.90 

158.8 


0.151 

C.450E-01 

C*2bQE^QS 

£.276 

31 7 .9 

0.421E-02 


27.68 

91 .59 

159.3 


0. 151 

C.500E-D1 

0.21 56-03 

0.509 

573.5 

0.23OE-02 


15.00 

72.43 

163.5 


0. 147 

H I HE AM * 0. 

4006-03 INCHES 

HI » 

0. 425E-0 3 INCHES 

OA/DR 


= Q.182E-01 

ALPHA 

= - 0. 100E-02 

RADIANS 



H2 - 

0.275E-03 INCHES 

OA/AIHIN) 

1 0.116 


( SUBSONIC FLOW) 


X 

H 

HACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE ( P ) 


FRICTION 

INCHES 

INCHE S 

NtHBER 

HT/StC 

LB-SEC2/FT4 


PSIA 

DEG F 



FACTOR 

C 

0.4256-0.3 

C.173 

200.1 

0.9396-02 


63 . 66 

96.6 7 

344.5 


0.697E-01 

C.5C0E-J2 

0.4206-03 

0. 182 

21 J-v) 

0.924E-02 


61.26 

96.33 

344.1 


0« 697 E- 01 

C.IC0F-01 

0.415E-03 

G.l92 

221.4 

0 . 88 36-02 


38.66 

95.92 

343*8 


0.698E-C1 

C.lbOF-OI 

C. 4106-03 

C.203 

234.9 

0.843E-02 


53. d4 

95.41 

343.5 


0.699E-01 

C.20Cf-0I 

Q.405E-03 

0 • 2 1 7 

251.1 

0*7976-02 


52.74 

94.75 

343*3 


0.699E-C1 

.2506-01 

C.400F-03 

0.235 

2 71.1 

0.747E-02 


49.29 

93.88 

343.2 


0.699E-01 

< 3CCF-0I 

C. J95fc-C 3 

0.236 

2S/.0 

0*6d9t-02 


45.40 

92.66 

343.2 


0.699E-01 

{ .3506-01 

C. 3906-0 * 

0.269 

332.2 

0*62 36-02 


40.90 

90.81 

343.6 


0.699E-01 

1.40CF-OI 

0.3B56-C3 

0.336 

38 5.0 

0.544E-02 


33.50 

87.67 

344.6 


0.697E-01 

C.4SCE-01 

C . 3 80F —0 a 

C. 42 1 

4 79. d 

0*44 lt-0*. 


28.46 

80.84 

347.4 


0.691E-01 

C.300F-UI 

C. »7 5E-G3 

0.776 

8 50.9 

0.232E-J2 


15.00 

39.75 

368.8 


0.651E-C1 

HIHFAM = U. 

5006-03 INCHES 

HI = 

0.525E-J3 INCHES 

DA/ DR 


= 0.175E-01 

ALPHA 

* -0.1OOE-O2 

RADIANS 


H2 = U.475C -03 INCHtS 0 A/ A ( M IN ) = 0.886E-01 


X 

H 

HACh 

VFLCCITY 

DENSITY 

PRESSURE 

TEMPERATURE 

RE ( P) 

friction 

INCHtS 

INCHE S 

NLMHtR 

HT/SEC 

Ld-SEC2/FT 4 

PSIA 

DEG F 


FACTOR 

C 

0.525F-U 3 

L . *40 

277 .2 

0. 9‘*ofc-02 

62.44 

93.60 

583.9 

0.41 IE- 01 

C.50GF-02 

C. 6 206-0 J 

C* *51 

2b9. / 

0.9 ldE-02 

60.17 

93.02 

583.4 

0.41 IE- 01 

£ .1 UOF-dl 

C. 5 l 46-03 

C • *fc4 

3 04 .0 

0.8/ 76-02 

57. 7j 

92.31 

583.1 

0. 4126-01 

C.l 5CF-0 1 

C. 5106-03 

U.* 74 

J2J . 7 

0. dddE-u2 

55.09 

91.44 

582.9 

0.412E-01 

C .2GGF-01 

C.505E-C > 

0.29c 

j4o .6 

0.79oc-02 

52.21 

90.35 

582.9 

3.412E- 01 

C.25CF-01 

C.500E-C3 

C.318 

464 ■ d 

0.7496-02 

49.02 

88.92 

583.2 

0.412E-01 

C.iDOF— Oi 

0. 49 56-03 

C. 345 

49 J .5 

0.697E-U2 

45.45 

66.98 

583.9 

0.41 IE-01 

£ . 15LF-IU 

C • 490C-0 3 

0* 2*2 

4 4C .4 

0.63/E-02 

41.34 

04.16 

5E5.3 

0*41 36- Cl 

C .4006 - .J i 

0.464t-O3 

0.435 

494.4 

0.5o6E-02 

36*42 

79.57 

508.3 

0.408E-C1 

C 50F — 0 I 

t • 4ttuF-C3 

0.5 2c 

55o .4 

0.9 /4C-02 

29.99 

70.39 

5S5.2 

0.403E-C1 

C.500F- H 

C. 4 i 5E-U3 

1 « OC0 

u « l LoE* C4 

G.2/OE-02 

15. UJ 

6. 66 7 

657.1 

3.365E-01 


HlftfAM = l.feOOF 

-0) INCHES 

HI = 
H 2 — 

J.626t-J3 INCHES 
U. 5756-03 INCHES 

Ua/uR 
LiA/AIM IN) 

= 0.169E-01 

= 0.706E-01 

ALPHA 

= - 0 • 100E-02 RADIANS 
(CHOKED FLOW) 

X 

H 

PALM 

VbLLO ITY 

DENSITY PRESSURE TEMPERATURE 

RE < P) 

FRICT ION 

Ini. he s 

1 Nf.Ht S 

nlhcEh 

F r/SEc 

L B-SEC2/HT4 

PSIA 

DEG F 


F AC TOP 

i 

0 • ft 2 66 — 0 j 

L. jo 1 

444.9 

0.9 jlt-U 2 

6 1.04 

90 .04 

fc 5 7. 5 

0. 2807-Cl 

C .40cF-i1/ 

C.u2cE-U3 

C. 313 

46 J • 3 

J.9 J Jt-32 

58.95 

89 .22 

857.2 

0.280f-Gl 

f . 1 l.LF-iJ 1 

C.ft l 6 F-(. 1 

0 . 2/t 

i / a • D 

0 • oti dE -02 

40. 70 

60 .2 3 

857.1 

0. 280E- Cl 

0 . 1 5 0 E - 01 

L • c l OH — uj 

C * J*»4 

4 94 • 0 

0.d32E-Ui 

54.24 

87.05 

85 7.2 

0.280 " - Cl 

l ) 

L • n u i 

C • 36 4 

* i 0-4 

0 . / Y**i. - J C 

51*65 

85.57 

85 7.6 

3.28 3! - 31 

C./SOt-DI 

t' • oinJt- -ii 4 

C • jo 7 

4 • O 

0 • 752E-U2 

4d. 75 

03.70 

058.6 

0.280 r -01 

I . titOF-.) i 

G . r»9 5t-u 4 

C • 4 1 / 

*♦ 7 4 • 1 

o. /U jL -02 

4 5*51 

81.21 

d60« 3 

0. 279F-01 

l . m.r-n 1 

C.ISDc-u i 

0.45 5 

Jit .J 

0 * 65 2t -02 

4l.o J 

77.71 

86 3. 3 

3. 278' - Cl 

C.4U JH- ) I 

0 . »M •»* » 

0 • ‘J J u 

4/7.0 

G «469c-U2 

4 /. 37 

72.29 

868.8 

0.276' -Cl 

l 

i . *> o u r - v 4 1 

L • 1 0 * 

O 74 • 1 

0 • 5 Ode -0 2 

3 1. 61 

62.19 

bd0.4 

0. 27i r - CL 

( . j 1 . 1 . r - (J 1 

■* 7 »f-u 4 

1 .Ccu 

O . 1 cofe. 04 

C. 32ol-o2 

ld.il 

6 « 66 7 

9 o 3 . 4 

3.25 3 l — Cl 


62 



MMEAM = 0*1?8£ 

-04 METERS 

HI = 

0*1846-04 METERS 

OA/OR 


* 0.413E-03 

ALPHA 

* -0.100E-02 RAO IANS 



H2 = 

u. 171E-04 METERS 

OA/AI MIN) 

* 0* 579E-01 


ICHOKEO FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

REIPI 

FRICTION 

METERS 

METERS 

NUMBER 

M/SEC 

KG /Mi 


N/M2 

OEG K 


FACTOR 

0 

0*1846-04 

0*354 

123*7 

4*717 


0*4116+06 

303*5 

0.115E+04 

0*2086-01 

G.127E-03 

0*1836-04 

0*368 

128*3 

4.573 


0*3986+06 

302*9 

0*1156+04 

0.208E-01 

0*2546-03 

0*1826-04 

C.383 

133*4 

4*420 


0.383E+06 

302*3 

0*1156+04 

0* 208E- 01 

0.381E- 03 

0*1806-04 

0*400 

139*4 

4*255 


0.368E+06 

301*4 

0*1156+04 

0.208E-01 

C.5C8F-03 

0*1796-04 

0*421 

14o*3 

4.077 


0.352E+06 

300*5 

0.115E+04 

0*2086-01 

0*6356-03 

0*1786-04 

0.445 

154*5 

3*882 


0.333E+06 

299.2 

0.116E+04 

0.2086-01 

0*7626—03 

0.1776-04 

0*476 

164*5 

3*665 


3.313E+06 

297*6 

0.116E+04 

0.207E-01 

0*8896—03 

0.1756-04 

0*515 

177.4 

3*419 


0.290E+06 

295*5 

0*1166+04 

0*2066-01 

0*1026—02 

0* 1746-04 

0*569 

195*0 

3*129 


0.262E+06 

292*2 

0.1176+04 

0* 205E- 01 

C.l 146-02 

0*1736-04 

0*656 

222*7 

2*755 


0.227E+06 

286*4 

0*1196+04 

0* 202E-01 

0* 1276—02 

0*1716-04 

1.000 

322*8 

1*912 


0.142E+06 

259.3 

0*1286+04 

0* 187E-01 

H1M64M ■ 0*203E 

-04 METERS 

HI * 

0*2106-04 METERS 

UA/OR 


* 0. 3976-03 

ALPHA 

* -0.100E-02 RADIANS 



H2 * 

0* 197E-04 METERS 

OA/AIMIN) 

* 3.4B5E-01 


ICHOKEO FLOW) 


X 

H 

MACH 

VELOCITY 

OENSITY 

PRESSURE TEMPERATURE 

RE I P ) 

FRICTION 

METERS 

METERS 

NUMBER 

M/SEC 

KG/M3 


N/M2 

OEG K 


FACTOR 

0 

0*2106-04 

0*401 

139*4 

4*637 


0*4016+06 

301*4 

0.146E+04 

0* 164E- 01 

C. 1276-03 

0.208E-04 

0*415 

144*2 

4*506 


0.389E+06 

300*8 

0.1466+04 

0* 164E-01 

0*2566—03 

0*2076-04 

0*430 

149*5 

4*366 


0*3766+06 

330*0 

D* 146E+04 

3*1646-01 

0*3616-03 

0*2066-04 

0*448 

155*5 

4*217 


0.362E+06 

299* 1 

0.146E+04 

0* 164E-01 

0*5066-03 

0*2046-04 

0*469 

162*4 

4*055 


0.347E+06 

298*0 

0* 147E+04 

0* 164E-01 

0*6356-03 

0*2036-04 

0*494 

170*6 

3*878 


0.330E+06 

296*6 

0*1476+04 

0* 163E-01 

0*762 E— 01 

0*2026-04 

0*524 

180*5 

3*683 


0*3126+06 

294*9 

0.147E+04 

0. 1636-01 

0*8896-03 

0*2016-04 

0*563 

193*0 

3*461 


0.291E+06 

292*6 

0.148E+04 

0* 162E-01 

C. 1026-02 

0*1996-04 

0*615 

209*8 

3*200 


0*2666+06 

289*2 

0* 149E+04 

0* 161E- 01 

0*ll4E-32 

C* 1486-04 

0*698 

235*6 

2*864 


0*2336+06 

283*5 

0* 151E+04 

0* 159E- Cl 

0*1276-02 

0.1576-04 

1.O00 

322*8 

2* 100 


0* 156E+06 

259*3 

0.162E+04 

0.148E-01 

MMEAM * 0*2296-04 METERS 

HI * 

0* 235E-U4 METERS 

OA/OR 


* 0.3B1E-03 

ALPHA 

* -0.100E-02 RADIANS 



H2 * 

0*2226-04 METERS 

OA/AI MIN) 

* 0*4126-01 


ICHOKEO FLOW) 


X 

H 

MACH 

VELOCITY 

OENSITY 

PRESSURE TEMPERATURE 

RE( PI 

FRICTION 

METERS 

METERS 

NUMBER 

M/SEC 

KG/M3 


N/M2 

OEG K 


FACTOR 

0 

0*2356— 0* 

0*441 

153*1 

4*561 


0.392E+06 

299.5 

0 • 178E+04 

0*1356-01 

0*1276-03 

0.234E-04 

0*455 

157*8 

4*441 


O.381E+06 

298*7 

0 • 178E+04 

0.135E-01 

0. 2546-03 

C.232E-04 

0*471 

163*1 

4*311 


0.3696+06 

297*9 

0.178E+04 

0. 135E-31 

0*3816-03 

0*2 3 IE— 04 

0*490 

169*1 

4.177 


0.356E+06 

296.9 

0.178E+04 

0* L35E-01 

0*5066—03 

0*2306-04 

0*511 

176*0 

4*029 


0.3426+06 

295.7 

0*1786+04 

0. 135E-01 

0*6356-03 

0* 229E— 04 

0.535 

184*1 

3*869 


0.327E+06 

294*3 

0* 179E+04 

0.134E-01 

0*7626-03 

0.227E-04 

0*565 

193*7 

3*691 


0.310E+06 

292*4 

0 • 179E+04 

0. 134E-01 

0*8896-03 

0*2266-04 

0*603 

205.7 

3*490 


0*2916+06 

290*1 

0.180E+04 

0*1336-01 

0*1026-02 

0*2256— 04 

0*653 

221*6 

3*253 


0*2686+06 

286*7 

0*1826+04 

0.132E-01 

0*1146-02 

0*2246-04 

0*731 

245*5 

2*948 


0*2386+06 

281*1 

0 • 184E+04 

0* 130E-01 

C. 1276-02 

0*2226-04 

1.000 

322*8 

2*252 


0* 1686+06 

259.3 

0.196E+04 

0.1225-01 

MHEiM » 0.254E-04 METERS 

HI « 

0*2608-04 METERS 

OA/OR 


* 0.3656-03 

ALPHA 

* -0*1006-02 RADIANS 



H2 = 

0*2486-04 METERS 

OA/AI MIN) 

« 0.354E-01 


ICHOKEO FLOW) 


X 

H 

MACH 

VELOCITY 

OENSITY 

PRESSURE TEMPERATURE 

REI P) 

FRICTION 

METERS 

METERS 

NUMBER 

M/SEC 

KG/M3 


N/M2 

OEG K 


FACTOR 

0 

C.26UE— 04 

0*477 

164.9 

4.491 


0.384E+06 

297*6 

0.210E+04 

0* 114E-01 

0*1276-03 

0*2596-04 

0*491 

169*6 

4.380 


0.373E+06 

296*8 

0.210E+04 

0* I14E-01 

0*2546— 03 

0*2586-04 

C* 507 

174.9 

4*263 


0.362E+06 

295.9 

0.210E+04 

0. 1146-01 

0*3816-03 

0*2576-04 

0*525 

180*8 

4*137 


0.350E+06 

294*8 

0.210E+04 

0* 114E-01 

C.5C66-03 

0.2556-04 

0*546 

187*5 

4*002 


0.337E+06 

293*6 

0.2116+04 

0* 1146-01 

0*6356—03 

0*2546-04 

0.570 

195*4 

3*855 


0.323E+36 

292* l 

0.211E+04 

0.114E-31 

0*7626-03 

0*2536-04 

0*599 

204*7 

3*693 


0.308E+06 

290*3 

0.212E+04 

0.113E-01 

0*8896-03 

0*2516-04 

0*636 

216*1 

3*509 


0.290E+06 

287*9 

0.213E+04 

0.113E-01 

C. 1026-02 

C.250E-04 

0*684 

231 *2 

3*292 


0.269E+06 

284*5 

0* 2156+04 

0.112E-01 

0*1146-02 

0*2496-04 

C* 757 

253*4 

3.014 


0*2426+06 

279*1 

0.218E+04 

0.1106-01 

0.1276-02 

0* 248E— 04 

1.000 

322*8 

2.375 


0.177E+06 

259*3 

0.230E+04 

0* 1046-01 

HIME4M - 0.279E-0* METERS 

HI * 

0*2866-04 METERS 

OA/OR 


= 0.350E-03 

ALPHA 

* -0.1306-02 RADIANS 



H2 * 

0*2736-04 METERS 

DA/AIMIN) 

* 0*3076-01 


ICHOKEO FLOW) 


X 

H 

MACH 

VELOCITY 

OENSITY 

PRESSURE TEMPERATURE 

REIP) 

FRICTION 

HFTFfiS 

METERS 

NUMBER 

M/SEC 

KG/M3 


N/M2 

OEG K 


FACTOR 

0 

0*2866-04 

0*500 

172*4 

4*443 


0.378E+06 

296*3 

0.2396+04 

0. 104E-01 

0*1276-03 

0*2846-04 

0*514 

177.1 

4*340 


0*3686+06 

295*5 

0*2396+04 

0* I04E-01 

0*2546-03 

0*28 36-04 

0.525 

182*2 

4.230 


0*3586+06 

294*6 

0* 240E+04 

0* 104E-01 

0*3816-03 

0*2826-04 

C. 547 

187*9 

4*113 


0*3476+06 

293.5 

0*2406+04 

0*1056-01 

0*5086—03 

0.281E-04 

0*567 

194*4 

3*987 


0.335E+06 

292*3 

0.240E+04 

0* 1056- 01 

0*6356-04 

0*2796-04 

0*591 

202*0 

3*850 


0.321E+ 06 

290*8 

0.241E+04 

0.1056-01 

0*7626-03 

0*2786-04 

0.019 

210.9 

3*699 


0* 307E+06 

289.0 

0.242E+04 

0*1056-01 

0*8896-03 

C. 277E-04 

0.654 

221*8 

3*527 


G.290E+06 

286*6 

0.243E+04 

0.1056-01 

0*1026-02 

0.2766-04 

C. 700 

236*1 

3*324 


0.270E+06 

283*4 

0.245E+04 

0.105E-01 

C • 1 146-02 

0.2746 -04 

0* 769 

257 *1 

3*062 


0. 2456+06 

278*2 

0.248E+04 

0* 105E-C1 

0. 1276-02 

0.2736-04 

1.000 

322 *8 

2*446 


0.182E+06 

259.3 

0.262E+04 

0* 105E-01 
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AM = 0.3U5E 

— 04 METEkS 

Hi = 

0. 31 IE-04 MEIERS 

DA/ JR 


= 0.334E-33 

ALPHA 

= -O.lOOE-02 RADIANS 



H2 = 

0.298E-04 METEftS 

DA/AI MINI 

= 0.268E-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE(P> 

FRICTION 

METERS 

METEftS 

NUMBER 

M/S6C 

KG/M5 


N/M2 

DEG K 


FACTOR 

0 

0.3116-04 

C.511 

1/6.1 

4.419 


0.37 5E + 06 

295.7 

0.265E+04 

0. 106E-01 

0.1276-03 

0.3106-04 

0.525 

160.7 

4.318 


0. 3656+06 

294.9 

0.265E+04 

0.106E-01 

C .2 5 46- 03 

0. 3096-04 

0.540 

185.B 

4.211 


0.355E+ 06 

293.9 

0.265E+04 

0. L06E-01 

0. *81E-03 

C.307E— C4 

0.556 

191.4 

4.097 


0.344E+06 

292.9 

0.266E+04 

0. 1066-01 

C.5066-U3 

0.3066-04 

C.57B 

197.6 

3.974 


0.333E+06 

291.6 

0.2666*04 

0. 106E-01 

C. 6356-03 

0.3056-04 

0.601 

205.3 

3.841 


0.320E+J6 

290.1 

0.267E+04 

0. 106E-01 

C. 7626-03 

0.304E— 04 

0.629 

214.0 

3.693 


U • 306E+ 06 

288.3 

0.268E+04 

0. 1066- 01 

C.B89E-03 

0. 3026-04 

0.663 

224.8 

3.526 


0.289E+C6 

266.0 

0.269E+04 

0. 106E-01 

0. 10/6-02 

C. 3016—04 

0.7C8 

236. / 

3.329 


0.2706*06 

282.7 

0.271E+04 

D« 136E- 01 

0.1146-02 

0. 3G0E-04 

C. 776 

259.2 

3.074 


0.2456*06 

277.7 

J.275E+04 

0. 1066-01 

C. 12 76-02 

0.2S86-G4 

1.000 

322 .8 

2.476 


0.1846*06 

259.3 

0.2896*04 

0. 107E-01 

:4M * 0.3106 

-04 METEftS 

Hi - 

0.337E-J4 METERS 

DA/ DR 


* 0.3186-03 

ALPHA 

* -0.100E-02 RADIANS 



H2 * 

0. 3246-04 METERS 

DA/AIMIN) 

* 0.2356-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE I P ) 

FRICTION 

MFTFRS 

METEftS 

NUMBER 

M/S6C 

KG/M3 


N/H2 

DEG K 


FACTOR 

0 

0.3376-04 

0.522 

179.6 

4.396 


0. 372E+06 

295. 1 

0.291E+04 

0.107E-01 

0. 1276—03 

0*3356-04 

0.535 

164.1 

4.297 


0.363E+06 

294.2 

0.2916*04 

0.137E-01 

0.254F-03 

0.3346-04 

0.551 

189.1 

4.193 


0.3536*06 

293.3 

0.292E+04 

O.iOTE-Oi 

C.3B1F-Q3 

0.3336-04 

0.566 

194.8 

4.081 


0.342E+06 

292.2 

0.2926+04 

0. 107E-01 

C.5C8F-03 

C. 3316-04 

0.568 

201.1 

3.961 


0.331E+06 

291.0 

0.293E+04 

0. 107E-01 

0. *356-03 

0.3306-04 

0.611 

208.5 

3.831 


0.3186+06 

289.5 

0.2936+04 

0.1076-01 

0.7626-03 

0.3296-04 

0.639 

217.1 

3.667 


0.3046*06 

287.7 

0.294E+04 

0.107E-01 

0.6896-03 

0.3286-04 

0.672 

227.7 

3.524 


0.289E+06 

285.3 

0.2966+04 

0. 107E-01 

0. 1026-02 

0.3266-04 

0.717 

241.4 

3.331 


0.270E+36 

282.1 

0.298E+04 

0.107E-01 

0.1 146-02 

0.3256-04 

0.784 

261.5 

3.083 


0.245E+06 

277. 1 

0.302E+04 

0. 107E-01 

0.1276-02 

0. 324E-04 

1.000 

322.8 

2.503 


0.1866+06 

259.3 

0.3186+04 

0. 105E-01 

•AM « 0.3566 

-04 METEftS 

HI * 

0. 362 E— 04 METEftS 

DA/OR 


* 0.3026-03 

ALPHA 

-0.1006-02 RAOIANS 



H2 * 

0.349E-04 METEftS 

DA/AIMIN) 

* 0.2076-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

REIP) 

FRICTION 

ME I FAS 

METEftS 

NUMBER 

M/SEC 

KG/MJ 


N/N2 

OEG K 


FACTOR 

0 

0.3626-04 

0.5 34 

183.8 

4.368 


0 • 3696*06 

294.3 

0. 3196*04 

0. 1056-01 

0.1276-03 

0. 3616-04 

0.548 

188.3 

4.272 


0.3606*06 

293.5 

0.319E+04 

0. 105E-01 

0.2546- J3 

0. J59F-Q4 

0.564 

193.3 

4.170 


0.3506*06 

292.5 

0.320E+04 

0.105E-01 

0.3616-03 

0. 3386-04 

0.561 

198.8 

4.062 


Q.340E+06 

291.4 

0.3206+04 

0. 105E-01 

C. 5086-03 

0. 3576-04 

0.601 

205.1 

3.945 


0.3296* 06 

290.2 

J.321E+04 

0. 1056-01 

0.6356-03 

0.3566-04 

0.623 

212 .3 

3.819 


0.3166+06 

288.7 

0.3216*04 

0. 105E-01 

0.7626-03 

0. 3546-04 

0.650 

220.8 

3.680 


0.303E+06 

286.8 

0.323E+04 

0. 105E-01 

0.6696— 03 

0. 3536-04 

0.684 

231.2 

3.522 


0.2886*06 

284.5 

0.J24E+04 

0.105E-01 

0.10/6-02 

0.3526-04 

0.727 

244.5 

3.336 


0.2696+06 

281.3 

0.3276+04 

0.105E-01 

C.l 146-02 

0. 3516-04 

0. 792 

264.0 

3.097 


0.246E+06 

276.4 

0.33 1E+04 

0.1046-01 

C • 1 2 76-02 

0. 3496-04 

1.000 

322.6 

2.53d 


0.L89E+06 

259. 3 

0.347E+04 

0.103E-01 

AM > 0.3616 

-04 METERS 

HI * 

0. 3876-04 METEftS 

DA/DR 


» 0.2866-03 

ALPHA 

- -0.1006-02 RAOIANS 




H2 » 

0. 3756—04 METERS 

OA/ A ( M IN ) 

» 0. 1836- 

01 


(CHOKEO FLOW) 


X 

H 

MACH 

VELOCITY 

OENSITY 

PRESSURE 

TEMPERATURE 

REIP) 

FRICTION 

METERS 

METERS 

NUMBER 

M/aEC 

KG/M3 


N/M2 


DEG K 


FACTOR 

C 

0.36 76-04 

0.547 

187.8 

4.340 


0. 3666*08 


293.6 

0.348E+04 

0.103E-01 

0. 127F-03 

0. 3866-04 

0.561 

192.3 

4.247 


3.357E+06 


292.7 

0.3486*04 

0. 103E-01 

0.254F-03 

0.365E-C4 

C.576 

197 .2 

4.149 


0-347E+06 


291.8 

0.348E+04 

0. 103E-01 

C.381F-03 

0. 3846-04 

0.593 

202.7 

4.044 


0.337E+06 


290.7 

0.349E+04 

0.103E-01 

C.5C86-U3 

0.3826-04 

0.612 

208.8 

3.931 


0.327E+36 


2B9.4 

3.349E+04 

0.103E-01 

0.635F— OH 

0.3816-0* 

0.635 

215.9 

3.B09 


0.3156*06 


287.9 

0.350E+04 

0. 1036-01 

L.762F-03 

C. 8806-04 

0.661 

224.2 

3.6/4 


0.302E+U6 


286. 1 

0.351E+04 

0.103E-01 

C. 8696-03 

0. 3 786—04 

0.694 

234.4 

3.52c 


0.287E+06 


263.8 

3.353E+04 

0.1026-01 

0.1026-02 

0.3776-04 

0.737 

247.4 

3.543 


0.269E+06 


280.6 

0.3566+04 

0. 102E-01 

0.1 14F-02 

0. 3766-04 

C.800 

266 . 3 

3.1lc 


0.2466*06 


275.8 

0.360E+04 

0. 1026-01 

C. 1276-02 

0.3750-04 

1.000 

322 .8 

2.572 


0.1916*06 


259. 3 

0.377E+04 

0.101E-01 

•AM ■ J • 4066 

-04 METEftS 

HI = 

0.4136-34 METERS 

DA/OR 


« 0.270E- 

-03 

ALPHA 

- -0. 100E-02 RADIANS 



H2 » 

0. 4006-04 METERS 

DA/AI MINI 

* 0.1626- 

-01 


(CHOKEO FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE 

TEMPERATURE 

REIP) 

FRICTION 

M6TENS 

MfcTEKS 

NUMBER 

M/SEC 

KG/Mi 


N/M2 


DEG K 


F ACTOR 

0 

0.413E-04 

0.55b 

191.o 

4.31* 


0 • J63E+06 


292.8 

0.3766+04 

0. 1016-01 

0.I27F-O3 

0.41 16-04 

0.572 

196.0 

4.224 


0.3546+06 


292. 3 

3.376E+04 

0. 101E-01 

0.2546-03 

C. 4106-04 

0.5B7 

2 00.8 

4.128 


0.3456+06 


291.0 

0.377E+04 

0. 101E-01 

C. 3616-03 

U. 4096—04 

0.604 

20o .2 

4.027 


0. 335E+C6 


269. 9 

0.377E+04 

0. 1016-01 

0.5086-03 

C .*0 BE— 04 

0 . 62 J 

212 .3 

3.91o 


0.3256*06 


288. 7 

3* 378E+04 

D.131E-01 

C.ftJSt- 03 

C • *086-04 

0.645 

219 .2 

3.799 


0.3136*06 


287. 2 

0.379E+04 

0.1016-01 

f ./02F-03 

0.4086-04 

0.6 72 

227.4 

3.669 


0.301E+06 


265.4 

0.3806+04 

0. 1016-01 

( .8896-03 

0.40*E-U4 

0. 7C4 

237.3 

3.522 


J.286E + 06 


283* 1 

3 • 382E+04 

3.1336- 31 

0.1026-02 

0.4036-0* 

C. 745 

250.0 

3.348 


0.2696+Ou 


260.0 

0.385E+04 

0. 1006-01 

L.l 146-02 

C.401E-04 

0 • toG 7 

2 68 • 4 

3.12 + 


0.247E+06 


275.3 

0.389E+04 

0. 100E-01 

L. 1 2 76-02 

C. 4006-04 

1.000 

322.8 

2*602 


0.194E+J6 


259.3 

0.438E+04 

3.9896-02 
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program Fnn ouasi-gne d imens ic nal flc« *ith fhicticn and area exfansicn 


SAMPLE PRCBLEM - AREA EXPANSION PROGRAM - u. S. UNITS 


INPLl CATA - 


R?. INCHES 
3-31 5 C 


CP.BTU/LBM-DEG R 
3*240 


K(LAHINAR) 

24.000 


FLCM length, inches 

0 . 


GAMMA 

1.400 


N(LAMINAR) 

1.00 


FLOW WIDTH. INCHES 


VISCOSITY, LB-SEC/ IN2 
0 


UPPER LIMIT RE (LAMINAR) 
2300.0 


LCSS COEF. 
I *00 


SPEED, RPM 


K{ TURBULENT) 
0.7900E-01 


V. FT/SEC 

200.00 


N( TURBULENT) 
0.2530 


CUT PH CATA 


MOLECULAR WEIGHT 
28.966 


LOWER LIMIT RE (TURBULENT) 
3000.0 


RP. INCHES 
3.3 15C 


FI CM length, inches 
0.0500 


P0/P3 

4.333 


GAMMA 

1.400 


POWER 

DIMENSIONLESS QUANTITIES 
PRESSURE DISTRIBUTIONS 


MEAN FLOW t> ICTH. INCHES 
20.6717 


VISCOSITY, LB-SEC/IN2 
J.275E-08 


GAS CONSTANT. LE-FT/LEM-R 
53.35221 


SPEED, RPM 
6966.1 


LCSS COEF • 

1.00 


V, FT/SEC 

200.00 


*************** *************** 

• / - CHCKED FLOW ♦ 

• ♦ - TRANSITION REGION • 

• T - TURBULENT FLOW • 

****** *•*««*♦*** ************** 



M MF AN ) 

Hi MIN) 

ALPHA 

Ml DOT) 

u 

RE ( R ) 

ST IFFNESS 

FORCE 

XC 


INCHES 

INCHES 

RAO IANS 

LB /MI N 

SCFM 


LB/IN 

LB 

INCHE S 


C.300E-0i 

C.275E-Q3 

— 0 • 10QE-32 

0.1 C4 

1.36 7 

94.74 

0. 42 3E*04 

33.34 

0. 194E-01 


C.40QF-03 

0.375F-03 

-C.10UE-02 

0.224 

2.935 

126.6 

0. 164E+04 

33.05 

0.196E-C1 

/ 

C.500E-03 

0.475E-C3 

-C. 100c— 02 

0.379 

4.955 

159.9 

-555.0 

33.01 

0« 2306- Cl 

4 

C. 6 OOF- 03 

0.575E-C3 

-C. 1006-02 

0.553 

7.291 

194.0 

0.118E+04 

32.95 

0.205E-01 

4 

C.700F-0 4 

0.675E-C4 

-t.i00c-02 

0.740 

9.687 

228.3 

0.26 JE+04 

32.75 

0.2 39E-01 

4 

C.BCCF- 03 

C.775E-04 

— C. 1C0E-U2 

0.534 

12.2 2 

262-4 

0.348E + 04 

32.45 

0-2126-01 

t 

1.900F-03 

C.,8 3 5E— C 3 

-C.10LE-02 

1. 131 

14.79 

296.3 

0.395E+04 

32.07 

0.2156-01 

4 

C.100E-02 

O.S 7 56—03 

-C. 1CGE-Q2 

1 .3*9 

17. 38 

330.0 

0* 415E+04 

31.66 

0. 217E-C1 

4* 

C.110E-02 

0. 107E-02 

-C. 1006-02 

1.509 

15.74 

363.8 

0 

31.45 

0.219E-01 

7* 

C.I20F-0P 

C. 1 1 7E-G2 

-C.ltUt-02 

1.669 

21.83 

396.6 

0 

31.26 

0.223E-01 

/♦ 

C • 1 3 OF— 02 

C. 1266-02 

-t. U0E-02 

i.oil 

23.9o 

429.2 

0 

31.07 

0.220E- 01 

/T 

C • 140E—02 

C.137E-0? 

-C . 1COE-4J2 

2.0CJ 

26.20 

461.8 

0.217E+04 

30.84 

0* 22 IE- Cl 

/T 

C .1506-02 

C . 1 4 7F-CP 

-C. 1006-02 

2.176 

28«4 i 

494.6 

0.208E+04 

30.63 

0.222E-C1 

41 

( . 16CE-02 

C. 15EE-CP 

-C. 1C0E-02 

2.35 1 

JC * 7 6 

527.2 

C.158E+C4 

30.42 

0.222E-01 


Hi HE AN ) 

h( HI N) 

ALPHA 

KNLOSEN 

LAMBDA 

X 4 * > 

P( *) 

FORCE 

XC 


INCHES 

INCHES 

RADIANS 

NUMBER 

INCHES 

INCHES 

PSIA 

( BAR ) 

( BAP ) 


C. 40CF-03 

0.275E-03 

-C. iLGfc-02 

0.2 70E-G2 

O.eiuc-Oo 

0.5J9E-01 

7. 169 

3*645 

0. 389 


C.4O0F-O3 

O.J 756-03 

-C. 14)06-02 

0.2C3E-C2 

0.81 lc-06 

0 .503c- 0 1 

11.26 

0.6 39 

0.393 

t 

C.5CGE-U3 

0.N75E-C3 

-C . iCUE-02 

U.161E-02 

O.BOot-U6 

0 • 5006- 0 1 

15.00 

3.039 

0.401 

4 

C.6LGF-03 

C.575F-C3 

-L . 1 L Oc —02 

0.1 33c- C2 

U.8U lt-Jo 

0. 3U0E- Ji 

18 .1 1 

0.638 

0.410 

4 

C.7CCE-0 3 

0 .o? 3E-C * 

-C . 14)06-02 

0.1 14E-C2 

C.79oE-06 

0. 5006-01 

20.64 

O.o34 

0.418 

4 

C .HCCF-C3 

C. 7756-03 

-C. 10JE-J2 

3.9SoE-Ci 

0. 797E-06 

J.500E- 31 

22 .67 

0.628 

0.424 

4 

4 .SCLE-O 1 

4 • o 3 5t— Ci 

-C. I00fc-0* 

w.oopE-CJ 

D.797L-UO 

U.500E-01 

24.3 1 

3.621 

0.430 

4 

c.icr.F-02 
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J .61 3 

0.434 

4 * 

L. 1 1UF-02 
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-c. icce-02 

0.3 <oE“ C 3 
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U. jU oE- 0 1 

<6 .4 l 
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0.938 

4* 

f .I2CE-U2 

4). 1 17E-U2 

-4). iCOc -02 
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C .dOOL-Ob 

C • 5G0E- 0 l 

26.7 2 

0.605 
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4 ♦ 
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0. ) 28E-G2 

-C. 1 LUC-02 
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C . dOiE-Oo 
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2 I .0 2 
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0. 440 

4 T 

C . 1 4GF- )P 
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-C . i LUfc— 02 

6. 37*,c-03 

0.6O3C-0O 

0 . 5o 0 E- 0 1 

27 .40 

0.597 

0.442 

/ T 

C. I50F-.V 

0. 14/P-JP 

-C. 10 06- J2 

O.i )ufc- 1/3 

u.8U4t-06 

0.3006-01 

27 .76 

0.593 

0.44 3 

4 T 

f . 1 6CH- )2 
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2o .0 d 
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0. 1168*04 


0. 207E- 01 

C .35GE— Q 1 

C.690E-C3 

0.515 

581.9 

0.663E-02 


42.06 

71.82 

0 • 116E+Q4 


0. 206E- 01 

C.4GQF-01 

0.665E-03 

G.565 

639.6 

0.607E-02 


38.06 

65.95 

0 • 1 17E+Q4 


0. 205E- 01 

t.<*50F-01 

G.fc.«0E-U3 

0.636 

730.5 

0.535E-02 


32.86 

55.59 

0.119E+04 


0.202E-01 

C .500E— 0 1 

0.675E-03 

1.000 

J.106E+04 

0-371E-02 


20.64 

6.667 

0.128E+04 


0. 187E-01 

HHF/M * 0.80GE-03 INCHES 

HI - 

0.825E-03 INCHES 

DA/OR 


= 0.156E-01 

ALPHA 

=* -0. IQ0E-02 

RADIANS 



H2 = 

0.J75E-03 INCHES 

DA/ A (MIN) 

* 0.485E-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

OENSITY 

PRESSURE TEMPERATURE 

RE ( P ) 


FRICTION 

INCHES 

INCHES 

NUMBER 

HT/SEC 

L B-SEC2/FT4 


PSIA 

DEG F 



FACTOR 

C 

0.8256-03 

C.4C1 

457.5 

C.900E-02 


58.19 

82.58 

0 . 146E+D4 


0« 164E- 01 

C.500E-02 

C.B20E-03 

0.415 

473.0 

0.874E-02 


56.42 

81.38 

0 • 146 £ *04 


0.164E-01 

C.l OGF— 0 1 

C« 8 1 5E—C3 

C.430 

490.3 

0.847E-02 


54.53 

79.99 

0. 146 E*04 


0. 164E-C1 

C.l 50E-GI 

C.810E-03 

0.448 

5 10.1 

0.818E-02 


52.50 

78.34 

0.146E*04 


0. 164E-01 

C.200E-01 

C.805E-Q3 

0.469 

532.9 

C.787E-02 


50.30 

76.36 

0. 147E*04 


0. 1646-01 

C.2S0E-01 

C.800E-03 

C.494 

559.8 

0. 752E-02 


47.89 

73.92 

0 • 147E*04 


0. 163E-01 

C-3G0E-01 

0.7956-03 

C.524 

592.3 

0.71 5E-02 


45.21 

70.80 

0 • 147E*04 


0. 163E-01 

C.350E- Oi 

0.790E-03 

0.563 

633.3 

0.67 2E-02 


42. 16 

66 .62 

0.148E*04 


0. 162E-C1 

C.400F-01 

0.785E-03 

0.615 

088.4 

0.62 IE— 02 


38.52 

60.56 

0 • 149E*04 


0- 16 IE- 01 

C.450E-01 

0. 780E-03 

0.698 

772.9 

0.556E-02 


33.80 

50.29 

0.151E+04 


0.159E-C1 

C.5QCE-01 

0. 775E-03 

1.000 

0.106E+04 

0.408E-02 


22.67 

6.667 

0. 162E*04 


0. 148E-01 

hi MF AM * 0.90GF-03 INCHES 

HI = 

C.S25E-03 INCHES 

DA/OR 


= 0.1 50E-01 

ALPHA 

= -0. 100E-02 

RADIANS 



H2 - 

0.875E-03 INCHES 

DA/AIMIN) 

= 0.412E-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE ( P ) 


friction 

INCHES 

INCHES 

NUMBER 

FT/SEC 

L 8-SEC2/FT 4 


PSIA 

DEG F 



FACTOR 

0 

C.925E-03 

0.441 

5C2.2 

0.885E-02 


56-87 

79.01 

0 • 178E*04 


Q.135E-01 

Q.5Q0E-02 

C.920E-U3 

0.455 

517.8 

0.862E-02 


55.23 

77.69 

0. 1788*04 


0. 135E- 01 

C.IOQE-Ol 

0.915E-G3 

0.471 

535.2 

0.837E-02 


53.49 

76*16 

0. 178E*Q4 


0. 135E-01 

C.l 5CE-01 

C.910E-03 

0.490 

554.9 

C. 8106-02 


51.62 

74.37 

0.178E+04 


0.135E-01 

C.200E-01 

0- 9056-03 

0.511 

577.5 

0.782E-02 


49.60 

72.25 

0. 1798*04 


0.134E-01 

C .2 5CE- 0 1 

C.900E-03 

0.535 

603.9 

Q.731E-02 


47.39 

69.65 

0.1796*04 


0.134E-01 

C.30CE-01 

C •895E-03 

0.565 

035.5 

0.716E-02 


44.94 

66.39 

0 • 180E*04 


0. 1346-01 

C.350E-01 

C.890E-03 

0.603 

674.9 

0.67 7E-02 


42.14 

62.09 

0 . 180E*04 


0. 133E-C1 

C.40CE-01 

C.885E-03 

0.653 

727.1 

0.63 IE-02 


38.82 

56.00 

0* 182E*04 


0. 132E-G1 

C.450E-O1 

C.68GE-03 

0.731 

805.6 

0.572E-02 


34.50 

46.00 

0.1846+04 


0. 130E- 01 

C.50CE-0A 

0.875E-03 

1.000 

0.106E+G 

0.4376-02 


24.31 

6 . 66 7 

0.196E+04 


0.122E-01 

hIMEAM * 0. 

100E-G2 INCHES 

HI = 

0. 103E-02 INCHES 

DA/DR 


= 0.144E-01 

ALPHA 

- -0* 100E-02 

RAO IANS 



H2 * 

G.SJ5E-03 INCHES 

JA/A( M IN) 

= 0. 354E-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE ( P ) 


FRICTION 

INCHES 

INCHES 

number 

FT/SEC 

L8-S6C2/FT4 


PSIA 

DEG F 



FACTOR 

G 

0. 103E-02 

0.477 

541 .0 

0.8716-02 


55.63 

75.65 

0.210E+04 


0.114E-01 

0.5Q0E-Q2 

0. L02E-02 

0.491 

556.5 

0.8506-02 


54. 12 

74.23 

0.210E+04 


0.114E-C1 

C • 1 OQE— 0 1 

0. 1 0 IE-02 

0.507 

573.7 

0.827E-02 


52.51 

72.60 

0.210E+04 


0. 114E-01 

C- 15GE- L)1 

C. 101E-G2 

0.525 

593.2 

0.803E-02 


50.78 

70.72 

0.211E+04 


0.114E-C1 

C.200E- Jl 

0. 100E-02 

0. 546 

615.3 

0.7776-02 


48.92 

68.50 

0.211E+04 


0.I14E-C1 

C.25GE-01 

C. 100E-02 

0.570 

640.9 

0.7486-02 


46.89 

65-81 

0.21 1E+04 


0. 113E- Cl 

C.4O0F-JI 

C.995E-03 

0.59 9 

671.5 

C.716E-02 


44.62 

62.48 

0.2126*04 


0.113E-01 

C.35QE-01 

0.990E-03 

0.636 

709.2 

0.6816-02 


42.05 

58.15 

0.213E+04 


0. 113E-Q1 

C.*OOF- 01 

0.985E-C3 

C.684 

758.5 

0.639E-02 


39.00 

52.12 

0.2156*04 


0. 112E-01 

C.45CF— 01 

0.980E-03 

0.757 

831 .5 

0.585E-02 


35.03 

42.46 

0 .2186*04 


0. HOE-01 

C.500F-,)1 

0.975E-03 

1.000 

0.1 06E+ C4 

0.4o 16-02 


25.64 

6.667 

0.23 1E+04 


0. 104E-01 

hi HE AM = 0. 

110E-02 INCHES 

HI - 

0.112E-02 INCHES 

DA/DR 


= 0.138E-01 

ALPHA 

= -0.100E-02 

RADIANS 



H2 = 

0. 107E-02 INCHES 

DA/A ( MIN) 

= 0.307E-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE ( P) 


FRICTION 

INCHES 

INCHES 

NUMBER 

f T/SEC 

L 8-SEC2/FT 4 


PSIA 

DEG F 



FACTOR 

C 

0. 1 12E-02 

C. 500 

565.8 

0.862E-J2 


54.81 

73.36 

0 . 2396*04 


0. 104E-01 

C «5G0E— 02 

C.112E-02 

0.514 

580.9 

0.842E-02 


53.39 

71.92 

0.2406+04 


0. 104E-01 

t« IGuF-Ol 

0.1 1 1 E— 02 

U.52S 

597.7 

C.821E-02 


51.88 

70.27 

0.240E+04 


0. 1055- Cl 

C.150F-01 

C.l 1 IE-02 

0.547 

616.5 

J. 7986-02 


50.26 

68.37 

0.240E+04 


0. 1055-01 

C .2U0F-QI 

C.l 10E-U2 

0 . 567 

6 J 7 .9 

0«77 4E-02 


48.52 

66.14 

0 .2406+04 


0. 1055-01 

C.25GE-01 

C. 1 10E-02 

C.591 

662.0 

0 • 747E-02 


46 • 62 

63.46 

0.241 E+04 


0. 105E-C1 

C.3GCF-GI 

0. 109E-02 

0.619 

691.8 

0*7 l6c-0 2 


44.50 

6 0.17 

0 .242 E ♦ 04 


0 » 1 05E- 01 

C.350E-01 

C. IC9E-02 

C.C54 

?2 7 • 7 

0.08 4E— 0^ 


42.09 

55.93 

0 . 243 E + 04 


0* 105E-01 

i .400F-»M 

C. 1C8E-02 

C. 7U0 

1 74.5 

0 • 645l-02 


->9.22 

50.0 7 

0 • 245 E + 04 


0. 105E-C1 

C.45GF-01 

0. 1G8F-U2 

0 • 765 

84J.5 

J. 5946-02 


35.46 

40.79 

). 248604 


3*1356-31 

C.*>OOF- 11 

C. IC7E-02 

1.000 

0.1 Uo£* C4 

0*4 756-02 


26.41 

6.66 7 

0.2626*04 


0. 1055- Cl 


67 


MMFAM = J.: 

L2GE-02 INCHES 

HI = 

0. 1*26-02 INCHES 

OA/OR 


= 0. 1516-01 

ALPHA 

= -0. 1006-02 RADIANS 



H2 = 

0.117E-O2 INCHES 

0 A/ A 1 H IN ) 

= o. 2686-01 


ICHOKED FLOW) 


K 

H 

MACH 

VELOCITY 

OEMS ITY 

PRESSURE TEMPERATURE 

RE l P ) 

friction 

1 NCHFS 

INCHES 

NUMBER 

FT/SEC 

L B-SEC2/FT 4 


PSIA 

UEG F 


FACTOR 

f. 

0. 1226-02 

0.511 

57 7.8 

0.857E— 02 


54.43 

72.22 

0.265E*04 

0. 106E-CI 

C.5CGF-J2 

0. 1226-02 

0.525 

592.0 

0.838E-02 


53.01 

70.75 

0.2656*04 

0. 106E-C1 

c. iooe-oi 

0. 121E-Q2 

0.540 

609.4 

C.817E-02 


51.53 

69.09 

0.266E+04 

0. 10o6- Oi 

C. 1506-01 

C. 1216-02 

0 .558 

o2d *0 

0.7956-02 


49.95 

67.18 

3.2666*04 

0. 106E-01 

C.200F-U1 

C. 1206-U2 

O.57o 

649.1 

0.77 IE-02 


48.25 

64.9 3 

0.266E*04 

0. 1066-01 

c.25ce-oi 

C. 1206-02 

0.601 

673.4 

0. 7456-02 


46.39 

o2.26 

0.267E*04 

0. 106E-01 

C .300F-01 

0.1196-02 

0.629 

702.2 

0.717E-02 


44.33 

58.96 

0.2686*04 

0. 106E-C1 

C. 35GF-01 

0.1 19E-02 

0.663 

737.5 

0. 6846-02 


41.98 

54.74 

0.2696*04 

0. 106E-C1 

L .40CF-UI 

0. 1 186-U2 

0. 7C8 

763.3 

0.646E-02 


39.18 

48.94 

0.2716*04 

0. 106E-CI 

C. 4506-01 

C. I 186-02 

0.776 

B50.6 

0.5966-02 


35.53 

39.79 

0.2 75 £♦ 04 

3. 106E- 31 

C .5006-01 

C.l 1 76-02 

1.000 

0 . 106E+ 04 

0. 4B0E-02 


26*72 

6.667 

0.2906*04 

0.107E-01 

H ( ME AM « 0. 

I JOE-02 INCHES 

HI - 

0.1336-02 INCHES 

OA/OK 


* 0. 125E-3 1 

ALPHA 

* >0.100E-02 RADIANS 



H2 * 

0.1286-02 INCHES 

OA/AI MINI 

* 0.235E-01 


(CHOKED FLOW 1 


X 

H 

MACH 

VELOCITY 

DENSITY 

PRESSURE TEMPERATURE 

RE I P ) 

FRICTION 

INCHES 

INCHES 

NUMBER 

FT/SEC 

L B-SEC2/FT4 


PSIA 

OEG F 


FACTOR 

C 

C. 1336-02 

0.522 

589.2 

0.6536-02 


54.00 

71.11 

0.2916*04 

0. 1076-01 

C.500F-G2 

0. 132E-G2 

0.535 

604.1 

0.8346-02 


52.64 

69.63 

Q.292E*04 

3. 1376- Cl 

C. 1006-01 

0.1316-02 

0.551 

620.6 

0.8136-02 


51.19 

67.95 

0 .2926*04 

0. 107E-01 

C. 1506-01 

0.1316-02 

0.568 

639.0 

0.7926-02 


49.65 

66.02 

0.2926*04 

0.1076- Cl 

C. 2006-01 

0. 1 306—02 

0.588 

6 59.9 

0. 7696-02 


47.98 

63-76 

0.2 93 £*04 

0.1376-01 

0. 2506-01 

C. 1306-02 

0.611 

683.9 

0.743E-02 


46.17 

61.07 

0.2946*04 

0. 107E-01 

C.iOOF-Gl 

0.129F-Q2 

0.639 

712.3 

0.7156-02 


44. 15 

57.78 

0 .2956*04 

0. 107E-C1 

C . 3506-01 

0.1296-02 

0.b72 

747.0 

0.684E-02 


41.85 

53.56 

0.2966*04 

0. 1376-31 

C .40CF— 01 

0. 1286— 02 

0.717 

792.0 

0. 6466-02 


39.12 

47.80 

0 .2986*04 

0. 1076-Ci 

C. 4506-01 

0. 1286-02 

0.784 

657.9 

C .5986-02 


35.56 

38.75 

0.302E*04 

0.1076-01 

C. 5006-01 

0. 1276-02 

1.000 

0. I0o£*04 

0.4866-02 


27.02 

6.667 

0.318E*04 

3.105E-31 

MM64N) « 0. 

1406-02 INCHES 

HI* 

0.1426-02 INCHES 

OA/OR 


* 0.1196-01 

ALPHA 

* -0. 100E-02 RADIANS 



H2 * 

0.1376-02 INCHES 

OA/AI MIN) 

* 0-2076-01 


(CHOKEO FLOW) 


X 

H 

MACH 

VELCCITY 

OENSITY 

PRESSURE TEMPERATURE 

RE I P ) 

FRICTION 

INCHES 

INCHES 

NUMBER 

FT/SEC 

L B-SEC2/FT4 


PSIA 

DEG F 


FACTOR 

C 

0.1426-02 

0.534 

602.9 

0.847E-02 


53.52 

69.75 

0.319E*04 

0. 105E-01 

C. 5006-02 

0.142E-U2 

0.548 

617.7 

0.8296-02 


52.19 

68.24 

0.3206*04 

0. 105E-01 

C.100F-01 

0.14 IE— 02 

0.564 

6 34.1 

0.8096-02 


50.78 

66.54 

0.3206*04 

0. 105E-01 

C. 1506-01 

C.141E-02 

0.581 

652.3 

0.788E-02 


49.28 

64.59 

0.3206*04 

0. 105E-01 

C.200E-01 

Cm 1406-02 

0.601 

672.9 

0.766E-02 


47.66 

62.32 

0.321E*04 

3.135E-01 

C. 2506-01 

C. 1406-02 

0.623 

696 .6 

0. 7416-02 


45.90 

59.62 

0.3226*04 

0.105E-C1 

Cm 3G0F-O1 

C. 1396-02 

0.650 

724.4 

0.7146-02 


43.94 

56.32 

0*3236*04 

0. 105E-C1 

C. 3506-01 

C. 1396-02 

0.684 

758.4 

0.6836-02 


41.72 

52.13 

0.3246*34 

3.135E-31 

C. 4006-01 

0.1386-02 

0.727 

802.3 

Q. 64 76-02 


39.08 

46.42 

0.327E*04 

0. 104E-01 

C.450F- 0 1 

0.138E-02 

0.792 

866 *2 

0.6016-02 


35.64 

37.56 

0.331E*04 

0. 1046-01 

0.50 OF- 01 

0.1376-02 

1.000 

0.1 06E+C4 

0.493E-02 


27.40 

6.667 

0.3476*04 

0. 133E-01 

MMFAM * ). 

1506-02 INCHES 

HI * 

0.1526-02 INCHES 

DA/Ok 


- 0. 1126-01 

ALPHA 

* -0.1006-02 RADIANS 



H2 * 

0.1476-02 INCHES 

OA/AI MIN) 

* 0.1836-01 


(CHOKED FLOW) 


X 

H 

MACH 

VELOCITY 

OENSITY 

PRESSURE TEMPERATURE 

RE IP) 

FRICTION 

INCHES 

INCHES 

NUMBER 

FT/SEC 

L8-SEC2/FT4 


PSIA 

DEG F 


FACTOR 

0 

0. 1526-02 

0.547 

616 .2 

0.8426-02 


53.04 

68.40 

0.3486*04 

0. 103E-01 

C .50 0 6— 02 

C. 1526-C2 

0.561 

630.8 

0. 8246-02 


51.76 

66.88 

3.3486*04 

0.1036-01 

C.IGOF-Oi 

C • 1 5 16-02 

0.5 76 

646.9 

0.8056-02 


50.39 

65.17 

0. 348E*04 

0. 103E-C1 

C. 1506-01 

0. 1516-02 

C. 593 

664 .9 

0. 785E-02 


48.93 

63.21 

0.3496*04 

0. 1036-01 

C. 2006-01 

C. 1506-02 

0.612 

685.2 

0.7636-02 


47.37 

60.93 

0.3496*04 

0. 1336-01 

C.250F-QI 

C. 1506-02 

0.635 

7C8.4 

0.7396-02 


45.66 

58.23 

0.3506*04 

0. 1036- Cl 

C.3006-GI 

0. 1496-02 

0.661 

735 • 7 

0.713E-02 


43.76 

54.95 

0. 352E*04 

0. 1036-01 

C.350F-01 

0.1496-02 

C.6S4 

769.0 

0.6836-02 


41.61 

50.79 

0.3536*04 

0. 1026-C1 

C. 4006-01 

0. 1486-02 

0.737 

811.7 

0*6496-02 


59.05 

45.16 

0.356E*04 

0. 102E-01 

C.45C6-01 

0. 1486-02 

C.800 

8 73.7 

0.6046-02 


35.73 

36.48 

0.3606*04 

3. 132E-31 

C. 5006-01 

0.1476-02 

1.000 

0. 106 Eh 04 

0.499t— 02 


27.76 

6.66 7 

0.378E+04 

0. 1016-01 

MMFAM * 0. 

1606-02 INCHES 

Hi * 

0.1636-02 INCHES 

UA/OR 


* 0.106E-01 

ALPHA 

* -0. 1006-02 RADIANS 



H2 * 

0.1586-02 INCHES 

OA/AI MIN) 

* 0.162E-J1 


(CHOKED FLOW) 


X 

H 

MACH 

VELOC ITY 

DENSITY 

PRESSURE TEMPERATURE 

RE ( P ) 

FRICTION 

INCHES 

INCHES 

NUMBER 

FT/SEC 

L0-SEC2/FT4 


PSIA 

DEG F 


FACTOR 

C 

0. 1 6 3 E — 0 2 

0.558 

628.5 

C. 8376-02 


32.00 

67.12 

0.3766*04 

0. 101E-C1 

C.5GCF-02 

0. 162E-02 

u.572 

8*»3 .0 

0. 8206-02 


51.35 

65.60 

0.3776*04 

0. 1016-01 

C. tuOE-tf l 

C. 1 6 IE — U2 

0.587 

658.9 

C.aOiE -02 


50.02 

63.87 

0.3776*04 

0.1316-01 

C. 1506- Ji 

0.1 6 IE— 02 

0.6 04 

676.5 

0.781E-02 


48.61 

61.91 

0.378E*04 

0. 10 IE- Cl 

C.20CF-01 

0. 160E-02 

0.623 

696 .5 

0. 7606-02 


4 7.08 

59.63 

0.378E*04 

0. 101E-C1 

C. 2506-01 

C . 1 <>06-02 

0.645 

714.3 

0.73/6-02 


45.43 

56.94 

0.3796*04 

0. 101E-G1 

C . iOCE- ■) l 

0. Ir> 96-02 

0.672 

?«*6. 1 

0.712E-02 


43.59 

53.67 

0.38 1 E*04 

0. 10 IE - 01 

C. x506~0l 

C. 159F-02 

C. 704 

7 78.6 

0.6856-02 


41.50 

49.55 

0.3826*04 

0. IOOE-OI 

C. 4006-01 

C. 158E-02 

0. 7*5 

820.3 

0.6506-02 


39.03 

44.00 

0.385E*04 

0. 1006- 01 

t 5C.F- U 1 

C. 1 586 -02 

0.607 

680 .4 

0.6 J6t-02 


35.81 

35.49 

3. 390 E* 04 

0. 1006 -Cl 

L.5GCF-01 

0. 1 *> 76-02 

l .OOU 

0 . 1 QcjE ♦ 04 

0. 5056-02 


28.08 

6.66 7 

0.408E *04 

O.908E-C2 
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TABLE I. - PARAMETERS CALCULATED BY AREAX 


Parameter 

Formula 3 

Units 



SI 

U.S. 

Gas constant 

^ Universal gas constant 

JAg-K 

lbf-ft 

Molecular weight 

lbm-°R 

Seal surface area 

A=»(r 2 2 -r2) 

m 2 

in. 2 

Flow length 

AR = Rg - R 1 

m 

in. 

Flow width (mean) 

w - 2 ’<rrj 

m 

in. 

Mass flow rate 

M = Whpu 

kg /sec 

lbm/min 

Volume flow rate 

Q = CjM 

scms 

scfm 

Reynolds number due 
to rotational flow 

b Re(r) = pVh ^/jT 



Knudsen number 

k„, 2 ' 96M »» 

Re(p) 



Mean free path 

A = Kn x h 

m 

in. 

Viscosity of air 
(Sutherland's law) 

CgT 1 ' 5 
T + c 3 

2 

N-sec/m 

lbf-sec/in. 2 

Power 

pAV 2 /h 

W 

hp 

Apparent temperature 
rise due to power 
dissipation 

c . x Power 
AT - ’ 

Mc p 

K 

°R 

Shear heat 

H = Cg X Power 

W 

Btu/min 

Torque 

Power /Speed 

N-m 

ft-lb 

Seal opening force 

/•AR 

F = W J (P - P 3 )dx 

N 

Ibf 

Center of pressure 

x = — / (P - P^X dx 

C F J 0 

m 

in. 

Axial film stiffness 

S = -dF/dh 

N/m 

lbf/in. 


a Constants in equations are as follows (for Si and U.S. units, respectively): 
^ = 5.051554X10' 3 (13.083); Cg = 1.4591X10" 6 (1 . 57639X10 -10 ) ; 
c 3 = 110.33333 (198.6); c 4 = 1.0 (42.42); and c g = 1.0 (42.42). 

' 3 Where p is density at midseal, JJ is viscosity at midseal, and V is 
mean rotational velocity. 




TABLE n. - VARIABLES IN NAMELIST /SEAL/ 


Variable 

name 

Description 

Units 


SI 

U.S. 

RINNER 

Inner radius of seal 

m 

in. 

ROUTER 

Outer radius of seal 

m 

in. 

RDIFIN 

Flow length 

m 

in. 

MOLWT 

Molecular weight 



CP 

Specific heat 

J/kg-K 

Btu/lbm-°R 

MUIN 

Reservoir viscosity. The pregram will calculate MUIN 
for air but not for other gases. 

o 

N-sec/m 

lbf-sec/in. ^ 

GAMMA 

Ratio of specific heats 



SPEED 

Rotational velocity 

rps 

rpm 

CAPV 

XLAM 

Seal -face speed 

Exponent in friction factor - Reynolds number relation 
for laminar flow (eq. (22)) 

m/sec 

ft/sec 



XTURB 

Exponent in friction factor - Reynolds number relation 
for turbulent flow 



CONLAM 

Constant in friction factor - Reynolds number relation 
for laminar flow 




Constant in friction factor - Reynolds number relation 
for turbulent flow 

RELAM 



Maximum Reynolds number for laminar flow 



RETURB 

Minimum Reynolds number for turbulent flow 



PWRSKP 

Logical variable. If it is set to .TRUE. , calcula- 
tions involving power are omitted. 



NRMSKP 

Logical variable. If it is set to .TRUE. , normal- 
ized values of F and x c will be omitted. 



PRSSKP 

Logical variable. If it is set to .TRUE. , printout of 
distributions across face of seal of P, T, p, u, M, 
f, and Re will be omitted. 



PLTSKP 

Array of logical variables. If any element is set to 
.TRUE. , the corresponding plot will be omitted: 
PLTSKP(l) applies to plot of power against h. 
PLTSKP(2) applies to plot of x e against h. 
PLTSKP(3) applies to plot of F against h. 
PLTSKP(4) applies to plot of P against x. 
PLTSKP(5) applies to plot of T against x. 
PLTSKP(6) applies to plot of p against x. 
PLTSKP(7) applies to plot of M against x. 
PLTSKP(8) applies to plot of Reynolds number 
against x. 

PLTSKP (9) applies to plot of mean friction 
factor against x. 



NSI 

Numerical indicator for units. NSI = 1 means SI 
units. NSI = 2 means U.S. units. 






TABLE HI. - VARIABLES IN NAMELIST /HDATA/ 


Variable 

Description 

Units 

name 


SI 

U.S. 

HMEAN 

Mean film thickness 

m 

in. 

ALPHA 

Tilt angle 

rad 

rad 

NJ 

Number of film thicknesses 


— 

JDONE 

Number of film thicknesses for which 
cards were punched in previous 
running of program 

— — ” 



TABLE IV. - VARIABLES IN NAMELIST/RESDAT/ 


Variable 

Description 

Units 

name 


SI 

U.S. 

piin 

Gas pressure at inner radius 

N/m 2 

lbf/in. 2 

P2IN 

Gas pressure at outer radius 

N/m 2 

lbf/in. 2 

PRIN 

Ratio of sealed-gas pressure to ambient pressure, 

V P 3 



TOIN 

Sealed-gas temperature (upstream reservoir temper- 
ature), T Q 

K 

°F 

LOSS 

Entrance velocity loss coefficient 



INCODE 

Input code. For running many cases with one 
loading of the program, the input code tells 
the program what new data are expected for 
the next case: 

INCODE = 1 means a new title card is expected. 

INCODE = 2 means new SEAL data are expected. 

INCODE = 3 means new HDATA data are expected. 

INCODE = 4 means new RESDAT data are expected. 






TABLE V. - INPUT DATA FOR SAMPLE PROBLEM 


2-6 


7-12 


$seal| 

GAMP 

CONTRl 

pltsk| 

$HDA 
2 . 
ALPH 
tRES 
INCC 


. 54 * 


DiAI 


ISAMPLE 
RINNEp 
1 . 4 , 
|B=. 079 j 
= 9 *F , 
HME 
- 5 , 2 . 
A|= 20 * 

T Pll 
= 1 * 




EAN 


13-18 


PROBL) 
= .08 
!muin=g| 
REL 
N S I = l| 
= .76211 
) 794 E 
bO l, 

N = 4481 | 




AW 


19-24 


Ah 


JEM - 
31 , 

|., SPEtE 
= 2300 | 
WI 
lE- 5 , 1.1 
3 . 048 |l 
|J= 14 , 
59 . 22 , 


MCI 


idtC 


25-30 


EA EX 
UTER = 
D= 0 . , 
., RET 
=0. i 
| 0 16 E -5 
E- 5,3 

Ijdone 

P 2 1 N=| 


□N 


PANSI 
1 . 08420)1 
CAPV 
|jRB= 30 (l 


31-36 




PROG|R 
RD 
60 . 96 , 

00., P* 


,1.270(1 

I302E- 

14 % 
1034231 


E- 5 , 1 . 

3.55 


37-42 


p 24 E- 5 | 

- 5,3 




. 59 , P 


43-48 

AM - 

I N= 0 . 
XLAM = 
IRSKP = 


RIN=0. 


49-54 




I UN 
MOLHlT 
\L., XT 

f. prs 


s 

= 28 . 9 | 
|URB = . 
SKP = F , 


,1.778 

B 10 E -5 


P- 5 , 2 . 
, 4 . 


55-60 


P> 6 , CP 
cun! 

NRMSK 


2P 


064E 


p 32 C- 5 | 
-5, 


TOI N(= 31 1 • 1 


61-66 


1004 
PAM= 24 | 
= F, 


6*C 


, 2 .2 8 < 

0 . , 


MU. 


67-72 


Loss=i 


73-80 


16 , 

• , 


6fe-5, 


SA 

$SE ALj 
MU I N 
| CONTRjl 

pltsk| 

$HDA 
1 . 3 E 
$RE S 


MPLE Pf| 
RINN 
.. G 
|8 = . 079 j 
]p= 9 *T,' 

hmeaJi 
| 3 , 1 . 4 Ef 
PI 


m 


OAT 


lOBLEM 
= 3 . 2 E|! 
A|MMA= 1 •[ 
REL 
NS I = 
|N=. 3 E 
3,1 
= 65 ., 


3 R = 


IC= 


1 - AREA 
5 , ROL 
K» CAP 
2300 

$ 

3 , . 4 E — | 
3,1 
P 2 1 N=| 


AM 


9E 


XPAN 
|T E R = 3 

V= 2 C 0 .| 
RE 


|S I ON 
315 , 
SPEE) 
T|UR 6=3 


PROGRAM 
IFIN 
0 = 0 ., 


FID I 


qoo. 


- u. 

0., Ml 
XLAM = 1 
WRSXP=) 


S. UN Ij 

(dth=o 

XT 

NRH 


TS 

MGLWl 
8 =. 25 j 
SKP=F , 


1*1 


|T =28 
CUNLl 
PRSSK 


66 , CPI 
|AM = 24 .| 
|P = F, 


= .? 4 , 


3 , . 5 E ■ 
f>E- 3 , 
115 . , 


, 3 ,. 7 E -1 
PHA= 
PhIN=OJ, TOI 


|3 , . 6 E 

0 • t 




|3, . 8E- 

*-.oo| 

100 . 


2to 


Nh 


, 3 ,. 9 E 
1 , NJ 
LOSS 


3 . 1 . 0 E 
14 , JO 
1 ., IN 


:]- 3 , 1 . l|E - 3 » 1 • 
'DNE = 14 i $ 

1 C ODE = ll * 


2E-3 , 




TABLE VI. - PARAMETERS IN COMMON BLOCK/ARRAYS/ 


Array 

Symbol 

Array 

dimension 

Description 

1 

X 

(ID 

Distance across face of seal 

2 

p 

(11,20) 

Pressure 

3 

M 

(11,20) 

Mach number 

4 

u 

(11,20) 

Leakage flow velocity (x -direction) 

5 

T 

(11,20) 

Temperature 

6 

P 

(11,20) 

Density 

7 

Re(p) 

(11,20) 

Pressure-flow Reynolds number 

8 

7 

(11,20) 

Mean friction factor 

9 

h 

(11,20) 

Film thickness 

10 

x c 

(20) 

Dimensionless center of pressure 

11 


(20) 

Power 

12 

F 

(20) 

Force 

13 

S 

(20) 

Axial film stiffness 

14 

F 

(20) 

Pressure profile factor 

15 

h l 

(20) 

Film thickness at entrance 

16 

h 2 

(20) 

Film thickness at exit 

17 

a 

(20) 

Tilt angle 

18 

h 

m 

(20) 

Mean film thickness 

19 

X* 

(20) 

Choking flow length 

20 

p* 

(20) 

Choking pressure 

21 

A* 

(20) 

Area at point of choking 



TABLE Vn. - PARAMETERS IN COMMON BLOCK/CONSTS/ 


Word 

Symbol 

Description 

1 

y 

Ratio of specific heats 

2 

AR 

Flow length 

3 

R 0 

SIGN 

Radius of seal at entrance 

If SIGN > 0, flow is radially outward; if SIGN < 0, 
flow is radially inward. 

4 

n l 

Exponent in friction factor - Reynolds number relation 
for laminar flow 

5 

n t 

Exponent in friction factor - Reynolds number relation 
for turbulent flow 

6 

k l 

Constant in friction factor - Reynolds number relation 
for laminar flow 

7 

k t 

Constant in friction factor - Reynolds number relation 
for turbulent flow 

8 

(Re), 

Upper limit on Re for laminar flow 

9 

(Re) t 

Lower limit on Re for turbulent flow 

10 

T 0 

Sealed-gas-reservoir temperature 

11 

P 0 

Sealed-gas-reservoir pressure 

12 

P 3 

Ambient pressure 

13 

P tol 

Pressure tolerance to stop iteration on exit pressure 
for subcritical flow 

14 

3t 

Gas constant 

15 

C L 

Entrance velocity loss coefficient 

16 

^0 

Sealed-gas-reservoir viscosity 

17 

77 

3.1415927 

18 

« u (2) 

Universal gas constant 

19 

CNVT(11, 2] 

Constants needed for internal calculations 

20 

NSI 

Numerical indicator for system of units being used: 
NSI = 1 means SI units; NSI = 2 means U.S. units. 




TABLE Vm. - CONSTANTS FOR INTERNAL CALCULATIONS 


Word 
in array 
CNVT 

Variable to which 
constant applies 

1 

Sutherland's law 

2 

Sutherland’s law 

3 

Speed 

4 

Mass flow rate 

5 

Standard volume flow 

6 

Reynolds number 


Word 
in array 
CNVT 

Variable to which 
constant applies 

7 

Power 

8 

Shear heat 

9 

Torque 

10 

Density 

11 

Velocity 


TABLE IX. - PARAMETERS IN COMMON BLOCK /TRAYS/ 


Array 

Symbol 

Array 

dimension 

Description 

1 

N 

(1) 

Number of grid points 

2 

X 

(101) 

x-grid points 

3 

M 

(101) 

Mach number at each x 

4 

h 

(101) 

Film thickness at each x 

5 

NJ 

(1) 

Index of film thickness for which 
solution is being found 

6 

r* 

(1) 

Temperature at choking 



M M + dM 

Figure I, - Model and notation of sealing faces, including control volume for quasi-one-dimensional flow 
with area change. 










Figure 4. - Heirachy of subprogram calls. 
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meters (3. 315 in. ). 



Mean film Flow regime 
thickness, 

^nv 

pm (mil) 

o 2.5 10,1) Subcritical 

□ 5.1 (.2) Critical 

& 7.6 (.3) Supercritical 

d 12.7 (.5) Supercritical 

o 12.7 (.5) Supercritical - Parallel-film, 

variable-area analysis 

2.5 (.1) Subcritical 1 Parallel-film, 

12.7 (.5) Supercritical] constant-area analysis 



U 40 80 120 160 200 240 

Pressure, P, psia 


Figure 8. - Results obtained using variable-area approximate 
analysis for pressure distributions. Positive 2-milliradian 
tilt; conditions represent subcritical, critical, and super- 
critical flow; mean film thicknesses, 2.5, 5.1, 7.6, and 
12.7 micrometers (0.1, 0.2, 0.3, and 0.5 mil); sealed-gas- 
reservoir pressure, 148 N/cm 2 (215 psia); exit pressure, 

10. 3 N/cm‘ (15 psia); sealed-gas-reservoir temperature, 

700 K (800° F); inner radius, 8. 300 centimeters (3. 265 in. ); 
outer radius, 8. 410 centimeters (3. 315 in. ). Also shown 
are comparable parallel-film cases. 
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